Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,334,705 C→T 6.8% intergenic (+24/‑204) ispB → / → sfsB octaprenyl diphosphate synthase/malPQ operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,334,7050CT 6.8% 389.6 / 10.1 131intergenic (+24/‑204)ispB/sfsBoctaprenyl diphosphate synthase/malPQ operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (59/62);  new base T (6/3);  total (66/65)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

Only 100 of 136 total aligned reads displayed.

GGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGATG  >  NC_000913/3334557‑3334854
                                                                                                                                                    |                                                                                                                                                     
ggACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAgcg                                                                                                                                                      <  2:727890/150‑1 (MQ=255)
ggACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAgcg                                                                                                                                                      >  2:719296/1‑150 (MQ=255)
ggACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAgcg                                                                                                                                                      >  1:1262308/1‑150 (MQ=255)
      gTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCggg                                                                                                                                                <  2:1576346/150‑1 (MQ=255)
       tCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCCCATCCCGCGCGCTGCGCGGGa                                                                                                                                               <  2:996349/150‑1 (MQ=255)
           cGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGGTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCGCTGCGCGGGATgag                                                                                                                                           <  2:669899/150‑3 (MQ=255)
            gTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCGCTGCGCGGGATgagg                                                                                                                                          <  1:182539/150‑4 (MQ=255)
               ccGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCCTCGGCCTCGCGCACATCGCTGTTCACCGCGATCGTTAATCACCTCCCCTCATCACGCGCAGCGCGCGGGATGGTGCCa                                                                                                                                       >  2:42533/1‑150 (MQ=255)
                 gAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCACCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGt                                                                                                                                     >  2:678029/1‑150 (MQ=255)
                   ggAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTaa                                                                                                                                   >  1:764145/1‑150 (MQ=255)
                   ggAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTaa                                                                                                                                   >  2:140031/1‑150 (MQ=255)
                      agaagCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACACCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGtt                                                                                                                                >  2:331455/1‑150 (MQ=255)
                         agCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCa                                                                                                                             >  2:1389462/1‑150 (MQ=255)
                             gACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAGGTTCCATaaa                                                                                                                         <  2:1670667/150‑1 (MQ=255)
                             gACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATaaa                                                                                                                         <  2:1229795/150‑1 (MQ=255)
                               cAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATaaaac                                                                                                                       >  1:234085/1‑149 (MQ=255)
                                     cATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAACCACttt                                                                                                                 >  2:890725/1‑149 (MQ=255)
                                      aTCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTAt                                                                                                                <  1:990915/150‑1 (MQ=255)
                                        cGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTc                                                                                                              >  1:1229188/1‑150 (MQ=255)
                                         gCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCCTCGGCCTCGCGCACATCGCTGTTCACCGCGATCGTTAATCCCCTCCCCTCATCCCGCGCCGCGCGCGCGATGGTTCCAGTAAGTTCCATATCAACACGTATTcg                                                                                                             >  1:305951/1‑149 (MQ=255)
                                          cAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGGTCCAGTAAGTTCCATAAAACCACCTATTGAg                                                                                                            >  2:1066965/1‑150 (MQ=255)
                                           aGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGCTCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGGTCCATAAAAACACTGATTCCGc                                                                                                           >  1:68978/1‑150 (MQ=255)
                                           aGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCCCAGCGCGCCGGGTGGTTACAGTAAGTTCCATAAAAACACTTATTCAGc                                                                                                           >  2:1269412/1‑150 (MQ=255)
                                            gCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTGCAGTAAGTTCCATAAAAACACATATTCAGct                                                                                                          >  2:536585/1‑150 (MQ=255)
                                            gCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGct                                                                                                          <  1:1445577/150‑1 (MQ=255)
                                                    ggTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGCTGGTTCCAGTAAGTTCCATAAATCCCCTGATGCAGCGCtatcccc                                                                                                  >  2:318291/1‑146 (MQ=255)
                                                     gTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCtata                                                                                                     >  1:998619/1‑146 (MQ=255)
                                                      tGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACAcc                                                                                                >  2:1464009/1‑150 (MQ=255)
                                                        cTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTg                                                                                              <  2:1090744/150‑1 (MQ=255)
                                                         tCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGa                                                                                             <  1:584511/150‑1 (MQ=255)
                                                          cccGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGaa                                                                                            >  1:1194342/1‑150 (MQ=255)
                                                          cccGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGaa                                                                                            >  2:278797/1‑150 (MQ=255)
                                                               acacCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATa                                                                                       <  2:907387/150‑1 (MQ=255)
                                                                cacCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAg                                                                                      <  1:1710186/150‑1 (MQ=255)
                                                                 acCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAgg                                                                                     <  1:1017857/150‑1 (MQ=255)
                                                                 acCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAgg                                                                                     <  1:1464009/150‑1 (MQ=255)
                                                                  ccccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCCCCGCGGGCGGGGTGGTTGCCGTAAGTTGCAAAAAAACCAATATTAAGCTCTATAACACCGGTAAAATAtgg                                                                                    >  2:1323047/1‑150 (MQ=255)
                                                                   cccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATGCCGCCCAGAGCGCGGGAGGGGTCCAGGAAGTTCCCTAAAAACACTTATTCAGCTCAATAACAACTGAAAAATAGGGa                                                                                   >  2:850786/1‑150 (MQ=255)
                                                                   cccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGa                                                                                   <  2:1262308/150‑1 (MQ=255)
                                                                    ccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTGCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAt                                                                                  >  1:1402286/1‑150 (MQ=255)
                                                                     cTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATGCCCCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACCCGTATTAAGCTCGTTAACACCAGAAAAATAGGGAtg                                                                                 >  1:1268160/1‑150 (MQ=255)
                                                                       tGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgtg                                                                               <  2:1123048/150‑1 (MQ=255)
                                                                        ggCGAGAAGCACTCATCGGCCGCGCGCACATCGCTGTTCAACGCGTTCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGTATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGCCATGTGc                                                                              <  2:201380/150‑1 (MQ=255)
                                                                           gagaAGCACTCATCGGCCTCGCGCACATCGCGGGTCGACGCGATCGTTAAGCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATc                                                                           <  2:1093377/150‑1 (MQ=255)
                                                                             gtcgCACTCAGCGGCCTCGCGCACATCGCTGATCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCCGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGGGCATCtt                                                                         <  2:257179/147‑1 (MQ=255)
                                                                              aaGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCttt                                                                        <  1:454563/150‑1 (MQ=255)
                                                                               aGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTta                                                                       <  1:1018423/150‑1 (MQ=255)
                                                                               aGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTta                                                                       >  2:983170/1‑150 (MQ=255)
                                                                               aGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTta                                                                       <  1:54021/150‑1 (MQ=255)
                                                                               aGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTta                                                                       <  1:801747/150‑1 (MQ=255)
                                                                                  aCTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTtatat                                                                    <  2:1118112/150‑1 (MQ=255)
                                                                                  aCTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGGTCCATAAAAACACTTATTCAGCTCTATAACCCCGGAAAAATAGGGATTTGCATATTtatat                                                                    >  1:451835/1‑150 (MQ=255)
                                                                                     cATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCCGCGCGCGGGATGGTTCCAGGAAGTTCCATAAAAACACTTATTCCGCTCTATTACCCCGGAAAAATAGGGATGTGGATCTTTTTATAtt                                                                 >  1:1099165/1‑150 (MQ=255)
                                                                                       tCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTTAGTTCCAGAAAAACACTTATTCAGCTCTATAACACCGGAAAAATAGGGATGTGCCTCTTTATATATTCt                                                               >  1:136516/1‑150 (MQ=255)
                                                                                        cGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATAAATTCTg                                                              <  2:699551/150‑1 (MQ=255)
                                                                                            cTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGGATCTTTATATAATCTGAata                                                          >  2:141492/1‑150 (MQ=255)
                                                                                              cgcgCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCCGCGCGCGGGATGGTTACAGTAAGTTACATAAAAAAACATATTAAGCTCTATAACACCAGAAAAATATGGATGTGCATCTTTATATATTATGAATAtt                                                        >  2:1410704/1‑150 (MQ=255)
                                                                                                cgcACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCCGCGCGCGGGGTGGTTCCAGTAAGTTGCATAAAAACACTTATTGAGGTCTTTAACACCTGAAAAATAGGGATGTGGATCTTTATTTATTATGAAttttga                                                      >  1:573585/1‑148 (MQ=255)
                                                                                                cgcACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCATCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTca                                                      <  1:994205/150‑1 (MQ=255)
                                                                                                cgcACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTca                                                      <  2:1662240/150‑1 (MQ=255)
                                                                                                cgcACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTca                                                      >  1:843046/1‑150 (MQ=255)
                                                                                                cgcACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCGATAACACCTGAAAAATAGGGATGTGCATCTTTATTTATTATGAATAttaa                                                      >  1:678235/1‑148 (MQ=255)
                                                                                                  cacaTCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGa                                                             >  1:902184/1‑141 (MQ=255)
                                                                                                  cacaTCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGa                                                             <  2:902184/141‑1 (MQ=255)
                                                                                                  cacaTCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTcaca                                                    <  1:590482/150‑1 (MQ=255)
                                                                                                   acaTCGCTGTTCAACGCGATCGTTAATCCACTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCacac                                                   >  1:1215680/1‑150 (MQ=255)
                                                                                                      tCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACActct                                                >  1:241271/1‑150 (MQ=255)
                                                                                                         cTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTa                                             <  1:719296/150‑1 (MQ=255)
                                                                                                           gTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCAT‑‑AAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAgg                                         >  2:1362607/1‑150 (MQ=255)
                                                                                                            ttCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCCCAGCGCGCGGGATGGTGCCAGTAAGTTCCATATAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTAAAATATTCTGAAAATTCACACTCTTTACAg                                          >  2:925598/1‑150 (MQ=255)
                                                                                                             tCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCAAGTAAGTTCCATTAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAGGTGCATCTTTATATATTCTGAATATTCACACTCTTGACAgg                                         >  1:1167504/1‑150 (MQ=255)
                                                                                                                actcgATCGTTAATCCCCTCCCATCATACATCGCATCGAGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTAATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAAc                                      <  1:102756/147‑1 (MQ=255)
                                                                                                                 cgcgATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACt                                     <  2:843046/150‑1 (MQ=255)
                                                                                                                  gcgATCGTTAATCCCCTCCCCTCATCCCGCGCCGCGCGCGGGATGGGTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACtt                                    >  1:548514/1‑150 (MQ=255)
                                                                                                                   cgATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTAAGCTCTATAACACCAGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACttt                                   >  1:1593526/1‑150 (MQ=255)
                                                                                                                    gATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACtttt                                  <  1:4560/150‑1 (MQ=255)
                                                                                                                    gATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACtttt                                  <  1:1515129/150‑1 (MQ=255)
                                                                                                                     aTCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACttttt                                 <  1:1389462/150‑1 (MQ=255)
                                                                                                                       cGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTag                               <  2:891991/150‑1 (MQ=255)
                                                                                                                       cGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTag                               <  1:1697531/150‑1 (MQ=255)
                                                                                                                        gTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTaga                              <  1:278797/150‑1 (MQ=255)
                                                                                                                         ttAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTagag                             <  1:1567004/150‑1 (MQ=255)
                                                                                                                               ccctcccctcATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTACAGGAACTTTTTAGAGCAATAgg                      <  2:604792/150‑1 (MQ=255)
                                                                                                                                cctcccctcATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACACTTTACAGGAACTTTTTAGAGCAATAgg                      <  2:817707/150‑1 (MQ=255)
                                                                                                                                cctcccctcATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATATTTATATATTCTGAATATTAACACTCTTTACAGTAACTTTTTAGAGCAATAgg                      >  2:677397/1‑150 (MQ=255)
                                                                                                                                      accATCCCGCGCGCTGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCa                >  1:599088/3‑150 (MQ=255)
                                                                                                                                       tcATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAg               >  2:954821/1‑150 (MQ=255)
                                                                                                                                        cATCCCGCGCGCTGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAgg              >  2:657591/1‑150 (MQ=255)
                                                                                                                                        cATCCCGCGCGCTGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCATTACAGGAACTTTTTAGAGCAATAGGCCATCAgg              >  2:651616/1‑150 (MQ=255)
                                                                                                                                         aTCCCGCGCGCTGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTAATCAGCTCTATAACACCAGAAAAATAGGGAGGTGCAACTTTATATATTAGCAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAggt             >  1:288035/1‑149 (MQ=255)
                                                                                                                                         aTCCCGCGCGCTGCGCGGGATGGTTCCAGTAAGTTCCATAAAAAAACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATAAAATCTGAATATTGACACTCTTTACAGGAACTTTTTAGAGCTATAGGCCATCAGGa             >  1:49500/1‑150 (MQ=255)
                                                                                                                                         aTCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGa             >  2:1697641/1‑150 (MQ=255)
                                                                                                                                          tCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGATCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAg            <  1:920135/150‑1 (MQ=255)
                                                                                                                                           cccGCGCAGCGCGCGGGATGGTTCCAGTAATTTCCATAAATACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGt           <  1:1128470/150‑1 (MQ=255)
                                                                                                                                           cccGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACTCCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGATCTTTTTAGAGCAATAGGCCATCAGGAGt           <  1:612905/150‑1 (MQ=255)
                                                                                                                                           cccGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGt           <  1:1092494/150‑1 (MQ=255)
                                                                                                                                           cccGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGt           <  2:136785/150‑1 (MQ=255)
                                                                                                                                           cccGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCAATCAGGAGt           >  2:225548/1‑150 (MQ=255)
                                                                                                                                            ccGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGta          >  1:633624/1‑150 (MQ=255)
                                                                                                                                                  aGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGa    <  2:550016/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTGCTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGCGCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGATG  >  NC_000913/3334557‑3334854

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: