Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,226,115 A→C 5.8% intergenic (‑35/‑485) minC ← / → ycgJ inhibitor of FtsZ ring polymerization/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,226,1150AC 5.8% 1547.4 / 43.4 510intergenic (‑35/‑485)minC/ycgJinhibitor of FtsZ ring polymerization/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (215/258);  new base C (22/14);  total (237/273)
Fisher's exact test for biased strand distribution p-value = 8.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 517 total aligned reads displayed.

CCTGAGCGATTTTGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGTCAATCAAC  >  NC_000913/1225966‑1226262
                                                                                                                                                     |                                                                                                                                                   
ccTGAGCGATTTTGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCa                                                                                                                                                     <  1:2236710/150‑1 (MQ=255)
ccTGAGCGATTTTGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCa                                                                                                                                                     <  2:8169312/150‑1 (MQ=255)
        aTTTTGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCga                                                                                                                                             <  1:2440136/150‑1 (MQ=255)
        aTTTTGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCga                                                                                                                                             <  1:5576265/150‑1 (MQ=255)
             gTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGatgatg                                                                                                                                        >  2:6950602/1‑150 (MQ=255)
              tCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCGCGCTGGCgatgatta                                                                                                                                       <  1:750681/150‑3 (MQ=255)
              tCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGt                                                                                                                                       <  1:3381709/150‑1 (MQ=255)
              tCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGt                                                                                                                                       >  2:3990497/1‑150 (MQ=255)
                ttCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCCGCGCGCGATGATGTTc                                                                                                                                     >  1:3245244/1‑150 (MQ=255)
                ttCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGCTGTTc                                                                                                                                     >  2:7879873/1‑150 (MQ=255)
                ttCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTc                                                                                                                                     >  2:1507955/1‑150 (MQ=255)
                 tCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTcc                                                                                                                                    >  1:1014344/1‑150 (MQ=255)
                   cAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGa                                                                                                                                  >  2:5170974/1‑150 (MQ=255)
                    aGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGaa                                                                                                                                 <  2:636054/150‑1 (MQ=255)
                    aGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGaa                                                                                                                                 >  1:4502310/1‑150 (MQ=255)
                    aGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGaa                                                                                                                                 <  1:4292559/150‑1 (MQ=255)
                            aTGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATaa                                                                                                                         <  2:9102679/150‑1 (MQ=255)
                            aTGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTCGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATaa                                                                                                                         >  2:8680262/1‑150 (MQ=255)
                             tGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAg                                                                                                                        <  1:4906170/150‑1 (MQ=255)
                               gATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGc                                                                                                                      <  1:7608641/150‑1 (MQ=255)
                                  aaCCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATg                                                                                                                   >  2:4837534/1‑150 (MQ=255)
                                  aaCCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATg                                                                                                                   <  1:294024/150‑1 (MQ=255)
                                      ttAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTtat                                                                                                               >  1:9583643/1‑150 (MQ=255)
                                            tCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTg                                                                                                         <  1:5121161/150‑1 (MQ=255)
                                              tGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGa                                                                                                       >  2:9399049/1‑150 (MQ=255)
                                                ccTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGAtt                                                                                                     >  2:9633781/1‑150 (MQ=255)
                                                  tCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATtat                                                                                                   >  2:6807550/1‑150 (MQ=255)
                                                    aTGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTatat                                                                                                 >  2:7468866/1‑150 (MQ=255)
                                                    aTGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTatat                                                                                                 >  1:6430126/1‑150 (MQ=255)
                                                     tGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATAtt                                                                                                <  1:805859/150‑1 (MQ=255)
                                                      gCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTg                                                                                               >  1:8272864/1‑150 (MQ=255)
                                                       cAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGa                                                                                              <  2:8831383/150‑1 (MQ=255)
                                                       cAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGa                                                                                              <  2:3834259/150‑1 (MQ=255)
                                                        aGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAg                                                                                             <  2:2794337/150‑1 (MQ=255)
                                                        aGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAg                                                                                             <  1:8557254/150‑1 (MQ=255)
                                                               ccACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGGTATAGTCTGGATTATATTGAGGCAAGTc                                                                                      >  2:2613534/1‑150 (MQ=255)
                                                                cacaGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCa                                                                                     >  1:8019990/1‑150 (MQ=255)
                                                                   aGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCAccc                                                                                  <  1:9408119/150‑1 (MQ=255)
                                                                     aTAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTc                                                                                <  1:5922003/150‑1 (MQ=255)
                                                                     aTAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTc                                                                                <  2:3269820/150‑1 (MQ=255)
                                                                     aTAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATGAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTc                                                                                <  2:2123647/150‑1 (MQ=255)
                                                                          gTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCAtt                                                                           >  2:4267536/1‑150 (MQ=255)
                                                                           tGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCAttt                                                                          <  1:2806756/150‑1 (MQ=255)
                                                                            gAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTa                                                                         >  2:2981046/1‑150 (MQ=255)
                                                                               gCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTc                                                                      >  1:3619189/1‑150 (MQ=255)
                                                                               gCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTc                                                                      <  1:1780555/150‑1 (MQ=255)
                                                                                cTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCa                                                                     >  2:7966017/1‑150 (MQ=255)
                                                                                  aCTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCaga                                                                   <  2:3153232/150‑1 (MQ=255)
                                                                                  aCTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCaga                                                                   <  1:1475956/150‑1 (MQ=255)
                                                                                   cTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCagag                                                                  <  1:7386782/150‑1 (MQ=255)
                                                                                        tttAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTaaaa                                                             <  1:3099028/150‑1 (MQ=255)
                                                                                        tttAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTaaaa                                                             <  2:7142802/150‑1 (MQ=255)
                                                                                                cGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTACTCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTc                                                     <  2:3407436/150‑1 (MQ=255)
                                                                                                 gATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCt                                                    >  1:4812214/1‑150 (MQ=255)
                                                                                                 gATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCt                                                    <  2:1396626/150‑1 (MQ=255)
                                                                                                 gATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCt                                                    >  2:1036545/1‑150 (MQ=255)
                                                                                                    tGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCtgtg                                                 <  1:8906795/150‑1 (MQ=255)
                                                                                                         tgtTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATaaa                                            <  2:7762348/150‑1 (MQ=255)
                                                                                                           tttGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATg                                          >  1:1115266/1‑150 (MQ=255)
                                                                                                           tttGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATg                                          >  1:3930624/1‑150 (MQ=255)
                                                                                                            ttGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATgg                                         <  2:544334/150‑1 (MQ=255)
                                                                                                             tGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTCTAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGgc                                        >  1:6599229/1‑150 (MQ=255)
                                                                                                             tGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGgc                                        <  1:905614/150‑1 (MQ=255)
                                                                                                             tGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCAAGTTATAGTCTGGATTATATTGAGGCAAGTCACACTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATgtc                                        >  1:7668811/1‑148 (MQ=255)
                                                                                                              gACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGgcg                                       >  2:3571889/1‑150 (MQ=255)
                                                                                                               aCATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGgcgc                                      >  1:2128808/1‑150 (MQ=255)
                                                                                                                 aTCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTg                                    <  2:4136143/150‑1 (MQ=255)
                                                                                                                  tCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGa                                   <  1:6033136/150‑1 (MQ=255)
                                                                                                                    cTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGAtt                                 >  1:3867049/1‑150 (MQ=255)
                                                                                                                      ggCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCa                               <  1:3457626/150‑1 (MQ=255)
                                                                                                                        ccTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAATGTCTATTCTGTGATAAATGGCGCTGATTCATa                             >  1:6947974/1‑150 (MQ=255)
                                                                                                                        ccTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATa                             <  1:5023637/150‑1 (MQ=255)
                                                                                                                        ccTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATa                             <  2:7614665/150‑1 (MQ=255)
                                                                                                                         cTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAg                            >  1:7853893/1‑150 (MQ=255)
                                                                                                                         cTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAg                            <  1:9648871/150‑1 (MQ=255)
                                                                                                                         cTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTAACGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAg                            <  2:8112300/150‑1 (MQ=255)
                                                                                                                             cTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGTTTATATTGAGGCAAGTCACCCTCCCTTTTATTAAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTa                        <  1:2498956/150‑1 (MQ=255)
                                                                                                                              tCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTaa                       >  1:2305247/1‑150 (MQ=255)
                                                                                                                                aaTTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTaaaa                     <  2:9222638/150‑1 (MQ=255)
                                                                                                                                  ttAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGAGTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCACAGCTTaaaaaa                   >  2:7351556/1‑150 (MQ=255)
                                                                                                                                   tAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAAt                  <  2:5676595/150‑1 (MQ=255)
                                                                                                                                    aGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATa                 <  2:6146911/150‑1 (MQ=255)
                                                                                                                                    aGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATa                 >  2:8405802/1‑150 (MQ=255)
                                                                                                                                     gCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATAc                <  2:1237518/150‑1 (MQ=255)
                                                                                                                                     gCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATAc                <  1:1571463/150‑1 (MQ=255)
                                                                                                                                      cTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGATGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATAcc               >  1:4627253/1‑150 (MQ=255)
                                                                                                                                      cTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATAcc               >  2:6536754/1‑150 (MQ=255)
                                                                                                                                      cTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATAcc               <  1:413221/150‑1 (MQ=255)
                                                                                                                                      cTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATAcc               >  1:724182/1‑150 (MQ=255)
                                                                                                                                        aTTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCt             >  1:1244764/1‑150 (MQ=255)
                                                                                                                                           aaTCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAATTCACCCTCACATTAATTCAGAGTAAAAGTCTATACTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATaacatagt          >  2:2228730/1‑143 (MQ=255)
                                                                                                                                           aaTCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGt          <  2:8299427/150‑1 (MQ=255)
                                                                                                                                            atcatcGCCCGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGtc         >  2:4465450/1‑150 (MQ=255)
                                                                                                                                            atcatcGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGtc         <  1:5287940/150‑1 (MQ=255)
                                                                                                                                            atcatcGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGtc         <  1:6801322/150‑1 (MQ=255)
                                                                                                                                             tcatcGCGCGCTGGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGATTAAAAAATACCCTTGtca        >  1:8775094/1‑150 (MQ=255)
                                                                                                                                             tcatcGCCCGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGtca        >  2:4282743/1‑150 (MQ=255)
                                                                                                                                              catcGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGtcaa       <  2:9010657/150‑1 (MQ=255)
                                                                                                                                                tcGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGtcaatc     >  2:1860998/1‑150 (MQ=255)
                                                                                                                                                   ccAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGTCAATCAAc  <  2:7030426/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCTGAGCGATTTTGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGTCAATCAAC  >  NC_000913/1225966‑1226262

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