Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,460,484 T→A 5.6% intergenic (‑120/+38) nrdG ← / ← nrdD anaerobic ribonucleotide reductase activating protein/anaerobic ribonucleoside‑triphosphate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,460,4840TA 5.6% 379.2 / 6.5 124intergenic (‑120/+38)nrdG/nrdDanaerobic ribonucleotide reductase activating protein/anaerobic ribonucleoside‑triphosphate reductase
Reads supporting (aligned to +/- strand):  ref base T (62/55);  new base A (2/5);  total (64/60)
Fisher's exact test for biased strand distribution p-value = 2.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

Only 100 of 129 total aligned reads displayed.

GTCGACAGGATAGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATCCGCACACGCGGCGA  >  NC_000913/4460335‑4460633
                                                                                                                                                     |                                                                                                                                                     
gTCGACAGGATAGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGttt                                                                                                                                                       <  2:980154/150‑1 (MQ=255)
gTCGACAGGATAGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGttt                                                                                                                                                       <  1:879514/150‑1 (MQ=255)
     cAGGATAGTACTGAGGATAATTCATTATTTCGCCTGGCTTATCAGACCGCTTCTAGGCTGCATACGGTATTCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAAGTGCCAAATGTCGGCCGGATAAGGAGATTAcgcc                                                                                                                                                  >  1:822797/1‑150 (MQ=38)
       ggATAGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTAcgccgc                                                                                                                                                >  2:947751/1‑150 (MQ=255)
         aTAGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATGCGGCGTAAAcgccttat                                                                                                                                              <  2:1477599/150‑5 (MQ=38)
          tAGTACTGATGATACTTCCTTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCACTGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATc                                                                                                                                             <  2:609242/150‑1 (MQ=255)
          tAGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATGCGGCGTAAAcgccttatc                                                                                                                                             <  2:1464104/150‑6 (MQ=38)
           aGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                            <  2:497653/150‑1 (MQ=38)
             tACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCTTCCGGTACTCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                                                                          <  2:1124081/150‑1 (MQ=38)
              aCTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGACAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                                                                         >  1:772926/1‑150 (MQ=255)
              aCTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTACCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                                                                         <  1:985076/150‑1 (MQ=255)
                 gatgatAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAtt                                                                                                                                      >  2:983557/1‑150 (MQ=255)
                 gatgatAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAtt                                                                                                                                      <  2:880163/150‑1 (MQ=255)
                 gatgatAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAtt                                                                                                                                      >  1:778160/1‑150 (MQ=255)
                 gatgagAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGTTAATCAATGCCTGATGCGACGCTGTCGCGTATTATCAGACCTGGAGCAAATGCAAAATGTAGGCCGGCTAAGGCGTATACGACGCATCCGGCAtt                                                                                                                                      >  1:671564/1‑150 (MQ=38)
                  atgatAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGcctac                                                                                                                                     <  2:1564153/150‑5 (MQ=38)
                   tgatAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTg                                                                                                                                    <  1:1313326/150‑1 (MQ=255)
                    gatAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGCTGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGc                                                                                                                                   <  1:1729556/150‑1 (MQ=255)
                    gatAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGc                                                                                                                                   <  2:145133/150‑1 (MQ=255)
                      tAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCtt                                                                                                                                 <  2:543126/150‑1 (MQ=255)
                         ttCATTATTTCGCCTGTCTTATCAGTCCGTGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGc                                                                                                                              <  1:141328/150‑1 (MQ=255)
                         ttCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGc                                                                                                                              >  2:1070573/1‑150 (MQ=255)
                          tCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCt                                                                                                                             >  1:1189633/1‑150 (MQ=255)
                               attTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGCCCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCagag                                                                                                                        >  2:525138/1‑150 (MQ=255)
                               attTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCagag                                                                                                                        <  1:1460794/150‑1 (MQ=255)
                               attTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCagag                                                                                                                        <  2:1630647/150‑1 (MQ=255)
                               attTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCagag                                                                                                                        >  2:844673/1‑150 (MQ=255)
                               attTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCagag                                                                                                                        >  1:153222/1‑150 (MQ=255)
                                   cGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAAAGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGACGGATAAGTCGTTTACGCCGCATACGGCATTTGCTTAGCTCAGAGTGaa                                                                                                                    >  2:282060/1‑150 (MQ=255)
                                      cTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGAt                                                                                                                 >  1:1316667/1‑150 (MQ=255)
                                      cTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTTGCTTAGCTCAGAGTGAAGAt                                                                                                                 >  2:339233/1‑150 (MQ=255)
                                       tGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGAtt                                                                                                                >  2:1314565/1‑150 (MQ=255)
                                       tGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCAGGCGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCTGCATTTGCTTAGCTCAGAGTGAAGAtt                                                                                                                >  2:1113630/1‑150 (MQ=255)
                                          tttctCAGTCCGAGTGTAGGCTGCATACGGTAATCAATGCCTTATGCGACGCTGTCGAGTCTTATCAGACATGGAGCAATTTCCAAATGTAGGCCTGATAAGGCGTTTACGCCGCATCAGGCATTTGCTTAGCTCAGAGTGAAGATTAAc                                                                                                             <  1:1113699/146‑1 (MQ=38)
                                           ttATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAAcc                                                                                                            >  2:1730996/1‑150 (MQ=255)
                                                 gTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTg                                                                                                      >  1:488955/1‑150 (MQ=255)
                                                 gTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTg                                                                                                      >  1:782201/1‑150 (MQ=255)
                                                   ccGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGcc                                                                                                    <  2:837709/150‑1 (MQ=255)
                                                        tgtAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTc                                                                                               >  2:1599077/1‑150 (MQ=255)
                                                        tgtAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTc                                                                                               >  2:633514/1‑150 (MQ=255)
                                                         gtAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTcc                                                                                              >  1:435050/1‑150 (MQ=255)
                                                          tatgCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTccc                                                                                             <  2:124558/147‑1 (MQ=255)
                                                           aGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGCGAAGATTAACCTATCTGCCCATTcccc                                                                                            >  1:1727341/1‑150 (MQ=255)
                                                              cTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCCCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTACCCaaa                                                                                         >  1:405909/1‑150 (MQ=255)
                                                               ggCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAAt                                                                                        <  2:1217723/149‑1 (MQ=255)
                                                                gCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATg                                                                                       >  2:551110/1‑150 (MQ=255)
                                                                gCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCAGCCGGCATTTGCTTAGCTCAGCGTGAAGATTAACCTATCTGCCCATTCCCCAAATg                                                                                       >  2:1555215/1‑150 (MQ=255)
                                                                   tCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGttt                                                                                    <  2:687855/150‑1 (MQ=255)
                                                                    ccGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATACGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTa                                                                                   >  1:156007/1‑150 (MQ=255)
                                                                     cGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTaa                                                                                  >  1:1173881/1‑150 (MQ=255)
                                                                      ggTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAAc                                                                                 <  1:636432/150‑1 (MQ=255)
                                                                        tAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAAcgc                                                                               <  2:1526336/150‑1 (MQ=255)
                                                                             aaTGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGATGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCTgcgc                                                                          <  2:1074048/150‑1 (MQ=255)
                                                                                gCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTa                                                                       <  2:26973/150‑1 (MQ=255)
                                                                                      tGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAActtct                                                                 >  2:448071/1‑150 (MQ=255)
                                                                                      tGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAActtct                                                                 >  2:762381/1‑150 (MQ=255)
                                                                                         gACGCTGTCGCGTCTTATCAGACCTGGAGCCATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCCTTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTcc                                                              <  1:919724/150‑1 (MQ=255)
                                                                                         gACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTcc                                                              <  1:633514/150‑1 (MQ=255)
                                                                                            gCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCAGTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGc                                                           <  1:97105/150‑1 (MQ=255)
                                                                                             cTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCt                                                          >  2:41169/1‑150 (MQ=255)
                                                                                              acacgcgTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCtt                                                         >  2:1006633/4‑150 (MQ=255)
                                                                                                 cgcgTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATAAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTAcc                                                      >  1:219796/1‑150 (MQ=255)
                                                                                                    gTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGc                                                   <  1:542191/150‑1 (MQ=255)
                                                                                                      cTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACAAGCGt                                                 >  1:1516749/1‑150 (MQ=255)
                                                                                                       ttATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGCGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGtt                                                <  2:1727341/150‑1 (MQ=255)
                                                                                                        tATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTCGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTa                                               <  1:1175730/150‑1 (MQ=255)
                                                                                                        tATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTa                                               >  1:1665937/1‑150 (MQ=255)
                                                                                                         aTCAGACATGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTaa                                              <  1:386099/150‑1 (MQ=255)
                                                                                                           cAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAAc                                            <  1:983557/150‑1 (MQ=255)
                                                                                                            aGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAAc                                                                     >  2:118830/1‑124 (MQ=255)
                                                                                                            aGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAAc                                                                     <  1:118830/124‑1 (MQ=255)
                                                                                                            aGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACg                                           <  1:1200924/150‑1 (MQ=255)
                                                                                                              cccTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTCAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGa                                         <  1:857072/149‑1 (MQ=255)
                                                                                                                cTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCTCTTAAATTCTTCCTGCTTACCAGCGTTAAACGGACg                                       >  2:150627/1‑150 (MQ=255)
                                                                                                                cTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACg                                       <  2:488955/150‑1 (MQ=255)
                                                                                                                 tGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGt                                      >  2:526691/1‑150 (MQ=255)
                                                                                                                 tGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGt                                      <  1:383584/150‑1 (MQ=255)
                                                                                                                  ggAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCATATTGAATATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTg                                     <  1:282060/150‑1 (MQ=255)
                                                                                                                  ggAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTg                                     <  2:1189633/150‑1 (MQ=255)
                                                                                                                  ggAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTg                                     <  2:221974/150‑1 (MQ=255)
                                                                                                                     gCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCAt                                  >  1:1184343/1‑150 (MQ=255)
                                                                                                                       aaTTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATcc                                >  1:540524/1‑150 (MQ=255)
                                                                                                                        aTTGCCAAATGTAGGCTGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCgg                              <  1:1320041/150‑1 (MQ=255)
                                                                                                                          tGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCggg                             >  1:308327/1‑150 (MQ=255)
                                                                                                                           gCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGGGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGGCGTGCATCCGGGc                            >  1:150515/1‑150 (MQ=255)
                                                                                                                            ccAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCt                           >  2:562563/1‑150 (MQ=255)
                                                                                                                            ccAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCt                           <  1:947751/150‑1 (MQ=255)
                                                                                                                                  aaatgCCGGATGCGGCGTAAACGCCTTATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTaa                     >  1:1428695/5‑150 (MQ=255)
                                                                                                                                   tAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTaaa                    <  1:1620697/150‑1 (MQ=255)
                                                                                                                                    aGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAat                   >  2:38323/1‑150 (MQ=255)
                                                                                                                                      gCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAatat                 >  1:504720/1‑150 (MQ=255)
                                                                                                                                      gCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAatat                 >  2:491656/1‑150 (MQ=255)
                                                                                                                                          gATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATCCGc             >  1:733675/1‑150 (MQ=255)
                                                                                                                                             aaGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATCCGcaca          >  2:1471003/1‑150 (MQ=255)
                                                                                                                                               ggCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATc                <  1:1327121/142‑1 (MQ=255)
                                                                                                                                               ggCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATc                >  2:1327121/1‑142 (MQ=255)
                                                                                                                                               ggCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATCCGCACAcg        >  1:931938/1‑150 (MQ=255)
                                                                                                                                                  gTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTCCCTAAATAGCCGCACACgcgg     >  2:1272672/1‑150 (MQ=255)
                                                                                                                                                   tttACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATCCGCACACgcggc    <  1:1730996/150‑1 (MQ=255)
                                                                                                                                                     tACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATCCGCACACGCGGCGa  >  1:1650627/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTCGACAGGATAGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATCCGCACACGCGGCGA  >  NC_000913/4460335‑4460633

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