Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,226,115 A→C 6.7% 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.7% 271.1 / 7.0 99intergenic (‑35/‑485)minC/ycgJinhibitor of FtsZ ring polymerization/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (52/39);  new base C (5/2);  total (57/42)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

TGAGCGATTTTGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGTCAATCAACC  >  NC_000913/1225968‑1226263
                                                                                                                                                   |                                                                                                                                                    
tGAGCGATTTTGTCTTCCAGCGCCTGATGGATACCGTTAGGTTCTGCCTCAGGCAGATGAACCACAGTTAACGTGAAGCTAATGCATTTAAGCACGATTGGCGGTTTTGACAGCGTGGCGTTAGGAAATTAGCTTTTAATCTGCGCCCgc                                                                                                                                                    >  2:1258477/1‑150 (MQ=255)
     gATTTTGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGcgcgcg                                                                                                                                               >  1:601930/1‑150 (MQ=255)
     gATTTTGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGcgcgcg                                                                                                                                               >  1:1679805/1‑150 (MQ=255)
       ttttGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCgat                                                                                                                                             <  1:1010824/150‑1 (MQ=255)
          tGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCGCGCTGGCgatgat                                                                                                                                          <  1:736093/150‑1 (MQ=255)
           gTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCGCGCTGGCgatgatt                                                                                                                                         <  2:720191/150‑2 (MQ=255)
             cTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGtt                                                                                                                                       >  2:406129/1‑150 (MQ=255)
             cTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGAGCGCGATGATGtt                                                                                                                                       >  2:1638632/1‑150 (MQ=255)
              ttCCAGCGCCTTATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTc                                                                                                                                      <  1:623967/150‑1 (MQ=255)
              ttCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTc                                                                                                                                      <  2:466192/150‑1 (MQ=255)
              ttCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTc                                                                                                                                      <  1:1734107/150‑1 (MQ=255)
               tCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTcc                                                                                                                                     >  2:1029486/1‑150 (MQ=255)
               tCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTGTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTcc                                                                                                                                     <  2:294063/150‑1 (MQ=255)
                ccAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGAGTGTTTGACATCATGTCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCg                                                                                                                                    >  1:608516/1‑150 (MQ=255)
                 cAGCGCCTCATTGAAAACCTTCGGTTCTGCCTCATCCAGATGAACCACAGATAAACTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTGGACATCCTGGCCTTACTCAATTCGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGa                                                                                                                                   <  1:825667/150‑1 (MQ=255)
                  aGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGaa                                                                                                                                  >  2:496888/1‑150 (MQ=255)
                  aGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGaa                                                                                                                                  >  1:454606/1‑150 (MQ=255)
                  aGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAAACACACATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATACTGGCCTAACTCAATTAGCTATTAATCATCGCCAGCGCGAGATGATGTTCCGaa                                                                                                                                  >  2:460082/1‑150 (MQ=255)
                     gcCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGAc                                                                                                                               <  1:1611800/150‑1 (MQ=255)
                     gcCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGAc                                                                                                                               <  1:1254655/150‑1 (MQ=255)
                       cTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACGGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACta                                                                                                                             >  2:1137368/1‑150 (MQ=255)
                        tGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACtat                                                                                                                            >  2:1723469/1‑150 (MQ=255)
                        tGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACtat                                                                                                                            >  1:651935/1‑150 (MQ=255)
                          aTGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACGGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATaa                                                                                                                          <  1:1137368/150‑1 (MQ=255)
                          aTGGATAACCTTAGGTTCTGCCTCATGCAGATGAACAACAGATAAAGTTATGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATaa                                                                                                                          <  2:236685/150‑1 (MQ=255)
                           tGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAg                                                                                                                         <  1:594160/150‑1 (MQ=255)
                             gATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGc                                                                                                                       >  1:1301526/1‑150 (MQ=255)
                              aTAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCa                                                                                                                      <  1:1486043/150‑1 (MQ=255)
                              aTAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCa                                                                                                                      <  1:1406008/150‑1 (MQ=255)
                                 aCCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGt                                                                                                                   >  1:606931/1‑150 (MQ=255)
                                 aCCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGATTGACATCATGTCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGt                                                                                                                   >  1:1551483/1‑150 (MQ=255)
                                  ccTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGtt                                                                                                                  >  2:871478/1‑150 (MQ=255)
                                   cTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTta                                                                                                                 >  1:353125/1‑150 (MQ=255)
                                             gCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGAt                                                                                                       >  2:987969/1‑150 (MQ=255)
                                             gCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGAt                                                                                                       >  1:1582755/1‑150 (MQ=255)
                                             gCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGAt                                                                                                       >  1:58773/1‑150 (MQ=255)
                                             gCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGAt                                                                                                       >  2:1718902/1‑150 (MQ=255)
                                                 cATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATtata                                                                                                   >  1:704371/1‑150 (MQ=255)
                                                   tGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTTATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATAtt                                                                                                 <  1:1550390/150‑1 (MQ=255)
                                                     cAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGAAGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGa                                                                                               <  2:601930/150‑1 (MQ=255)
                                                      aGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAg                                                                                              <  2:1413490/150‑1 (MQ=255)
                                                      aGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAg                                                                                              >  1:491363/1‑150 (MQ=255)
                                                       gATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAgg                                                                                             <  2:39495/150‑1 (MQ=255)
                                                            aCCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGt                                                                                        >  2:430514/1‑150 (MQ=255)
                                                               acaGATAAAGTGAAGCTACTGCCTTTACGCTCGATTGGCGTTTTTGACATCCTGGCCTTACTCAATTAGCTATTAATAATCGCCAGCGCGCGATGATGTTCCGAAGACTATAATGCATGTTATAGTCTGGATTATATTGAGGCAAGTCAc                                                                                     >  1:85641/1‑150 (MQ=255)
                                                               acaGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCAc                                                                                     <  2:454606/150‑1 (MQ=255)
                                                                  gATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCATGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACACt                                                                                  >  1:1020764/1‑150 (MQ=255)
                                                                   aTAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTc                                                                                 <  1:1666728/150‑1 (MQ=255)
                                                                    tAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGAGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCTGCGCGCGATGATGTTCCGAAGACTTTAAGGCATGTTATAGTCTGGATTATTTTGAGGCAAGTCACCCTcc                                                                                <  2:755650/150‑1 (MQ=255)
                                                                      aaGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCa                                                                              >  1:406028/1‑150 (MQ=255)
                                                                      aaGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCa                                                                              >  1:1672949/1‑150 (MQ=255)
                                                                       aGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTTTAGTCTGGATTATATTGAGGCAAGTCGCCCTCCCAt                                                                             >  1:868321/1‑150 (MQ=255)
                                                                       aGTGAAGCAACTGCCTTTAAGCTCGATGGGCGTGTATGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGGCTGGATTATATTGCGGCAAGTCCCCCTCCCAt                                                                             >  2:1210522/1‑150 (MQ=255)
                                                                         tGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCAttt                                                                           <  2:754508/150‑1 (MQ=255)
                                                                          gAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTa                                                                          >  1:202233/1‑150 (MQ=255)
                                                                           aaGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTAt                                                                         >  2:1625885/1‑150 (MQ=255)
                                                                            aGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTAtt                                                                        >  1:1402223/1‑150 (MQ=255)
                                                                             gCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTc                                                                       <  2:1413438/150‑1 (MQ=255)
                                                                             gCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTc                                                                       <  1:756822/150‑1 (MQ=255)
                                                                               tACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCATTTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCag                                                                     <  2:608516/150‑1 (MQ=255)
                                                                               tACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGATCCGAAGACTATAAGGCAGGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCag                                                                     <  1:195466/150‑1 (MQ=255)
                                                                                aCTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCaga                                                                    >  2:344462/1‑150 (MQ=255)
                                                                                 cTGCCTTTACGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCagag                                                                   <  2:828168/150‑1 (MQ=255)
                                                                                 cTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCagag                                                                   <  2:295409/150‑1 (MQ=255)
                                                                                     cTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTaaa                                                               >  2:1366417/1‑150 (MQ=255)
                                                                                        tAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAATTCACCCTCCCATTTATTCAGAGTAAAAGt                                                            >  1:1242946/1‑150 (MQ=255)
                                                                                         aaGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTc                                                           >  1:731272/1‑150 (MQ=255)
                                                                                          aGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCt                                                          >  1:860227/1‑150 (MQ=255)
                                                                                                   ggCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGAt                                                >  1:946840/1‑150 (MQ=255)
                                                                                                      gtgtTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTATGTGATaa                                              >  1:955154/1‑150 (MQ=255)
                                                                                                       tgtTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATaaa                                             >  1:354737/1‑150 (MQ=255)
                                                                                                         tttGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATg                                           <  1:4379/150‑1 (MQ=255)
                                                                                                           tGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGgc                                         <  2:606931/150‑1 (MQ=255)
                                                                                                           tGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGgc                                         <  1:146604/150‑1 (MQ=255)
                                                                                                            gACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGgcg                                        <  2:873091/150‑1 (MQ=255)
                                                                                                             aCATCCTGGCCTTACTCAATTAGCAATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGgcgc                                       <  1:1546091/150‑1 (MQ=255)
                                                                                                              cATCCTGGCCTTACTCAATAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTg                                     <  2:20111/150‑1 (MQ=255)
                                                                                                                 ccTGGCCTTACTCAATTATATATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGAt                                   <  1:708888/150‑1 (MQ=255)
                                                                                                                         tACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCt                           >  1:1418390/1‑150 (MQ=255)
                                                                                                                         tACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCt                           >  1:1336348/1‑150 (MQ=255)
                                                                                                                           cTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTa                         <  2:58773/150‑1 (MQ=255)
                                                                                                                             cAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTaaa                       >  2:421086/1‑150 (MQ=255)
                                                                                                                             cAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTACCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTaaa                       >  2:40604/1‑150 (MQ=255)
                                                                                                                               aTTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTaaaaa                     >  2:1444288/1‑150 (MQ=255)
                                                                                                                                   gCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATAc                 <  1:496888/150‑1 (MQ=255)
                                                                                                                                    cggaacatcatcGCGCGCTGGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATAcc                >  2:669796/7‑150 (MQ=255)
                                                                                                                                     tATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATAccc               >  2:576283/1‑150 (MQ=255)
                                                                                                                                        tAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTg            >  2:944529/1‑150 (MQ=255)
                                                                                                                                        acatcatcGCGCGCTGGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTCTATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTg            >  2:1340615/3‑150 (MQ=255)
                                                                                                                                         catcatcGCGCGCTGGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTATATTCTGTGATAAATGTCGCTGATTAATATCTTAAAAAATAACCTTGt           >  1:324503/2‑150 (MQ=255)
                                                                                                                                         aaTCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGt           <  2:1679805/150‑1 (MQ=255)
                                                                                                                                          atcatcGCGCGCTGGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGta          >  2:1690302/1‑149 (MQ=255)
                                                                                                                                          atcatcGCCAGCGCGCGATGATGTTACGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGtc          <  1:611750/150‑1 (MQ=255)
                                                                                                                                             atcGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGtcaat       >  1:921468/1‑150 (MQ=255)
                                                                                                                                             atcGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGtcaat       >  1:794020/1‑150 (MQ=255)
                                                                                                                                             atcGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGtcaat       >  2:693175/1‑150 (MQ=255)
                                                                                                                                                gCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGtcaatcaa    <  1:1029486/150‑1 (MQ=255)
                                                                                                                                                 ccAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATAACCTTTTCAATAAAc   >  1:918692/1‑150 (MQ=255)
                                                                                                                                                  cAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGTCAATCAAcc  <  2:353125/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGAGCGATTTTGTCTTCCAGCGCCTGATGGATAACCTTAGGTTCTGCCTCATGCAGATGAACCACAGATAAAGTGAAGCTACTGCCTTTAAGCTCGATTGGCGTGTTTGACATCCTGGCCTTACTCAATTAGCTATTAATCATCGCCAGCGCGCGATGATGTTCCGAAGACTATAAGGCATGTTATAGTCTGGATTATATTGAGGCAAGTCACCCTCCCATTTATTCAGAGTAAAAGTCTATTCTGTGATAAATGGCGCTGATTCATAGCTTAAAAAATACCCTTGTCAATCAACC  >  NC_000913/1225968‑1226263

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