Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,042,411 C→A 5.7% intergenic (+118/+78) ygfZ → / ← yqfA iron‑sulfur cluster repair protein, plumbagin resistance/hemolysin III family HylIII inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,042,4110CA 5.7% 709.8 / 24.0 193intergenic (+118/+78)ygfZ/yqfAiron‑sulfur cluster repair protein, plumbagin resistance/hemolysin III family HylIII inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (90/92);  new base A (3/8);  total (93/100)
Fisher's exact test for biased strand distribution p-value = 2.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

Only 100 of 193 total aligned reads displayed.

GAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGC  >  NC_000913/3042264‑3042558
                                                                                                                                                   |                                                                                                                                                   
gAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCGCTTTTTGGCATCAGGTATTTATGGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCAc                                                                                                                                                   <  1:1739137/150‑1 (MQ=255)
  gCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGc                                                                                                                                                 <  1:2157407/150‑1 (MQ=255)
  gCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGc                                                                                                                                                 <  2:1183065/150‑1 (MQ=255)
         ccGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATcc                                                                                                                                          >  2:1990168/1‑150 (MQ=255)
           gTATTCGCTCGAATAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATTCGGCGTGAACGCCTTATCCga                                                                                                                                        <  1:2184106/150‑2 (MQ=38)
           gTATTCGCTCGAAGTGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGTTGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATTCGGCGTGAACGCCTTATCCga                                                                                                                                        <  1:2184311/150‑2 (MQ=255)
           gTATTCGCTCGAAGATTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATTCGGCGTGAACGCCTTATCCga                                                                                                                                        <  1:650346/150‑2 (MQ=38)
           gTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATTCGGCGTGAACGCCTTATCCga                                                                                                                                        <  2:1363988/150‑2 (MQ=39)
           gTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATTCGGCGTGAACGCCTTATCCga                                                                                                                                        <  2:1675939/150‑2 (MQ=39)
           gTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATTCGGCGTGAACGCCTTATCCga                                                                                                                                        <  1:2326964/150‑2 (MQ=39)
                cGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGtt                                                                                                                                   <  2:222248/150‑1 (MQ=255)
                 gCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTtg                                                                                                                                  >  2:2055395/1‑150 (MQ=255)
                   tCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTtgtg                                                                                                                                >  1:2388106/1‑150 (MQ=255)
                     gAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCt                                                                                                                              <  2:1080943/150‑1 (MQ=255)
                     gAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCt                                                                                                                              >  2:798678/1‑150 (MQ=255)
                       agagTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTcc                                                                                                                            <  1:580966/150‑1 (MQ=255)
                             aaTCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGcc                                                                                                                      <  2:2149727/150‑1 (MQ=255)
                             aaTCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGcc                                                                                                                      <  2:1321985/150‑1 (MQ=255)
                              ctctcTTTATCGCATCAGGCATTTATCGCCGGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCt                                                                                                                     <  2:219156/149‑1 (MQ=255)
                              aTCTCTTTATCGCATCAGGCATTTATCGCCTGAGGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCt                                                                                                                     <  1:504909/150‑1 (MQ=255)
                               tctcTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTg                                                                                                                    >  2:358818/1‑150 (MQ=255)
                               tctcTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTg                                                                                                                    >  2:358774/1‑150 (MQ=255)
                               tctcTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTg                                                                                                                    >  2:2115463/1‑150 (MQ=255)
                                 tctTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGAt                                                                                                                  >  1:1880525/1‑150 (MQ=255)
                                  ctTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCt                                                                                                                     <  2:682297/146‑1 (MQ=255)
                                  ctTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATg                                                                                                                 <  2:1970332/150‑1 (MQ=255)
                                         gCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCt                                                                                                          <  2:1547918/150‑1 (MQ=255)
                                            tCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGgc                                                                                                       >  1:2010653/1‑150 (MQ=255)
                                              aGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGgcgc                                                                                                     >  2:2257734/1‑150 (MQ=255)
                                              aGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGgcgc                                                                                                     <  1:1224723/150‑1 (MQ=255)
                                              aGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGgcgc                                                                                                     <  1:552539/150‑1 (MQ=255)
                                              aGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGgcgc                                                                                                     <  1:84656/150‑1 (MQ=255)
                                               ggCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGcgcg                                                                                                    >  1:1479714/1‑150 (MQ=255)
                                                 cATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTc                                                                                                  <  2:920202/150‑1 (MQ=255)
                                                     tATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                              <  1:1548040/150‑1 (MQ=255)
                                                     tATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTGATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                              <  2:1775557/150‑1 (MQ=255)
                                                      aTCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATc                                                                                             >  2:157556/1‑150 (MQ=255)
                                                       tCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCa                                                                                            >  1:1537611/1‑150 (MQ=255)
                                                       tCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCa                                                                                            >  1:1662076/1‑150 (MQ=255)
                                                       tCGCCTGATGCGACGCTGGCGCGTCATATCAAGTATGGGATTTGATGCCTTTATAAAAAAAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAAAACGGCAGTTGTGCTCCGATTCCAGATTCTACGCTGGCGCGTCTTAtaa                                                                                            >  1:403756/1‑148 (MQ=37)
                                                        cGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAg                                                                                           <  1:1869570/150‑1 (MQ=255)
                                                          ccTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCt                                                                                                          >  1:585759/1‑133 (MQ=255)
                                                          ccTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCt                                                                                                          <  2:585759/133‑1 (MQ=255)
                                                          ccTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGc                                                                                         >  2:297064/1‑150 (MQ=255)
                                                          ccTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGc                                                                                         >  2:1952508/1‑150 (MQ=255)
                                                          ccTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGc                                                                                         <  2:1618088/150‑1 (MQ=255)
                                                          ccTGAGGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGc                                                                                         <  1:2309540/150‑1 (MQ=255)
                                                                gCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACaa                                                                                   >  2:528857/1‑150 (MQ=255)
                                                                gCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACaa                                                                                   >  2:1091628/1‑150 (MQ=255)
                                                                 cGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACaaa                                                                                  <  2:1527200/150‑1 (MQ=255)
                                                                 cGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACaaa                                                                                  <  1:1091628/150‑1 (MQ=255)
                                                                  gACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAg                                                                                 >  1:2337758/1‑150 (MQ=255)
                                                                  gACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAg                                                                                 >  1:2338087/1‑150 (MQ=255)
                                                                     gCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATCAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCa                                                                              <  2:2388106/150‑1 (MQ=255)
                                                                     gCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCa                                                                              >  2:1933717/1‑150 (MQ=255)
                                                                      cTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCAt                                                                             >  2:2291282/1‑150 (MQ=255)
                                                                         gcgcGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATAcc                                                                          >  1:1043872/1‑150 (MQ=255)
                                                                         gcgcGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATAcc                                                                          >  1:76008/1‑150 (MQ=255)
                                                                           gcgTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCa                                                                        >  1:1529167/1‑150 (MQ=255)
                                                                             gTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCAtt                                                                      >  1:954256/1‑150 (MQ=255)
                                                                             gTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCAtt                                                                      <  2:1118312/150‑1 (MQ=255)
                                                                             gTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCAtt                                                                      <  1:1629050/150‑1 (MQ=255)
                                                                                ttATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACg                                                                   >  1:383039/1‑150 (MQ=255)
                                                                                  aTCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGcc                                                                 >  1:1430442/1‑150 (MQ=255)
                                                                                   tCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCt                                                                >  2:474624/1‑150 (MQ=255)
                                                                                         cTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCaa                                                          <  2:383039/150‑1 (MQ=255)
                                                                                          tGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAat                                                         <  2:1289971/150‑1 (MQ=255)
                                                                                          tGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAat                                                         <  1:60038/150‑1 (MQ=255)
                                                                                            ggATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCActat                                                       >  1:1934086/1‑150 (MQ=255)
                                                                                            ggATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAatat                                                       >  1:1804210/1‑150 (MQ=255)
                                                                                              aTTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATAc                                                     >  1:720253/1‑150 (MQ=255)
                                                                                                   aTGCCTTTATCACAACAGGGACGAGGACCGTAGATCGGATAAGGCGTTCACGCCGAATACGGCAGTTGTGCTCCGATGCCTTATTCGACGCTGGCGCGTCTTATCAGGCCTACAAACGCATACCCATTTCGCCTGCCAAATTTACAAATa                                                <  2:1302607/150‑1 (MQ=37)
                                                                                                         ttATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGc                                          >  2:858010/1‑150 (MQ=255)
                                                                                                          tATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGcc                                         <  2:1430887/150‑1 (MQ=255)
                                                                                                              acaacaGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAATGGCATACCCATTACGCCTGCCCAATATTCTAATAGATCGCCAGaa                                     >  2:1969633/1‑150 (MQ=255)
                                                                                                              acaacaGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGaa                                     >  2:394511/1‑150 (MQ=255)
                                                                                                                     ggACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGca                              <  1:2257734/150‑1 (MQ=255)
                                                                                                                         gAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGcacaca                          <  1:358774/150‑1 (MQ=255)
                                                                                                                         gAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGcacaca                          <  1:65920/150‑1 (MQ=255)
                                                                                                                         gAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGcacaca                          <  1:2055395/150‑1 (MQ=255)
                                                                                                                         gAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGcacaca                          <  1:2115463/150‑1 (MQ=255)
                                                                                                                          aGGCCCGTAGATCGGATAAGTCGTTCACGCCGAATCCGGCAGTTGTGCTCCTATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCTCCAGAAAGTGGCacacac                         <  2:403756/150‑1 (MQ=255)
                                                                                                                          aGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCacacac                         <  2:1537611/150‑1 (MQ=255)
                                                                                                                            gCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTa                       >  2:830483/1‑150 (MQ=255)
                                                                                                                            gCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTa                       <  2:727187/150‑1 (MQ=255)
                                                                                                                            gCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTa                       >  1:560224/1‑150 (MQ=255)
                                                                                                                                gTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTAccgc                   >  1:388857/1‑150 (MQ=255)
                                                                                                                                   tgccGGATTCGGCGTGAACGCCTTATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGa                >  2:1742404/4‑150 (MQ=37)
                                                                                                                                   gATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGa                <  1:798678/150‑1 (MQ=255)
                                                                                                                                    gccGGATTCGGCGTGAACGCCTTATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAg               >  2:818002/3‑150 (MQ=38)
                                                                                                                                    aTCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTAcc                     >  1:2240338/1‑144 (MQ=255)
                                                                                                                                    aTCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAg               >  1:559583/1‑150 (MQ=255)
                                                                                                                                    aTCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAg               >  1:914935/1‑150 (MQ=255)
                                                                                                                                     tCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGc              >  2:1877993/1‑150 (MQ=255)
                                                                                                                                     tCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGc              >  2:2025970/1‑150 (MQ=255)
                                                                                                                                            aGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGc       >  1:1085538/1‑150 (MQ=255)
                                                                                                                                              gCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCg     >  1:1463150/1‑150 (MQ=255)
                                                                                                                                              gCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCg     >  1:1296088/1‑150 (MQ=255)
                                                                                                                                              gCGTTCACACCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCg     >  2:2328827/1‑150 (MQ=255)
                                                                                                                                                 ttCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGc  <  1:276557/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GAGCCGCTGCCGTATTCGCTCGAAGAGTAAATCTCTTTATCGCATCAGGCATTTATCGCCTGATGCGACGCTGGCGCGTCTTATCATGTCTGGGATTTGATGCCTTTATCACAACAGGGACGAGGCCCGTAGATCGGATAAGGCGTTCACGCCGAATCCGGCAGTTGTGCTCCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGC  >  NC_000913/3042264‑3042558

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