Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,743,306 C→T 5.4% Q192* (CAG→TAG)  yiaK → 2,3‑diketo‑L‑gulonate reductase, NADH‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,743,3060CT 5.4% 237.3 / 2.3 74Q192* (CAG→TAG) yiaK2,3‑diketo‑L‑gulonate reductase, NADH‑dependent
Reads supporting (aligned to +/- strand):  ref base C (31/39);  new base T (2/2);  total (33/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

TGCCGCCGTGGGGCGCAAAAGAGTGTCGCATAGGCACTAACCCGCTGATCGTCGCCATTCCTTCCACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTC  >  NC_000913/3743163‑3743454
                                                                                                                                               |                                                                                                                                                    
tGCCGCCGTGGGGCGCAAAAGAGTGTCGCATAGGCACTAACCCGCTGATCGTCGCCATTCCTTCCACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTc                                                                                                                                                 >  7:53740/1‑149 (MQ=255)
  ccgccgTGGGGCGCAAAAGAGTGTCGCATAGGCACTAACCCGCTGATCGTCGCCATTCCTTCCACGCCGATCCCCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTccc                                                                                                                                               >  2:22784/1‑149 (MQ=255)
  ccgccgTGGGGCGCAAAAGAGTGTCGCATAGGCACTAACCCGCTGATCGTCGCCATTCCTTCCACGCCGATCACCATGGTCGATATGGCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTcc                                                                                                                                                <  2:286483/148‑1 (MQ=255)
    tccgTGGGGCGCAAAAGAGTGGCGCATAGGCACTAACCCGCTGATCGTCGCCACTCCTTCCAAGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCgg                                                                                                                                             <  6:214016/148‑1 (MQ=255)
      cgTGGGGCGCAAAAGAGTGTCGCATAGGCACTAACCCGCTGATCGTCGCCATTCCTTCCACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTc                                                                                                                                           >  5:218383/1‑149 (MQ=255)
          gggCGCAAAAGAGTGTCGCATAGGCACTAACCCGCTGATCGTCGCCATTCCTTCCACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGAtg                                                                                                                                       <  8:52302/149‑1 (MQ=255)
          gggCGCAAAAGAGTGTCGCATAGGCACTAACCCGCTGATCGTCGCCATTCCTTCCACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGAt                                                                                                                                        <  8:90155/148‑1 (MQ=255)
             cgcAAAAGAGTGTCGCATAGGCACTAACCCGCTGATCGTCGCCATTCCTTCCACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGAtggtg                                                                                                                                    <  5:77701/149‑1 (MQ=255)
             cgcAAAAGAGTGTCGCATAGGCACTAACCCGCTGATCGTCGCCATTCCTTCCACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGAtggtg                                                                                                                                    <  1:47074/149‑1 (MQ=255)
                aaaaGAGTGTCGCATAGGCACTAACCCGCTGATCGTCGCCATTCCTTCCACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCt                                                                                                                                 >  1:155161/1‑149 (MQ=255)
                         tCGCATAGGCACTAACCCGCTGATCGTCGCCATTCCTTCCACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTtgatgat                                                                                                                         >  1:108718/1‑148 (MQ=255)
                         tCGCATAGGCACTAACCCGCTGATCGTCGCCATTCCTTCCACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTtgatgat                                                                                                                         >  3:157534/1‑148 (MQ=255)
                           tcATAGGCACTAACCCGCTGATCGTCGCCATTCCTTCCACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTCGCAGGTCGTCGGCTCCCGGTCGATGGTGGCTTTGATGATGAg                                                                                                                      >  7:243685/2‑149 (MQ=255)
                           gCATAGGCACTAACCCGCTGATCGTCGCCATTCCTTCCACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAg                                                                                                                      <  1:193774/149‑1 (MQ=255)
                           gCATAGGCACTAACCCGCTGATCGTCGCCATTCCTTCCACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGGGGCTTtgatgatga                                                                                                                       <  4:325272/148‑1 (MQ=255)
                            cATAGGCACTAACCCGCTGATCGTCGCCATTCCTTCCACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGAGGATGAgg                                                                                                                     >  4:213482/1‑149 (MQ=255)
                             aTAGGCACTAACCCGCTGATCGTCGCCATTCCTTCCACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAggg                                                                                                                    <  1:60014/149‑1 (MQ=255)
                                   aCTAACCCGCTGATCGTCGCCATTCCTTCCACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAAttt                                                                                                              >  3:257916/1‑149 (MQ=255)
                                     tAACCCGCTGATCGTCGCCATTCCTTCCACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTtgatgatga                                                                                                                       >  3:258462/1‑138 (MQ=255)
                                     tAACCCGCTGATCGTCGCCATTCCTTCCACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTtgatgatga                                                                                                                       <  4:258462/138‑1 (MQ=255)
                                          cGCTGATCGTCGCCATTCCTTCCACGCCGATCACCATGGACGATATGTCGATGTCGATGTTCTCTTACGGCCTGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCGCGGTCGATGGTGGCTATGATGAGGAGGGCAATTTGACCa                                                                                                         >  3:300121/1‑147 (MQ=255)
                                                      ccATTCCTTCCACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGt                                                                                            >  2:89378/1‑148 (MQ=255)
                                                        aTTCCTTCCACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTAt                                                                                         <  5:15966/149‑1 (MQ=255)
                                                        aTTCCTTCCACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAAATTGACCAAAGAACCTGGCGTTAt                                                                                         >  7:180685/1‑149 (MQ=255)
                                                           ccttccACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATc                                                                                        >  7:111163/1‑147 (MQ=255)
                                                            cttccACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCgag                                                                                     <  3:330336/149‑1 (MQ=255)
                                                             ttccACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCgag                                                                                     >  5:324458/1‑148 (MQ=255)
                                                                cACGCCGATCACCATTGTCGATATGTCGATGTCTATGTTCTCTTACGTCATGTTAGAAGTTAAGCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGATGGCAATTTGACCAAAGAACCTGGCGTTATCGagaag                                                                                  <  8:253224/148‑1 (MQ=255)
                                                                cACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTTGCTTTGATGATGAGGGCAATTTGACCAAAGCACCTGGCGTGATCGagaag                                                                                  >  3:13668/1‑148 (MQ=255)
                                                                   gCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATc                                                                              <  2:122824/149‑1 (MQ=255)
                                                                   gCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTGACGGCATGTTAGAAGTTAACCGTCTGGCAGGTAGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGAGTTATCGAGAAGAATc                                                                              >  6:113705/1‑149 (MQ=255)
                                                                   gCCGATCACCATGGGCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTTTCAGTCTGGCAGGTCTTCAGCTCACGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATc                                                                              <  2:107542/149‑1 (MQ=255)
                                                                       aTCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATcgccg                                                                          >  2:76624/1‑149 (MQ=255)
                                                                         cACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCa                                                                        >  8:135367/1‑149 (MQ=255)
                                                                              tGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGGCGTCAGCTCCCGGACGATGGTGGCTTTGATGATGAGGGCAATTTGCCCAAAGAACCTGGCGTGATCGCGAAGAATCTCCGCAttttt                                                                   >  6:202788/1‑148 (MQ=255)
                                                                                   gATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGa                                                               <  2:155161/148‑1 (MQ=255)
                                                                                    atatGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATg                                                             <  4:157534/149‑1 (MQ=255)
                                                                                    atatGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATg                                                             <  2:108718/149‑1 (MQ=255)
                                                                                     tatGGCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATgg                                                            <  6:325189/149‑1 (MQ=255)
                                                                                         tcgatgtcgatgTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCg                                                                                             <  6:228291/112‑1 (MQ=255)
                                                                                         tcgatgtcgatgTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCg                                                                                             >  5:228291/1‑112 (MQ=255)
                                                                                         tcgatgtcgatgTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAATAATCGCCGCATTTTGCCGATggg                                                           <  1:22784/146‑1 (MQ=255)
                                                                                         tcgatgtcgatgTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTa                                                        >  7:169434/1‑149 (MQ=255)
                                                                                          cgatgtcgatgtTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTTAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTAc                                                       >  1:349887/1‑149 (MQ=255)
                                                                                          cgatgtcgatgtTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTTAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTAc                                                       >  1:344923/1‑149 (MQ=255)
                                                                                          cgatgtcgatgtTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTAc                                                       <  3:213482/149‑1 (MQ=255)
                                                                                          cgatgtcgatgtTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGGCAGCTCCCGGTCGATGGTGGCTTGGATGATGACGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTAc                                                       <  6:328837/149‑1 (MQ=255)
                                                                                           gatgtcgatgtTCTCTTACGGCATGTTAGAAGTTACCCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACt                                                      >  1:81718/1‑149 (MQ=255)
                                                                                             tgtcgatgtTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATc                                                                                        >  6:304314/1‑113 (MQ=255)
                                                                                             tgtcgatgtTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATc                                                                                        <  5:304314/113‑1 (MQ=255)
                                                                                                 gatgtTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTGGATGGTGGGGTTGATGATGTGGGCTATTTCGCCTAAGTCCCGGGCGTTATCGCGAAGAATCGCCGGATTTTGCGGAGGGGCTCCTGGTAAg                                                >  6:239680/1‑149 (MQ=255)
                                                                                                 gatgtTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCaaa                                                                                                       >  5:22566/1‑94 (MQ=255)
                                                                                                 gatgtTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCaaa                                                                                                       <  6:22566/94‑1 (MQ=255)
                                                                                                 gatggTCTCTTACGTCATGTTAGAAGTTAACCGTGTGGCAGGTCGGCAGCGCCCGGGCGATGGGGGCTTTGATGATGACGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTAGTGGaaa                                                 <  2:252632/148‑1 (MQ=255)
                                                                                                    gtTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGtt                                             >  3:325927/1‑149 (MQ=255)
                                                                                                     ttCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTTAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTc                                            <  2:349887/149‑1 (MQ=255)
                                                                                                     ttCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTTAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTc                                            <  2:344923/149‑1 (MQ=255)
                                                                                                     ttCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCTTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGtt                                             <  6:218383/148‑1 (MQ=255)
                                                                                                       ctctTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCTGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCt                                           <  8:243685/148‑1 (MQ=255)
                                                                                                        tctTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTgg                                         >  3:308511/1‑149 (MQ=255)
                                                                                                          ttACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTgg                                         <  4:257916/147‑1 (MQ=255)
                                                                                                            aCGGCATGTTAGAAGTTAACCGTCTTGCAGGTCGTCAGCTCCCGGTCGGTGGTGGCTTTGATTATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCa                                       <  3:310588/147‑1 (MQ=255)
                                                                                                            aCGGCATGTTAGAAGTTAACCGGCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGc                                        <  4:309023/146‑1 (MQ=255)
                                                                                                                  tggtAGAAGTTACCCGGCTGGCAGGTCTTCAGCGCCCGGTCGCTGGGGGCTTTGATGATGAGGGCAAGTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGGCGAtt                               <  5:202788/146‑1 (MQ=255)
                                                                                                                  tGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGAtt                               >  6:181268/1‑149 (MQ=255)
                                                                                                                  tGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGAtt                               <  4:308511/149‑1 (MQ=255)
                                                                                                                  tGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGAtt                               <  1:89378/149‑1 (MQ=255)
                                                                                                                  tGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGAtt                               <  1:141169/149‑1 (MQ=255)
                                                                                                                    ttAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGAt                                                              >  1:294663/1‑116 (MQ=255)
                                                                                                                    ttAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGAt                                                              <  2:294663/116‑1 (MQ=255)
                                                                                                                    ttAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATtg                              <  1:201064/148‑1 (MQ=255)
                                                                                                                       gAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGc                           <  1:76624/148‑1 (MQ=255)
                                                                                                                        aaGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTg                         >  3:256884/1‑149 (MQ=255)
                                                                                                                         aGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCAACTGGAAAGGTTCTGGCATGTCGATTGTGCTgg                        <  2:59659/149‑1 (MQ=255)
                                                                                                                                             gtcAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCAt                                      <  6:191150/115‑1 (MQ=255)
                                                                                                                                             gtcAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCAt                                      >  5:191150/1‑115 (MQ=255)
                                                                                                                                              tcAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCAGGTCGATTGTGCTGGATATGATCGCTACTCTCCttt   >  5:306866/1‑149 (MQ=255)
                                                                                                                                               cAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTc  >  8:197332/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
TGCCGCCGTGGGGCGCAAAAGAGTGTCGCATAGGCACTAACCCGCTGATCGTCGCCATTCCTTCCACGCCGATCACCATGGTCGATATGTCGATGTCGATGTTCTCTTACGGCATGTTAGAAGTTAACCGTCTGGCAGGTCGTCAGCTCCCGGTCGATGGTGGCTTTGATGATGAGGGCAATTTGACCAAAGAACCTGGCGTTATCGAGAAGAATCGCCGCATTTTGCCGATGGGCTACTGGAAAGGTTCTGGCATGTCGATTGTGCTGGATATGATCGCTACTCTCCTTTC  >  NC_000913/3743163‑3743454

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑