Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,227,968 T→G 6.6% V56V (GTT→GTG ygjJ → putative periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,227,9680TG 6.6% 585.3 / 8.7 185V56V (GTT→GTGygjJputative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base T (98/72);  new base G (2/13);  total (100/85)
Fisher's exact test for biased strand distribution p-value = 9.92e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 190 total aligned reads displayed.

GCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTCGGCCTGATTTACGGCGAAGGCTGGACCGTGCAGGCCGACCACGACGGCAATG  >  NC_000913/3227823‑3228116
                                                                                                                                                 |                                                                                                                                                    
gCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTtggtg                                                                                                                                                  <  1:1798876/150‑1 (MQ=255)
 cAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTtggtgg                                                                                                                                                 <  1:1951077/150‑1 (MQ=255)
  agagCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTGGGTGGa                                                                                                                                                >  2:371019/1‑150 (MQ=255)
    agCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGata                                                                                                                                              <  2:2357424/150‑1 (MQ=255)
    agCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGata                                                                                                                                              <  2:3169645/150‑1 (MQ=255)
    agCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGata                                                                                                                                              <  1:169822/150‑1 (MQ=255)
     gCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGatat                                                                                                                                             <  2:712763/150‑1 (MQ=255)
     gCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGatat                                                                                                                                             <  2:2026941/150‑1 (MQ=255)
       tttACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGatatat                                                                                                                                           <  1:943838/149‑1 (MQ=255)
       aTTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGatatat                                                                                                                                           <  1:3830258/150‑1 (MQ=255)
         tACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATgg                                                                                                                                         >  2:2964293/1‑150 (MQ=255)
           cTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGAt                                                                                                                                       >  2:2461506/1‑150 (MQ=255)
           cTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGAt                                                                                                                                       >  2:2361030/1‑150 (MQ=255)
            ttGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCACTTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGCATTAGGGGGGGGTGGATATAGGGagg                                                                                                                                      >  2:3785382/1‑148 (MQ=255)
            ttGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGAtt                                                                                                                                      >  1:2646809/1‑150 (MQ=255)
               ctctGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTag                                                                                                                                   >  2:2828628/1‑150 (MQ=255)
                 ctGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTagag                                                                                                                                 <  1:818481/150‑1 (MQ=255)
                 ctGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTagag                                                                                                                                 <  1:252525/150‑1 (MQ=255)
                 ctGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTGGGTGGATATATGGATTTagag                                                                                                                                 <  1:1773716/150‑1 (MQ=255)
                  tGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAgg                                                                                                                                <  1:2268365/150‑1 (MQ=255)
                    ccGTTGTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGt                                                                                                                              <  2:3592821/150‑1 (MQ=255)
                       ttttGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTgag                                                                                                                           <  2:2302889/150‑1 (MQ=255)
                       ttttGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGACCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCAACTCCTTCAATTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTgag                                                                                                                           >  1:2366796/1‑150 (MQ=255)
                         ttGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGaa                                                                                                                         <  2:2267782/150‑1 (MQ=255)
                         ttGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGaa                                                                                                                         <  1:3506430/150‑1 (MQ=255)
                          tGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAAt                                                                                                                        <  2:1895626/150‑1 (MQ=255)
                           gCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATa                                                                                                                       >  2:3822162/1‑150 (MQ=255)
                            cTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACTTCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGTGGGTGAGAATaa                                                                                                                      >  1:304832/1‑150 (MQ=255)
                             tACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATaaa                                                                                                                     <  1:2230462/150‑1 (MQ=255)
                             tACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATaaa                                                                                                                     <  2:3744529/150‑1 (MQ=255)
                              aCGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAAc                                                                                                                    >  2:347354/1‑150 (MQ=255)
                               cGCCTTTTCTGCGGCAGGAGAAGAAGCTCGTCCGGCAGAACATGACGACACAAACACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACa                                                                                                                   <  2:3583812/150‑1 (MQ=255)
                                    tttCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAAt                                                                                                              >  2:4096269/1‑150 (MQ=255)
                                      tCTGCGGCAGGAGAAGAAGCACGTCCGGCATAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATAc                                                                                                            <  2:1095386/150‑1 (MQ=255)
                                        tGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAg                                                                                                          <  1:1815557/150‑1 (MQ=255)
                                           ggCTGGAGAAGAAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCTTTTTTACGGCTAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGa                                                                                                       <  1:3449003/150‑1 (MQ=255)
                                           ggCAGGAGAAGAAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGa                                                                                                       >  1:3214133/1‑150 (MQ=255)
                                              aGGAGAAGAAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACgg                                                                                                    >  1:3283396/1‑150 (MQ=255)
                                                gagaAGAAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTa                                                                                                  >  2:257044/1‑150 (MQ=255)
                                                 agaagaAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTAc                                                                                                 >  2:1342413/1‑150 (MQ=255)
                                                 agaagaAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTAc                                                                                                 >  2:3395860/1‑150 (MQ=255)
                                                  taagaagCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCTAATTAGGGTTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTAcc                                                                                                <  2:4032811/149‑1 (MQ=255)
                                                  gaagaaGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTAcc                                                                                                <  2:4032839/150‑1 (MQ=255)
                                                    agaagCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCta                                                                                              <  2:2888235/150‑1 (MQ=255)
                                                     gaagCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCtat                                                                                             >  2:3451044/1‑150 (MQ=255)
                                                        gCACGTCCGGCAGAACATGACGACACAAAAACACCCGAAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATAtt                                                                                          <  2:142807/150‑1 (MQ=255)
                                                              ccGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGt                                                                                    <  2:2008030/150‑1 (MQ=255)
                                                                  cAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGcc                                                                                >  1:3912553/1‑150 (MQ=255)
                                                                   aGAACATGACGACACAAAAACACCCGCATTTACCTCGACATCTCCTCCTTCATTTCTTTTTTACGGCGATTTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCt                                                                               <  1:2474574/150‑1 (MQ=255)
                                                                      aCATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAg                                                                                      <  2:326353/140‑1 (MQ=255)
                                                                      aCATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAg                                                                                      >  1:326155/1‑140 (MQ=255)
                                                                      aCATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGAGTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGga                                                                            >  2:1562311/1‑150 (MQ=255)
                                                                               cacaAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGGGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGTGAATATACATAAATACAGCGACGGTACCTATATTGTAGGTGGCGTGGAGATGAAGTa                                                                   >  1:3580739/1‑150 (MQ=255)
                                                                               cacaAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTa                                                                   >  1:681112/1‑150 (MQ=255)
                                                                               cacaAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTa                                                                   >  1:593545/1‑150 (MQ=255)
                                                                                acaAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTAc                                                                  >  1:893991/1‑150 (MQ=255)
                                                                                         cccGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTgg                                                         >  1:1380723/1‑150 (MQ=255)
                                                                                          ccGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGt                                                        >  1:646345/1‑150 (MQ=255)
                                                                                          ccGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACATATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGt                                                        >  1:279318/1‑150 (MQ=255)
                                                                                            gCAATTACCTCTACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTc                                                      <  2:3399450/150‑1 (MQ=255)
                                                                                            gCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAACCATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTc                                                      <  2:2742805/150‑1 (MQ=255)
                                                                                            gCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGGGGGTGGATATATGGATTTAGCGGGTGAGAATAAACATAACTACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTc                                                      <  1:1165090/150‑1 (MQ=255)
                                                                                               aTTCCCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTCGGc                                                   <  2:515444/150‑1 (MQ=255)
                                                                                               aTTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTTGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTCGGc                                                   <  1:1491760/150‑1 (MQ=255)
                                                                                               aTTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTCGGc                                                   <  2:1731993/150‑1 (MQ=255)
                                                                                               aTTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTCGGc                                                   <  2:3308574/150‑1 (MQ=255)
                                                                                               aTTACCACGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTTGGGGGTGGTGGTTATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTCGGc                                                   <  1:2895960/150‑1 (MQ=255)
                                                                                                tttccTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAATCATAAATACAGCTACGGTACGTATATTGATGGTGGCGTGGAGATGTAGTACGGCTCCTTGTTCGGcg                                                  >  2:851451/4‑149 (MQ=255)
                                                                                                   ccTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTCGGCCTGa                                               >  1:1154299/1‑150 (MQ=255)
                                                                                                    cTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTCGGCCTGAt                                              >  1:1037158/1‑150 (MQ=255)
                                                                                                    cTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTCGGCCTGAt                                              >  1:2021340/1‑150 (MQ=255)
                                                                                                    cTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTCGGCCTGAt                                              >  1:3708627/1‑150 (MQ=255)
                                                                                                    cTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCGGGAGATGAAGTACGGCTCCTGGTTCGGCCTGAt                                              >  1:1037096/1‑150 (MQ=255)
                                                                                                     tCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTCGGCCTGAtt                                             >  1:841357/1‑150 (MQ=255)
                                                                                                     tCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTCGGCCTGAtt                                             >  1:688044/1‑150 (MQ=255)
                                                                                                     tCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTCGGCCTGAtt                                             >  1:3842878/1‑150 (MQ=255)
                                                                                                             ttctCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTCGGCCTGATTTACGGCGa                                     <  2:1585573/150‑1 (MQ=255)
                                                                                                                tcCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTCGGCCTGATTTACGGCGAAgg                                  >  2:3064244/1‑150 (MQ=255)
                                                                                                                       tttCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTGCCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTCGGCCTGATTTACGGCGAAGGCTGGAcc                           >  2:2361430/1‑150 (MQ=255)
                                                                                                                       tttCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTCGGCCTGATTTACGGCGAAGGCTGGAcc                           >  1:1515709/1‑150 (MQ=255)
                                                                                                                        ttCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTCGGCCTGATTTACGGCGAAGGCTGGACCg                          >  1:3182465/1‑150 (MQ=255)
                                                                                                                        ttCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTCGGCCTGATTTACGGCGAAGGCTGGACCg                          >  1:3182319/1‑150 (MQ=255)
                                                                                                                        ttCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTCGGCCTGATTTACGGCGAAGGCTGGACCg                          >  1:3182275/1‑150 (MQ=255)
                                                                                                                           gTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTCGGCCTGATTTACGGCGAAGGCTGGACCGTGc                       >  2:675561/1‑150 (MQ=255)
                                                                                                                               ttACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTCGGCCTGATTTACGGCGAAGGCTGGACCGTGCAGGc                   >  2:475905/1‑150 (MQ=255)
                                                                                                                                tACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATCCAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTCGGCCTGATTTACGGCGAAGGCTGGACCGTGCAGGcc                  <  2:1221808/150‑1 (MQ=255)
                                                                                                                                 aCGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAGCATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTCGGCCTGATTTACGGGGAAGGCTGGACCGTGCAGGCCg                 <  2:3131697/150‑1 (MQ=255)
                                                                                                                                 aCGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTCGGCCTGATTTACGGCGAAGGCTGGACCGTGCAGGCCg                 >  2:3587934/1‑150 (MQ=255)
                                                                                                                                  cGGCGAATTAGGGGGGGGTGGATCTAGGGATTTAGAGGGTGAGAATAACCATACATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTCGGCCTGATTTACGGCGAAGGCTGGACCGTGCAGGCCGa                <  1:4025861/150‑1 (MQ=255)
                                                                                                                                   ggCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTCGGCCTGATTTACGGCGAAGGCTGGACCGTGCAGGCCGAc               >  2:151538/1‑150 (MQ=255)
                                                                                                                                   ggCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTCGGCCTGATTTACGGCGAAGGCTGGACCGTGCAGGCCGAc               >  2:1939962/1‑150 (MQ=255)
                                                                                                                                    gCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTCGGCCTGATTTACGGCGAAGGCTGGACCGTGCAGGCCGAcc              <  1:2714022/150‑1 (MQ=255)
                                                                                                                                    gCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTCGGCCTGATTTACGGCGAAGGCTGGACCGTGCAGGCCGAcc              >  1:1861001/1‑150 (MQ=255)
                                                                                                                                    gCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTCGGCCTGATTTACGGCGAAGGCTGGACCGTGCAGGCCGAcc              >  1:1799379/1‑150 (MQ=255)
                                                                                                                                    gCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTCGGCCTGATTTACGGCGAAGGCTGGACCGTGCAGGCCGAcc              <  1:1342413/150‑1 (MQ=255)
                                                                                                                                    gCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTCGGCCTGATTTACGGCGAAGGCTGGACCGTGCAGGCCGAcc              >  2:206208/1‑150 (MQ=255)
                                                                                                                                    gCGAATTAGGGGTGGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTCGGCCTGATTTACGGCGAAGGCTGGACCGTGCAGGCCGAcc              <  1:3395860/150‑1 (MQ=255)
                                                                                                                                        aTTAGGGGTGGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTCGGCCTGATTTACGGCGAAGGCTGGACCGTGCAGGCCGACCacga          <  1:1379699/150‑1 (MQ=255)
                                                                                                                                          tAGGGGTGGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTGTATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTCGGCCTGATTTACGGCGAAGGCTGGACCGTGCAGGCCGACCacgacg        <  1:3077851/150‑1 (MQ=255)
                                                                                                                                          tAGGGGGGGGGGGATATATGGATTGAGAGGGTGAGAATAAACATAACTACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTCGGCCTGATTTACGGCGAAGGCTAGACCGTGCAGGCCGACCacgacg        <  1:2827999/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GCAGAGCATTACTTGCTCTGCCGTTTTGCTACGCCTTTTCTGCGGCAGGAGAAGAAGCACGTCCGGCAGAACATGACGACACAAAAACACCCGCAATTACCTCGACATCTTCTCCTTCATTTCGTTTTTACGGCGAATTAGGGGTTGGTGGATATATGGATTTAGAGGGTGAGAATAAACATAAATACAGCGACGGTACCTATATTGAAGGTGGCCTGGAGATGAAGTACGGCTCCTGGTTCGGCCTGATTTACGGCGAAGGCTGGACCGTGCAGGCCGACCACGACGGCAATG  >  NC_000913/3227823‑3228116

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: