Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,040,032 T→G L246R (CTG→CGG)  cbdB → cytochrome bd‑II oxidase, subunit II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,040,0320TG96.8% 190.4 / ‑1.8 63L246R (CTG→CGG) cbdBcytochrome bd‑II oxidase, subunit II
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base G (32/28);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 2.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.80e-01

CATCTGCGTTCACAGCTGGCGACCAAACGCGCTGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATC  >  NC_000913/1039884‑1040179
                                                                                                                                                    |                                                                                                                                                   
cATCTGCGTTCACAGCTGGCGACCAAACGCGCTGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACggg                                                                                                                                                   >  1:1180599/1‑151 (MQ=255)
   cTGCGTTCACAGCTGGCGACCAAACGCGCTGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACgggtgg                                                                                                                                                >  1:484687/1‑151 (MQ=255)
     gCGTTCACAGCTTGCGACCAAACCCGCTGCACTCTTCGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGTCTGTGGGTCGGTATTGATGCCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACgggtg                                                                                                                                                 <  1:995215/148‑1 (MQ=255)
     gCGTTCACAGCTGGCGACCAAACGCGCTGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACGGGTGGCa                                                                                                                                              <  1:134307/151‑1 (MQ=255)
     gCGTTCACAGCTGGCGACCAAACGCGCTGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACGGGTGGCa                                                                                                                                              <  2:1033054/151‑1 (MQ=255)
       gTTCACAGCTGGCGACCAAACGCGCTGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACGGGTGGCAGt                                                                                                                                            <  1:20995/151‑1 (MQ=255)
       gTTCACAGCTGGCGACCAAACGCGCTGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACGGGGGGCAg                                                                                                                                             >  1:494560/1‑150 (MQ=255)
       gTTCACAGCTGGCGACCAAACGCGCTGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGGTACTGGCCGTGGGGTGGTATTGGTGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGTAACGGGTGGCa                                                                                                                                              >  2:378341/1‑149 (MQ=255)
        ttCACAGCTGGCGACCAAACGCGCTGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGGTAATGGCTGTGGGGCGGTATTGGTGGCTTTTTGCTGCTCGCCCAGGGTGGTTACGGGCCCTCCAATCCGTTAATGAAAAAGGGGGCAGTg                                                                                                                                           >  2:736725/1‑151 (MQ=255)
         tCACAGCTGGCGACCAAACGCGCTGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTg                                                                                                                                           <  1:953826/150‑1 (MQ=255)
             aGCTGGCGACCAAACGCGCTGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGGTACTGGGTGTGGGGCGGGATTGATGGCTTTGTTCTGCTCGGCCAAGATGCTAAAGGTCCTTTCAATCCGTTAATGAAACGGGTGGCCGTg                                                                                                                                           >  1:61187/1‑146 (MQ=255)
                tGGCGACCAAACGCGCTGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTCGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGGTGCTAACGGTCCTTCCAATCCGTTAATGAAACGGGTGGCAGTGCTACCCgg                                                                                                                                   >  2:310844/1‑151 (MQ=255)
                 ggCGCCCTAACGCGCTGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGTCTGTGTTTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCGATCCGTTAATGAAACGGTTGGCAGTGCTACCTGGt                                                                                                                                  <  2:1160779/151‑1 (MQ=255)
                 ggCGACCAAACGCGCTGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCCGTGGGGCGGTATTGATGGCTTTGTACTGCCCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAAAGGGTGGCCGGGCTACCTgg                                                                                                                                   >  2:581134/1‑150 (MQ=255)
                   cGACCAAACGCGCTGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGGATTTATGGCTTTTTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACGGGTGGCAGTGCTACCTGGTg                                                                                                                                 >  1:843672/1‑150 (MQ=255)
                    gACCAAACGCGCTGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGGCGGTATTGATGGGTTTGTACTGCTCGCCCCGGATGCTTACGGGTCTTCCAAACCGTTAATGAAACGGGGGGCAGTGCTACCTgg                                                                                                                                   >  2:681757/1‑147 (MQ=255)
                        aaaCGCGCTGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACGGGTGGCAGTGCTACCTGGTGCCTgg                                                                                                                            <  1:71581/150‑1 (MQ=255)
                              gcTGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACGGGTGGCAGTGCTACCTGGTGCCTGGATGaat                                                                                                                      <  2:1069914/150‑1 (MQ=255)
                                tGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACGGGGGGCAGGGCTACCTGGGGCCTGGGTGaataat                                                                                                                   >  1:865751/1‑151 (MQ=255)
                                tGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGGATTGATGGCTTTGTACTGCTCGCCCCGGGTGCTAACGGTCCTTCCAATCCGGTAATGAAACGGGGGGCAGTGCTACCTGGGGCCTGGATGgataa                                                                                                                    >  2:662126/1‑150 (MQ=255)
                                 gCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACGGGTGGCAGTGCTACCTGGTGCCTGGATGAATAAtt                                                                                                                  <  1:94049/151‑1 (MQ=255)
                                   aCTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGGATTGATGGCTTTGTTCTGCTCGCCCAGGATGCTAACGGGCCTTCCAATCCGTTAATGAAACGGGTGGCAGGGCTACCTGGGGCCTGGGTGAATAAtttt                                                                                                                >  2:342817/1‑151 (MQ=255)
                                   aCTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGGTACTGGCTGTGGGTCGGGATTGGTGGCTTTGTACTGCTCGCCCAGGGTGCTAACGGTCCTTCCAATCCGTTAATGGAACGGGTGGCAGTGCTACCTTGGGCCTGGGTTAATTAtttt                                                                                                                >  2:370117/1‑151 (MQ=255)
                                    cTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACGGGTGGCAGTGCTACCTGGTGCCTGGATGAATAAtttt                                                                                                                >  1:507539/1‑150 (MQ=255)
                                    cTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGGCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGGTAATGAAACGGGGGGCAGGGCTACCTGGTGCCTGGATGAATAAtttt                                                                                                                >  2:850867/1‑150 (MQ=255)
                                    cTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGATACTGGCTGTGGGTCGGGATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACGGGGGGCAGTGCTACCTTGTGCCTGGATGAATAAtttt                                                                                                                >  2:103404/1‑150 (MQ=255)
                                      gTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCCGGGTGCTAACGGTCCTTCCAATCCGTTAATGGAACGGGTGGCAGTGCTACCTGGGGCCTGGATGAATAAtttt                                                                                                                >  2:107609/1‑148 (MQ=255)
                                      gTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACGGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTc                                                                                                             <  1:368919/151‑1 (MQ=255)
                                      gTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGGATTGATTGCTTTTTACTTCTCGCCCAGGGTGCTAACGGTCCTTCCAATCCGTTTATTAAAAGGGGGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGt                                                                                                              >  2:851030/1‑150 (MQ=255)
                                        tGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACGGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCg                                                                                                            <  1:1161875/150‑1 (MQ=255)
                                         ggTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACGGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGa                                                                                                           <  2:139102/150‑1 (MQ=255)
                                           tGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGGATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACGGGTGGCCGTGGTACCTGGGGCCTTGATGAATAATTTTGt                                                                                                              >  2:978021/1‑145 (MQ=255)
                                             aTGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACGGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCg                                                                                                       <  1:974916/150‑1 (MQ=255)
                                             aTGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACGGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGc                                                                                                      <  2:109915/151‑1 (MQ=255)
                                                  gtgCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACGGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCg                                                                                                   <  1:572378/149‑1 (MQ=255)
                                                       ttttGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACGGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTtg                                                                                              <  2:1045383/149‑1 (MQ=255)
                                                              ggCGGGTTACTGGCTGTTGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAATCGGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCt                                                                                       <  2:660022/149‑1 (MQ=255)
                                                                   gTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACGGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCgct                                                                                <  2:1054955/151‑1 (MQ=255)
                                                                     tACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACGGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCgctg                                                                               <  2:296500/150‑1 (MQ=255)
                                                                     tACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACGGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCgctg                                                                               <  2:873891/150‑1 (MQ=255)
                                                                      aCTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACGGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCgctgct                                                                             >  1:104878/1‑151 (MQ=255)
                                                                         ggCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGctgctg                                                                            <  1:492723/149‑1 (MQ=255)
                                                                                   cGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACGGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTg                                                                 <  1:1142400/150‑1 (MQ=255)
                                                                                   cGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACGGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGc                                                                <  1:209574/151‑1 (MQ=255)
                                                                                    ggTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACGGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGcc                                                               >  2:149754/1‑151 (MQ=255)
                                                                                        ttGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACGGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCAtt                                                           >  1:270814/1‑151 (MQ=255)
                                                                                          gATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACGGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTg                                                          <  1:636439/150‑1 (MQ=255)
                                                                                                  tGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACGGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGt                                                  <  1:239688/150‑1 (MQ=255)
                                                                                                  tGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACGGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTg                                                 >  2:469294/1‑151 (MQ=255)
                                                                                                   gTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACGGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTGa                                                <  2:527990/151‑1 (MQ=255)
                                                                                                    tACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACGGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTGAt                                               >  2:907264/1‑151 (MQ=255)
                                                                                                     aCTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACGGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTGAt                                               <  1:423138/150‑1 (MQ=255)
                                                                                                         cTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACGGGGGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGATTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGGGATGGCg                                           >  1:884843/1‑150 (MQ=255)
                                                                                                           cGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACGGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTGATGGCGAtt                                        >  1:3943/1‑151 (MQ=255)
                                                                                                           cGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACGGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTGATGGCGAt                                         <  2:1166340/150‑1 (MQ=255)
                                                                                                                aGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACGGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTGATGGCGATTTATCg                                   >  2:753333/1‑151 (MQ=255)
                                                                                                                 ggATGCTAACGGTCCTTCCAATCCGTTAATGAAACGGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTGATGGCGATTTATCGt                                  >  1:656886/1‑151 (MQ=255)
                                                                                                                                  ccAATCCGTTAATGAAACGGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTggg                  <  1:1143905/150‑1 (MQ=255)
                                                                                                                                     aTCCGTTAATGAAACGGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGAtt              >  1:929503/1‑151 (MQ=255)
                                                                                                                                        cGTTAATGAAACGGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGAttttt           <  1:75915/151‑1 (MQ=255)
                                                                                                                                          ttAATGAAACGGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGAttttt           >  1:400336/1‑149 (MQ=255)
                                                                                                                                           tAATGAAACGGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGa         <  1:378358/150‑1 (MQ=255)
                                                                                                                                           tAATGAAACGGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGAt        >  2:1170151/1‑151 (MQ=255)
                                                                                                                                            aaTGAAACGGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCACGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATg       >  1:1140512/1‑151 (MQ=255)
                                                                                                                                               gAAACGGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGc     >  1:322930/1‑150 (MQ=255)
                                                                                                                                                 aaCGGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGcatc  >  1:633612/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CATCTGCGTTCACAGCTGGCGACCAAACGCGCTGCACTGTTGGTGATGCTGTGCTTTTTGCTGGCGGGTTACTGGCTGTGGGTCGGTATTGATGGCTTTGTACTGCTCGCCCAGGATGCTAACGGTCCTTCCAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTGATGGCGATTTATCGTGGTCGCCCGGGTTGGGGATTTTTGATGGCATC  >  NC_000913/1039884‑1040179

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: