Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,134,675 G→T 7.5% intergenic (+118/‑54) flgF → / → flgG flagellar basal‑body rod protein FlgF/flagellar basal‑body rod protein FlgG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,134,6750GT 7.5% 304.2 / 10.3 93intergenic (+118/‑54)flgF/flgGflagellar basal‑body rod protein FlgF/flagellar basal‑body rod protein FlgG
Reads supporting (aligned to +/- strand):  ref base G (44/42);  new base T (2/5);  total (46/47)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.67e-01

GCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCCAACAACCTGGCAAACGTCAGTACTA  >  NC_000913/1134526‑1134822
                                                                                                                                                     |                                                                                                                                                   
gCCTTGCCAACCAACTGCTGTCGTTGAGTTAATTGTAAGGATACATGACAAGTATAAGTTGCCCGATGCTCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGATCACTTGTAGGCCGGATAAGGCGTTTAcgccg                                                                                                                                                     <  1:1075164/150‑1 (MQ=255)
gCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTAcgccg                                                                                                                                                     >  2:129377/1‑150 (MQ=255)
     gCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                                <  1:512403/150‑1 (MQ=255)
      ccaaccaaCTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCg                                                                                                                                               >  1:1815412/1‑150 (MQ=255)
         accaacTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                            <  2:954187/150‑2 (MQ=255)
         accaacTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa                                                                                                                                            <  2:66890/150‑1 (MQ=255)
           caacTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCaag                                                                                                                                          <  2:1800588/150‑1 (MQ=255)
              ctgctgTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCaagaag                                                                                                                                       <  2:923942/150‑1 (MQ=255)
               tgctgTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGcctacaag                                                                                                                                      <  2:1307039/150‑8 (MQ=38)
                gctgTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCaag                                                                                                                                          <  2:523160/145‑1 (MQ=255)
                gctgTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCaag                                                                                                                                          >  1:523160/1‑145 (MQ=255)
                 ctgTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACa                                                                                                                                    >  1:647522/1‑150 (MQ=255)
                  tgTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACat                                                                                                                                   <  2:1183597/150‑1 (MQ=255)
                     cGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATg                                                                                                                                >  1:1596500/1‑150 (MQ=255)
                        tGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCAc                                                                                                                             <  2:1013143/150‑1 (MQ=255)
                        tGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCAc                                                                                                                             <  2:647522/150‑1 (MQ=255)
                        tGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCAc                                                                                                                             <  1:1828677/150‑1 (MQ=255)
                           gTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACttt                                                                                                                          >  2:242227/1‑150 (MQ=255)
                           gTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACttt                                                                                                                          >  2:581643/1‑150 (MQ=255)
                              aaTTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTc                                                                                                                       <  2:977767/150‑1 (MQ=255)
                              aaTTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTc                                                                                                                       <  1:914100/150‑1 (MQ=255)
                              aaTTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTc                                                                                                                       >  1:1285057/1‑150 (MQ=255)
                              aaTTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTc                                                                                                                       >  1:1093390/1‑150 (MQ=255)
                                ttGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCAc                                                                                                                     <  2:383798/150‑1 (MQ=255)
                                ttGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCAc                                                                                                                     <  1:1545764/150‑1 (MQ=255)
                                ttGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCAc                                                                                                                     <  1:974319/150‑1 (MQ=255)
                                  gAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTa                                                                                                                   >  1:1330417/1‑150 (MQ=255)
                                  gAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTa                                                                                                                   >  2:861855/1‑150 (MQ=255)
                                           cATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTAc                                                                                                          <  2:1643517/150‑1 (MQ=255)
                                           cATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTAc                                                                                                          <  1:1754583/150‑1 (MQ=255)
                                            aTGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACa                                                                                                         <  2:167314/150‑1 (MQ=255)
                                             tGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAg                                                                                                        >  2:1821193/1‑150 (MQ=255)
                                             tGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAg                                                                                                        >  1:81181/1‑150 (MQ=255)
                                             tGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAg                                                                                                        <  1:853086/150‑1 (MQ=255)
                                                cAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGAc                                                                                                                                     >  1:341181/1‑118 (MQ=255)
                                                cAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGAc                                                                                                                                     <  2:341181/118‑1 (MQ=255)
                                                      taAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTa                                                                                               >  2:830846/1‑150 (MQ=255)
                                                      taAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTa                                                                                               <  2:545665/150‑1 (MQ=255)
                                                      taAGTTGCCAGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGAGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATCCACTTTCTCACTAAACCACTACAGGACATTTTa                                                                                               >  2:416628/1‑150 (MQ=39)
                                                        atttGCCCGATGCGCAAGTTTATCGGGTCTATGTGGGCAATCTCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATg                                                                                             <  2:1285057/148‑1 (MQ=255)
                                                        aGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATg                                                                                             >  1:959760/1‑150 (MQ=255)
                                                        aGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATg                                                                                             <  1:1834447/150‑1 (MQ=255)
                                                            gCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCa                                                                                         >  2:1040284/1‑150 (MQ=255)
                                                             cccGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGGCATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAg                                                                                        >  1:1819420/1‑150 (MQ=255)
                                                              ccGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGt                                                                                       >  2:838642/1‑150 (MQ=255)
                                                               cGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGtt                                                                                      >  2:592004/1‑150 (MQ=255)
                                                                 aTGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCa                                                                                    <  2:111007/150‑1 (MQ=255)
                                                                   gcgcAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGGCATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCAtt                                                                                  <  2:1819420/150‑1 (MQ=255)
                                                                     gcAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTAt                                                                                >  1:619196/1‑150 (MQ=255)
                                                                     gcAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTAt                                                                                >  1:497441/1‑150 (MQ=255)
                                                                        aGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGa                                                                             <  2:1578052/150‑1 (MQ=255)
                                                                         gTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGAt                                                                            >  1:1535513/1‑150 (MQ=255)
                                                                          tttATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATc                                                                           <  1:1451154/150‑1 (MQ=255)
                                                                          tttATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATc                                                                           >  1:487676/1‑150 (MQ=255)
                                                                                 ggTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCaaaa                                                                    >  1:648610/1‑150 (MQ=255)
                                                                                  gTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAAc                                                                   <  1:1101819/150‑1 (MQ=255)
                                                                                     tATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACggg                                                                >  1:1800189/1‑150 (MQ=255)
                                                                                     tATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACggg                                                                >  2:122492/1‑150 (MQ=255)
                                                                                     tATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACggg                                                                >  2:896008/1‑150 (MQ=255)
                                                                                       tGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGcc                                                              <  2:1534955/150‑1 (MQ=255)
                                                                                       tGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGcc                                                              <  1:242227/150‑1 (MQ=255)
                                                                                        gggggCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACCGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCt                                                             >  2:694283/1‑150 (MQ=255)
                                                                                         ggggCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCtt                                                            >  1:1597772/1‑150 (MQ=255)
                                                                                          gggCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGATCATTATGGATCGCCAAAACGGGCCTTg                                                           >  1:1871884/1‑150 (MQ=255)
                                                                                                   cAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGc                                                  >  2:398568/1‑150 (MQ=255)
                                                                                                    aaTTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCa                                                 <  2:1462165/150‑1 (MQ=255)
                                                                                                     aTTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCaa                                                >  1:1547656/1‑150 (MQ=255)
                                                                                                      tttATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCaaa                                               <  2:959760/150‑1 (MQ=255)
                                                                                                       ttATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAAc                                              <  1:830846/150‑1 (MQ=255)
                                                                                                        tATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAAcc                                             <  1:1821193/150‑1 (MQ=255)
                                                                                                           cGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAat                                          >  1:760376/1‑150 (MQ=255)
                                                                                                                         cTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGc                            >  1:1120415/1‑150 (MQ=255)
                                                                                                                         cTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGc                            >  1:1459094/1‑150 (MQ=255)
                                                                                                                          ttGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGcc                           <  1:838642/150‑1 (MQ=255)
                                                                                                                           tgaaGGCCGAATAAGGCGTTTACGCCGCAACAGTCAAGAAGACATATGCACTTTGTAACTAATCCACAACAGGACAATTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCATGACGCCCAGCAAACCAATATGGACGTCATTCCca                          >  2:1630618/4‑150 (MQ=255)
                                                                                                                             tAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCcaac                        <  1:581643/150‑1 (MQ=255)
                                                                                                                               ggCCGGATAAGGCGATTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCcaacaa                      >  2:1966586/1‑150 (MQ=255)
                                                                                                                                gCCGGATGCGGCGTAAACGCCTTATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATgg                                                                              >  1:877879/1‑93 (MQ=37)
                                                                                                                                gCCGGATGCGGCGTAAACGCCTTATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATgg                                                                              <  2:884525/93‑1 (MQ=37)
                                                                                                                                gCCGGATGCGGCGTAAACGCCTTATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATgg                                                                              <  2:877879/93‑1 (MQ=37)
                                                                                                                                gCCGGATGCGGCGTAAACGCCTTATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATgg                                                                              >  1:884525/1‑93 (MQ=37)
                                                                                                                                  cGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCCAACAACCt                   <  2:619196/150‑1 (MQ=255)
                                                                                                                                  cGGATAAGGCGATTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCCAACAACCt                   <  1:1966586/150‑1 (MQ=255)
                                                                                                                                    gATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCCAACAACCTgg                 <  1:1040284/150‑1 (MQ=255)
                                                                                                                                     aTAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCCAACAACCTGGc                <  2:81181/150‑1 (MQ=255)
                                                                                                                                       aaGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCCAACAACCTGGCaa              >  2:439726/1‑150 (MQ=255)
                                                                                                                                        aGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCCAACAACCTGGCaaa             >  1:1618734/1‑150 (MQ=255)
                                                                                                                                        aGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCCAACAACCTGGCaaa             >  2:1395009/1‑150 (MQ=255)
                                                                                                                                        aGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCCAACAACCTGGCaaa             <  1:1664670/150‑1 (MQ=255)
                                                                                                                                         ggCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCCAACAACCTGGCAAAc            >  1:513186/1‑150 (MQ=255)
                                                                                                                                         ggCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCCAACAACCTGGCAAAc            >  1:809295/1‑150 (MQ=255)
                                                                                                                                         ggCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCCAACAACCTGGCAAAc            >  2:1535619/1‑150 (MQ=255)
                                                                                                                                            gTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCCAACAACCTGGCAAACGTc         <  1:592004/150‑1 (MQ=255)
                                                                                                                                              ttACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCCAACAACCTGGCAAACGTCAg       <  2:497441/150‑1 (MQ=255)
                                                                                                                                                aCGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCCAACAACCTGGCAAACGTCAGTa     <  1:122492/150‑1 (MQ=255)
                                                                                                                                                  gccgcATCCGGCAAGAAGACATATGCACTTTTTCACTAATCCACTACAGGACATTTTATGATCAGTTGATTATTTATCGCCAAAACTGGCCTTGACGCACAGCAAACCAATATGGTCGTCATTGCCAACAACCTTTCAAACGTAAGTACt   <  2:1815412/150‑1 (MQ=255)
                                                                                                                                                   ccgcATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAACAAACCAATATGGACGTCATTGCCAACAACCTGGCAAACGTCAGTACTa  <  2:1596500/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCCAACAACCTGGCAAACGTCAGTACTA  >  NC_000913/1134526‑1134822

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