Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,226,115 A→C 6.6% 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 6.6% 677.7 / 23.6 224intergenic (‑35/‑485)minC/ycgJinhibitor of FtsZ ring polymerization/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (102/105);  new base C (12/4);  total (115/109)
Fisher's exact test for biased strand distribution p-value = 6.75e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 228 total aligned reads displayed.

GCGATTTTGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGTCAATCAACCC  >  NC_000913/1225971‑1226264
                                                                                                                                                |                                                                                                                                                     
      ttGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCgatga                                                                                                                                            >  1:1002245/1‑150 (MQ=255)
       tGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCgatgat                                                                                                                                           <  2:898557/150‑1 (MQ=255)
         tCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGt                                                                                                                                         >  1:3735979/1‑150 (MQ=255)
         tCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGt                                                                                                                                         >  1:3899183/1‑150 (MQ=255)
           ttCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTc                                                                                                                                       >  2:3685865/1‑150 (MQ=255)
           ttCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTc                                                                                                                                       <  2:695238/150‑1 (MQ=255)
            tCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTcc                                                                                                                                      <  1:1933415/150‑1 (MQ=255)
             ccAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCg                                                                                                                                     >  2:4097081/1‑150 (MQ=255)
              cAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGa                                                                                                                                    >  1:1225633/1‑150 (MQ=255)
               aGCGCCTGATGGATAACCTGCGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGaa                                                                                                                                   <  2:3405194/150‑1 (MQ=255)
                 cgcCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGa                                                                                                                                 >  2:613367/1‑150 (MQ=255)
                   ccTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACt                                                                                                                               >  1:2071243/1‑150 (MQ=255)
                      gATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACtata                                                                                                                            >  2:1247488/1‑150 (MQ=255)
                       aTGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCGGCGCGCGATGATGTTCCGAAGACTATaa                                                                                                                           >  2:2673216/1‑150 (MQ=255)
                       aTGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATaa                                                                                                                           >  1:2443793/1‑150 (MQ=255)
                       aTGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATaa                                                                                                                           <  2:2759351/150‑1 (MQ=255)
                           aTAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCa                                                                                                                       <  1:2880552/150‑1 (MQ=255)
                                 ttAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTtat                                                                                                                 >  2:2815643/1‑150 (MQ=255)
                                  tAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTtata                                                                                                                <  2:668969/150‑1 (MQ=255)
                                      ttCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCt                                                                                                            >  1:1173069/1‑150 (MQ=255)
                                       tCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTg                                                                                                           >  2:489317/1‑150 (MQ=255)
                                       tCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTg                                                                                                           >  2:3403923/1‑150 (MQ=255)
                                          gCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCCGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGAt                                                                                                        >  1:1135747/1‑150 (MQ=255)
                                          gCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGAt                                                                                                        <  2:2307727/150‑1 (MQ=255)
                                          gCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGAt                                                                                                        >  1:3202526/1‑150 (MQ=255)
                                           ccTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGAtt                                                                                                       >  1:3035558/1‑150 (MQ=255)
                                           ccTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGAtt                                                                                                       >  1:3499609/1‑150 (MQ=255)
                                           ccTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGAtt                                                                                                       >  2:813165/1‑150 (MQ=255)
                                            cTCATGCAGATGAACCACAGATAAAGTGAAGCTTCTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATta                                                                                                      <  2:2474110/150‑1 (MQ=255)
                                             tCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATtat                                                                                                     >  2:3906334/1‑150 (MQ=255)
                                             tCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATtat                                                                                                     >  2:2958662/1‑150 (MQ=255)
                                              cATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATtata                                                                                                    <  1:523422/150‑1 (MQ=255)
                                              cATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATtata                                                                                                    <  2:1067632/150‑1 (MQ=255)
                                                tGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATAtt                                                                                                  <  2:3295732/150‑1 (MQ=255)
                                                 gCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTg                                                                                                 >  2:3911488/1‑150 (MQ=255)
                                                  cAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGa                                                                                                <  2:2440251/150‑1 (MQ=255)
                                                   aGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAg                                                                                               <  1:4001127/150‑1 (MQ=255)
                                                    gATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAgg                                                                                              <  1:1748081/150‑1 (MQ=255)
                                                            acaGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCAc                                                                                      >  1:333901/1‑150 (MQ=255)
                                                            acaGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCAc                                                                                      >  2:3076888/1‑150 (MQ=255)
                                                            acaGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTAGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCAc                                                                                      >  2:3962912/1‑150 (MQ=255)
                                                                aTAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTc                                                                                  >  2:3913697/1‑150 (MQ=255)
                                                                aTAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTc                                                                                  <  2:3759320/150‑1 (MQ=255)
                                                                       gAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTa                                                                           >  2:3471017/1‑150 (MQ=255)
                                                                        aaGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTAt                                                                          >  2:1953435/1‑150 (MQ=255)
                                                                         aGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTAtt                                                                         <  1:3274052/150‑1 (MQ=255)
                                                                          gCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTc                                                                        <  1:3065538/150‑1 (MQ=255)
                                                                           cTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCa                                                                       <  1:2919584/150‑1 (MQ=255)
                                                                           cTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCa                                                                       <  2:1872494/150‑1 (MQ=255)
                                                                           cTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCa                                                                       <  1:511201/150‑1 (MQ=255)
                                                                             aCTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCaga                                                                     <  2:3187462/150‑1 (MQ=255)
                                                                             aCTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCaga                                                                     <  2:1183610/150‑1 (MQ=255)
                                                                             aCTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCaga                                                                     <  1:1980081/150‑1 (MQ=255)
                                                                               tGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGt                                                                   >  1:2173959/1‑150 (MQ=255)
                                                                               tGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGt                                                                   >  1:2802665/1‑150 (MQ=255)
                                                                               tGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGt                                                                   >  2:2541806/1‑150 (MQ=255)
                                                                                 ccTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTaa                                                                 >  2:1476263/1‑150 (MQ=255)
                                                                                 ccTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTaa                                                                 >  1:3922384/1‑150 (MQ=255)
                                                                                     tAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGt                                                             <  1:4099657/150‑1 (MQ=255)
                                                                                             aTTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCtg                                                     <  2:2443793/150‑1 (MQ=255)
                                                                                              ttGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCtgt                                                    <  2:3900255/150‑1 (MQ=255)
                                                                                                ggCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTTAAAGTCTATTCTGTGa                                                  <  1:378841/150‑1 (MQ=255)
                                                                                                ggCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGa                                                  >  2:7133/1‑150 (MQ=255)
                                                                                                ggCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGa                                                  <  1:1245760/150‑1 (MQ=255)
                                                                                                ggCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGa                                                  <  1:2335971/150‑1 (MQ=255)
                                                                                                    tgtTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATaaa                                              <  2:1284732/150‑1 (MQ=255)
                                                                                                      tttGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATg                                            >  1:1868609/1‑150 (MQ=255)
                                                                                                        tGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGgc                                          <  1:569264/150‑1 (MQ=255)
                                                                                                         gACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGgcg                                         <  2:3202526/150‑1 (MQ=255)
                                                                                                          aCATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGgcgc                                        >  1:1162519/1‑150 (MQ=255)
                                                                                                             tCCTTGCATTACTCAATTAGCTATTAATAATAGCAAGCGCGCGATGTTGTTCCGAATAGTATAAGGCATGTTATAGTCTGGATTATATTGAGGGAAGTCACCCTCCCATTTATTCAGAGTAAAAGTGTATTGTGTGATAAATGGGGCTgt                                     >  2:1519337/1‑149 (MQ=255)
                                                                                                             tCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGa                                     <  1:3471017/150‑1 (MQ=255)
                                                                                                               cTGGCCTTACTCAATTAGCTATTAATCATCGCCCGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATCTTGTGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGAtt                                   >  2:2612780/1‑150 (MQ=255)
                                                                                                               cTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGAtt                                   >  2:2919008/1‑150 (MQ=255)
                                                                                                                  gCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCAt                                <  1:1534604/150‑1 (MQ=255)
                                                                                                                  gCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCAt                                >  2:2648484/1‑150 (MQ=255)
                                                                                                                  gCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCAt                                <  2:2071243/150‑1 (MQ=255)
                                                                                                                    cTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTCTTCTGTGATAAATGGCGCTGATTCATAg                              >  2:2836224/1‑150 (MQ=255)
                                                                                                                     ttACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGc                             >  1:934838/1‑150 (MQ=255)
                                                                                                                     ttACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGc                             >  1:2762860/1‑150 (MQ=255)
                                                                                                                     ttACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGc                             >  1:2127957/1‑150 (MQ=255)
                                                                                                                         tCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTaa                         >  1:1314036/1‑150 (MQ=255)
                                                                                                                               accTAATAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAACGCATGTTATAGTCTGGAATATATTGAGGCAAGTCACCCTCCCATATATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATa                   <  1:2673216/148‑1 (MQ=255)
                                                                                                                               aGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATa                   <  1:3750908/150‑1 (MQ=255)
                                                                                                                               aGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATa                   <  1:1088046/150‑1 (MQ=255)
                                                                                                                                 cTATTAATCATCGCCCGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATAcc                 >  1:3430869/1‑150 (MQ=255)
                                                                                                                                  tATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATAccc                >  2:936164/1‑150 (MQ=255)
                                                                                                                                  tATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATAccc                <  2:2733112/150‑1 (MQ=255)
                                                                                                                                  tATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATAccc                >  2:1688312/1‑150 (MQ=255)
                                                                                                                                  tATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATAccc                <  2:1002245/150‑1 (MQ=255)
                                                                                                                                   aTTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTCAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCt               >  2:2092389/1‑150 (MQ=255)
                                                                                                                                    aacatcatcGCGCGCTGGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCtt              >  1:2240804/4‑150 (MQ=255)
                                                                                                                                      aaTCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGt            <  2:1768085/150‑1 (MQ=255)
                                                                                                                                      aaTCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTATGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGt            <  2:3035827/150‑1 (MQ=255)
                                                                                                                                       atcatcGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGtc           >  2:1351356/1‑150 (MQ=255)
                                                                                                                                        tcatcGCGCGCTGGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGtca          >  1:1417767/1‑150 (MQ=255)
                                                                                                                                             gCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGtcaatcaa     >  2:3843295/1‑150 (MQ=255)
                                                                                                                                              ccAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGGTTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGTCAATCAAc    >  2:1047717/1‑150 (MQ=255)
                                                                                                                                              ccAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGTCAATCAAc    >  2:1546981/1‑150 (MQ=255)
                                                                                                                                                aGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGTCAATCAAccc  >  2:3504169/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
GCGATTTTGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGTCAATCAACCC  >  NC_000913/1225971‑1226264

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: