Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,151,691 C→T 15.2% intergenic (+45/+20) yegJ → / ← yegK uncharacterized protein/ser/thr phosphatase‑related protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,151,6910CT15.2% 121.3 / 3.6 69intergenic (+45/+20)yegJ/yegKuncharacterized protein/ser/thr phosphatase‑related protein
Reads supporting (aligned to +/- strand):  ref base C (35/22);  new base T (4/7);  total (39/30)
Fisher's exact test for biased strand distribution p-value = 1.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAAT  >  NC_000913/2151492‑2151924
                                                                                                                                                                                                           |                                                                                                                                                                                                                                                 
aatgtgggaaGGAGTGAAAAAAAAAGAAAGAATATTTTTCAGATTGGATGTAGCCGCGGGGGAATAAAATTTAAGGCGGTTATCCCATTTATCTTTTTTTGGTCACCTATCAAAAAGAAGAAGCCGCCGAGTTAAGAACAAATTTAGTCCCCTAACCAGCCTTTTTTTTTATTGTTTAAACGGGCGGTTGAATTTTTTTCTTCATCCGCCATGAGAAACTCACTCTGCCCACAGCGCCa                                                                                                                                                                                                         <  1:1213417/230‑1 (MQ=255)
            agagaaaaaaaaaaaaaaaaaaTTTTTCAGAGTGGTTGTATACGGGGGGGAAAAAAATTAACGGCGGTTATACAATTGATCTTTTTTTGGCCCCCTACCAAAAAGAAGAAGCCGACGAGATAAGAGAAAAATTATTTCGCTAACCCGCTTTTATTTTTTTTGTTAAAACGGGCGTTTGAATTGTTGTTTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGtt                                                                                                                                                                                              <  2:149351/224‑1 (MQ=255)
                  aaaaaaaaaaaaaaGTTATTCCAGAGTTGGTGTATCCGCGGGGGAAAAAAATTTAGGGCGGTTAAACTATTGATCTTTTATTGGCCCCTTATACAAAAAAAGAAGCCGACGAGAAAAGAGCAAAATTAGTTCGATAATCAGCCTTTATATTTTTTGTTTAAACTGGCGGTTGATTTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTatcat                                                                                                                                                                                       <  1:1484154/231‑1 (MQ=255)
                   aaaaaagaaagaAGATTTTTCAGCTGGGATGTATACGCGCGGGAAAAAAATTTCCGGGGTTTATCCCTTTGATCTTTTATTGGTCCCCTATCAAAAAGAAGAGGCCAACGAGCTAAGAGAAAAATTAGTTCGCAAATCAGCCTTTATATTTATTGTTAAACCGGGCGGTTGAATTGTTTTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTatcatc                                                                                                                                                                                      <  2:1203877/235‑1 (MQ=255)
                                            tgtaTGTATACGGGGGGGAATAAAATTTAGGGCGTTTATCCCATTTATCCTTTATTGCTCCCCTAACAAAAAAAAAAAGCCGACGAGCTAAGAGAAAAATTGGTTCGCAAATAGCCCTTTATATTTATTGTTTAAACGGGCGGTTGAATTGTTTTTTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCggg                                                                                                                                                           <  2:874035/236‑1 (MQ=255)
                                                                                tataCCATTGACCTTTTATTGGTCACCTATCAAAAAAAAGAAGCCGACGAGCAAAGAGAAAAATTATTCCGCTAATCAGCCTTTATATTTATTGTTAAACCTGGGGGTTAAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATaaa                                                                                                                      <  2:698095/239‑1 (MQ=255)
                                                                                        tGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGTATTGTTGTCTTCAACCGCCATGAGAAAATCACTCTGCCCACAGCGCCAGTGCCCACGTTTTATCATCATCCCTACGCTCGTTTACCGCCGCGGTGGACAAAAAATGCTTTTACAATTCTGGGAATAAAACGAggtg                                                                                                            >  1:558064/1‑239 (MQ=255)
                                                                                        tGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTTCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCCCGTTTACCGCCGGGCTGGACAAAAAATGCTTTAACAATTTAGGTAATAAATCGAgttg                                                                                                            >  2:985986/1‑239 (MQ=255)
                                                                                         gATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCTTGAGAAAACCACTCTTCCCCCACCGCCCGTTCCCACCGTTTATTATCATTCCTACGCCCGCTTTTCCCCGGGCGGGGCCAAAACTGCCTTAACAATTCAGCTAATAAATCGAgttgt                                                                                                           >  1:779101/1‑239 (MQ=255)
                                                                                               ttATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTTTTCCAACGCCATGAGAAAATCACTCTGCCCCCAGCGCCAGGGCCAACGGTTTTTCATCCTACGTAACCCCGTTTACCCCCGGGCTGGACAAAAAATCCTTTAACAATTCCGGTAAACAATCGAGTTGTTCCCGc                                                                                                     >  1:978505/1‑239 (MQ=255)
                                                                                               ttATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTTAAACCCCCTGAGAAACTCACTCTTCCCACAGCCCCACTGCCAACGTTTTATCATCATCCCTCCGCCCGTTTTCCGCCGGGGTGGGCAAAAAAAGCTTTTACAATTTCGGTTATAAATCGCGTTGTTCCttc                                                                                                     >  2:1114988/1‑237 (MQ=255)
                                                                                                 aTTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTTTTGTCTTCCACCCCCCTGAGAAACTTACTCTTCCCACAGCGCCAGTGCCCACCTTTTATCATCATCCCTACGCTCGTTTACCCCCGGGCTGGACAAAAACTGCTTTAACAAATCAGGtaatacactcagtttttccctcgt                                                                                                     >  1:1699584/1‑220 (MQ=255)
                                                                                                 aTTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCACGAGAAACTCCCTCTTCCCACAACGCCAGTGCCAACGCTTTATCATCATCCGTAAGCTCGTTTACCGCCCGGCTGGACAAAAACTGCTTTTACAATTCAGGTAATAAATTGAGTTGTTCCTTCGt                                                                                                   >  2:490302/1‑239 (MQ=255)
                                                                                                     gTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTTATTTTTGTCTTCCACCGCCATTAGAAAATCACTCTGCCCACAGCGCCAGTGCCAACGGTTTATCATCATCCGTACCCTCGCTTACCGCCCGGCTGGACAAAAAATGCTTTTACAAATCAGGGAATAAATCGCGTTGTTCCTGCggt                                                                                                  >  1:280196/1‑234 (MQ=255)
                                                                                                     gTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATTATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGGAATAAATCGAGGTGTTCCTGCGtt                                                                                                  >  2:753391/1‑236 (MQ=255)
                                                                                                       cACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTTTTGTCTTCCACCGCCCTGAGAAACTCACTCTGCCCACAGCGCCAGTGCCCACGTTTTATCCACATCCGTACGCTCGTTTACCGCCCGGCTGGACAAAAACTGCTTTTACAATTCCGGGAATAAAACGCGTTGTTCCTGCGTttctt                                                                                              >  1:1236316/1‑237 (MQ=255)
                                                                                                       cACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAAATGTTGTCTTCAACCGCCATGAAAAACTCACTCTGCCCACAACCCCAGTGCCAACGTTTTATCATCATCCCTACGCTCGGTTACCGCCCGGCGCGCCAAAAACCGCTTTACCAATTCAGGGAATAAATAAGGTTGTCCCCTCGTtgctg                                                                                              >  1:370670/1‑238 (MQ=255)
                                                                                                        aCCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTAATTTTTTTCTTTAACCGCCATGAGAAACTCACTCTGCCCACCGCGCCAGTGCCCACGTTTTATCCTCATCCGCACGCTCGTTTTCCGCCGGGCGTGCCAAAAACTTCTTTCACACTTCAGGGAATAAAATCAGGTGTTCCTCCGTttctgc                                                                                             >  1:1227508/1‑238 (MQ=255)
                                                                                                          cTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCCAGTTGTTCCTGCGTTGCTgccg                                                                                           >  2:630231/1‑238 (MQ=255)
                                                                                                          cTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAAATCACTCTGCCCACAACGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTTGGCAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGGTTTTTCTGCGTTGCTgccg                                                                                           >  2:1627048/1‑238 (MQ=255)
                                                                                                          ataacaaaaaaaaaaaGCCAACGAGATAAGAAAAAATTTATTTCTATAACCAGCTTTTATATTTATTGTTAAACCTGGCGTTTGATTTTTTGTTTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTgccgc                                                                                          <  2:521797/233‑1 (MQ=255)
                                                                                                            accaaaaaaaaaAGGCCGACGAGCAAAGATAAAAATGAGTCCTTAACTCAGCCTTTATTTTTATTTTTAAAACTGGGGTTTAATTTGTTGTCTTCACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCa                                                                                        <  2:62107/235‑1 (MQ=255)
                                                                                                            aTCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAACGCCAGGGCCAACGGTTTATCATCATCCCGACGCTCCGTTACCGCCCGGCCGGGCAAAAAATGCTTTAACAAATCAGGGAATAAAAAGGGGTGTTCCTGCGTTGCTgccccc                                                                                         >  1:1505699/1‑238 (MQ=255)
                                                                                                             tCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCCCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTTTCATCATCCGTACGCTCGTTTTCCCCCGGGCTGGGCCAAAACTGCTTTAaaaaat                                                                                                                                  >  2:1256588/1‑197 (MQ=255)
                                                                                                             tCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAAATTTTTTTTTCAACCCCCCTTAAAAACACACTCTGCCCACAACGACAATGCCAACGGTTTATCATCATCCGTAAGCTCGTTTACCGCCCGGACGGCCAAAAACTGCTTTAACAATATCGGAAATAAACCGAGTTGTTCCTTCGTTTCTGCCGCCat                                                                                       >  1:393885/1‑238 (MQ=255)
                                                                                                             tCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGGTGTATTTTTTTCTTCTAACCCCATCGAAAACTACCCCTGCGCAACGCGCCAGGGCCAAACGTTTATCATCAACCCTACGCTCGACTATCGCCCGGCCGGGCAAAAACTTCTTTTACAATACAGGTAAAACATCTAGTTGTTTCtcgctcgcagccgcccg                                                                                       >  1:859913/1‑223 (MQ=255)
                                                                                                               aaaagaaagaacccacctaGCTAAAGAAAAATTTGTTTCTTAAACAGCCCTTTTTTTTTTTTTTTTACGTGGGGGATTAATTTTTTTCTCCCACCCGCCATGAAAAACTCCCTCTTCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCTTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTc                                                                                     <  1:490302/221‑1 (MQ=255)
                                                                                                                       gaagCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGATTTATTGTTTTCTTCAACAGCCATGAGAAACTCACTTTGCCCACAGCGCCCGAGCCAACGGTTTATTAACCTTCGGACGCTCCGTTCACCCCGGGCCGGGCCAAACCTGGTTTTAAAAATCATGGAATAAATACGGGTTTTCCTTGCTTGCAGCCGCCCTTCCAttttaa                                                                             >  1:240416/1‑235 (MQ=255)
                                                                                                                        aagCCGGCGAGCTAAGAAAAAATTTATTTCCATAACCAGCTTTTTTTTTTTTTTTTTAAACGGGGGGTTAATTTGTTGTTTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGaaaa                                                                            <  2:512202/239‑1 (MQ=255)
                                                                                                                         ccccggggagaaaagaaaaaaaTTATCTCCCTAATCCCCTGTTATATTTTTTTTTAAACCGGGCGGTTGATTTTTTGTCTAAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGaaaaa                                                                           <  2:953603/221‑1 (MQ=255)
                                                                                                                         agCCGACGAGATAAGAGAAAAATTAGTTCGTTAATCAGCTTTTATATTTATTGTTTAACCGGGGGTTTAATTTTTTGTCTCAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGaaaaa                                                                           <  2:280196/239‑1 (MQ=255)
                                                                                                                              acgaGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGGTTAATTGTTTTCTTTAACCGCCATGAGAAACTCAATCTGCCCACAGCGCCAGTGCCCCCGATTTATCATCATCCGTACGCTCGCTTACCGCCCGGGCGGGACAAAAATTCTTTAACAAATCCGGGAAAAAAACGAGTTGTTCCTCTGTTTGTGCCGCCAGTCCCTTGaaaaaaaggg                                                                      >  1:1406802/1‑234 (MQ=255)
                                                                                                                                         tgtgaaaaagaaaaaaaaTCGGCCTTTTTTTTTTTTGTTAAAACGGGGGGTTAATTTTTTGTCTCCACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAgg                                                            <  2:994599/222‑1 (MQ=255)
                                                                                                                                          cAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTAATTGTTGTCTTTAACCGCCATGAGAAAATCACTCTGCCCACAACGCCAGTGCCAAAGTTTTTTCATCATCCGTACGCCCGTTTAACGCCCGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTTCCGCCCGTTCATTGAAAAAAGGGGGAAAGAAAgggc                                                          >  1:120551/1‑237 (MQ=255)
                                                                                                                                           aaaattattcccatatacAGCCTTTATTTTTTTTTTTAAACGGGGGGGTTAATTGTTTGTCTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGcac                                                         <  2:779101/222‑1 (MQ=255)
                                                                                                                                            aaaTTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTTTCTTCAAACGCCATGAGAAAATCACTCTGCCCACAACGCCAGTGCCAACGGTTTATCATCATCCGTACGCTCGTTTACCCCCGGGCTGGGCAAAAACTGCTTTTACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCCGGCCCTTGAAAAAAGGGGGAAAGAAAAGcaca                                                        >  2:955784/1‑239 (MQ=255)
                                                                                                                                             aaTTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGGATTGTTGTCTTCCAACGCCATGAGAAACTCACTCTGCCCCCAGCGCCAGTGCCAAAGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAAATGCTTTTACAATTCTGGTAATAAATCGAGTTGTTCCTCCGTTGCTGCCGCCCGTCCATTTAAAAAAGGGGGAAAAAAAAGCACAt                                                       >  1:1586408/1‑239 (MQ=255)
                                                                                                                                                tagtcctttatccagcttttatttttatttTTAAACCTGGTGTTTAATTTTTTTTCTCAACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATgag                                                    <  2:899347/221‑1 (MQ=255)
                                                                                                                                                tAGTCCTCTAATCAGCCTTTATTTTTTTTTTATAACGGGGCGTTTGATTTTTTTTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATgag                                                    <  2:120551/239‑1 (MQ=255)
                                                                                                                                                tAGTCCGCCAATCGCCCTTTATTTTTATTGTTAAACCTGGTGGTGAATTTTTTTTCTTAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTCCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATgag                                                    <  2:370670/239‑1 (MQ=255)
                                                                                                                                                tAGTACGATAACCGGCCGTAATATTTGTTGTTAAAACTGGGGGTTAATTGTTTTCCTAAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATgag                                                    <  2:823773/239‑1 (MQ=255)
                                                                                                                                                  gTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTTAACCCCCATGAGAAACTAAATCTGCCCCAAGCGCCCGTGCCAACGTTTTATAATCATCCGTTAGCTCGTTTACCGCCCGGCTGGACAAAAACTGCTTTAACAATTCAGGTAAAAAATCGAGTTGTTCCTGCGGTGCCGCCGCCCGTCCATTGAACAACCGGGTAAAGAAAAGCACATgagga                                                  >  2:640619/1‑239 (MQ=255)
                                                                                                                                                  gTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAAACGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCCACGTTTTATCATCATCCGTACGCCCGTTTACCGCCCGGCTGGACAAAAACTGCTTTAACAATTCCGGTAATAAATCGCGGTGTTCCTGCGTTGCTGCCGCCAGTCCCTTTAAAAAAGGGGTaaagaa                                                               >  2:923485/1‑226 (MQ=255)
                                                                                                                                                   ttCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTTTTGTCTTTAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCCGTGCCCACGGTTTTTCATCATCCGTACGCCCGTTTACCCCCGGGCCGGACCAAAACTGCTTTTAAAATTTAGGGAATAAATCGAGTTGTTCCTGCGTTGCTGCCCCCAGTCCATTTAAAAACGGGGTAAAGAAAGGCACATgaggag                                                 >  2:753207/1‑239 (MQ=255)
                                                                                                                                                     ctataacgaGCTGTTATTTTTATTGTTAAACCTGGGGGTTAAATTGTTGTCTACACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGaa                                               <  2:1550690/231‑1 (MQ=255)
                                                                                                                                                     cGCTAATCAGCCTTTATTTTTGTTTTTAAAACGGGCGGTTGATTTGTTTTCTTCACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATgag                                                    <  2:693092/234‑1 (MQ=255)
                                                                                                                                                             aGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCCACGTTTTATCATCATCCGTACGCCCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCCGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCa                                        >  2:933228/1‑238 (MQ=255)
                                                                                                                                                             aGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACACCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGGCAAAAACTGCTTTTACAATTCAGGGAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCCTTGAAAAAAGGGGGAAAGAAAGGCACATGAGGAGAATTATCCa                                        >  2:244877/1‑238 (MQ=255)
                                                                                                                                                              gcttttATTTTTTTTTTTTAAACTGGGGTTTAATTTTTTGTCTCAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGc                                      <  2:1505699/235‑1 (MQ=255)
                                                                                                                                                               tccttATTTTTTTTTTTTTAACGGGGCGTTTGATTGTTTGTTTCAACCCGCCATGAGAACCTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGc                                      <  1:630231/236‑1 (MQ=255)
                                                                                                                                                               cgtttATTTTTTTTGTTTAACTGGCTGTTGTATTTTTTACTTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCa                                     <  2:285245/236‑1 (MQ=255)
                                                                                                                                                               cgtttATATTTGTTTTTAAATCTGGTGTTTGAATTGTTGTTTCAACCCGCCATGAGAATCTCACTCTGCCCACAGCGCCAGTCCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGATCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCa                                     <  2:1630111/236‑1 (MQ=255)
                                                                                                                                                               ccGTTATATTTATTGTTTAAACTGGCGGTTGTATTTTTGTTTTCCACCGCCATGAGAAAATCACTCTGCCCACAGCGCCAGTGCCCACGTTTTATCATTATCCCTACCCCCGTTTACCGCCCGGCTGGGCAAAAACTGCTTTAACAATTCAGGGAATAAAACGAGTTGTTTCTCCGTTTCTGCCGCCCGTCCATTGAAAAACGGGGGAAAGAAAGCGACATGCGGGGAATTATc                                          >  1:693092/1‑234 (MQ=255)
                                                                                                                                                               agaaagatttttttttttAAAACGGGGGGTGAATTTGTTGCTTAAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGc                                      <  2:1070370/231‑1 (MQ=255)
                                                                                                                                                                  ttATATTTATTGTTTAAACTGGCGGTTTAATTGTTGTCTTCCACCGCCCTGAGAAACTCACTCTTCCCACAATCCCAGTGCCCACGGTTTTTTATAAACCGTAAGCCCGGTTACCGCCCGGGTGGGCAAAAACCGCTTTTAAAATTCAGGTAATAAATCGAGTTTGTTCCGTGTTGCTGCCCCCCGTCC‑‑TGAAAAACGGGGGAAAAAAAAGCCCCAGGGGAGAATTATCCAgactaatt                                >  1:967434/1‑232 (MQ=255)
                                                                                                                                                                    atatTTATTGTTTAAACTGGCGGTTTAATTGTTGTCTTCAACCCCCATGAGAAACTCACTCTGCCCACAGCGCCCGTGCCAACGTTTTATCCTCATCCCTACGCTCGTTTACCGCCGGGGCGGGCAAAAAATGATTTTACAATTAAGGTAATAAATCCAGTTTTTTCTGCCTTGGTGCCGCCAGTCCATTTGAAAACGGGGGAAAGAAAAGCAAACGAGgagactgctctttaaacaa                                 >  1:994599/1‑223 (MQ=255)
                                                                                                                                                                       gtttttttttAAACCGGGGGGTTGATTTGTTGTCCAAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTaa                              <  2:240416/232‑1 (MQ=255)
                                                                                                                                                                         tattGTTTAAACTGGCGGTTGAATTTTTGTCTTTCACCGCCATGAGAAACTCACTCTGCCCACCGCGCCAGTGCCCACGGTTTTTCATCATCCGTACGCTCGTTTACCCCCCGGGTGGACAAAAACTGCTTTTACAATTCAGGTAATAAATCGAGGTGTTTCTGCGTTGCTGCCGCCCGTCCATTGTAAAACGGGGTAAAGAAAAGCACATGAGGAGAATTATTCAGCATAttttaag                            >  1:1120763/1‑234 (MQ=255)
                                                                                                                                                                         tattGTTTAAACTGGCGGTTGAATTGTTTTCTTCCACCCCCCTGAGGAAATCACTCTTCCCACAACGCCCGTGCCAACCTTTTTTTATCATCCGTACGCTCGTTTTCCCCCGGGGCGGGCAAAAAACTCTTTTACAATTACAGGAATAAAACGAGTTGTTTCCGCTGTTCTGCCGCCAGTCCCTTTTAAAACGAGGTAAAGAAAGGCACATGAGGAGAAtttttccgccaattttacg                            >  1:843084/1‑221 (MQ=255)
                                                                                                                                                                            tGTTTAAACTGGCGGTTGAATTGTTGTTTTCAACCCCCCTGAGAAAATCACTCTTCCCACAGCGCCAGTGCCAAAGTTTTATCATCATACGTAAGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGGTGTTCCTGCGTTGTTGCCGCCAGTCCATTGAAAAAAGGGGTAAAGAAAGGCACATGAGGAGAAATATCCAGCATCTTTAAAGCCAg                        >  1:980512/1‑239 (MQ=255)
                                                                                                                                                                             tattAAACTGGGCGTTTAAGTGTTGTTCTCAACCCGCCATGAGAAACTCACTCTGCCCAAAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAg                        <  2:859913/237‑1 (MQ=255)
                                                                                                                                                                              tttAAACGTGGGGGTGAATTTTTTTCTTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGt                       <  2:885555/239‑1 (MQ=255)
                                                                                                                                                                               ttAAACTGGCGGTTGAATTGTTGTCTTTAACCGCCCTGAGAAACTCACTCTGCCCACAGCGCCCGTGGCAACGTTTTTTCATCATCCGTACCCCCGCTTTCCCCCGGGCGCGACAAAAACTGCTTTAACAATTCATGGAATAAATCGAGGTGTTCCTGCGCTTCTTCCGCCAGTCCATTTAAAAAAGGGGGAAAGAAAAGCACAAGAGGAGAATTATCCAACCTATTTAAAGCCAggc                      >  1:885555/1‑236 (MQ=255)
                                                                                                                                                                               ttAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAAATCACTCTGCCCCCAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCCAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGcaca                                                        >  2:581375/1‑204 (MQ=255)
                                                                                                                                                                                    gtGGCGGTAGATTTTTTGTCCTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACTTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGa                <  2:978505/238‑1 (MQ=255)
                                                                                                                                                                                       cgggtgttttattttttttctcCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGa                <  2:1120763/221‑1 (MQ=255)
                                                                                                                                                                                        ggggTTAAATTTTTTTCTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATa              <  2:1236316/236‑1 (MQ=255)
                                                                                                                                                                                                      ttgtCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTaaag                                                                 >  2:207845/1‑175 (MQ=255)
                                                                                                                                                                                                        gtCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAAt  >  1:474019/1‑236 (MQ=255)
                                                                                                                                                                                                        gtCTTCAACCGCAATGAGAAACTCACTCTAACCACAGCGCCAGTGCCATCGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGCCACAAACTGCTTTAACATTTCCGGAAATAAATCGAATTGTTCCTGAGTTGAAGCCGCCAGTCCATTGAAAAACGGGGTAACGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCTCTGGATAGCATCCGTAAAt  >  1:1538222/1‑236 (MQ=255)
                                                                                                                                                                                                           |                                                                                                                                                                                                                                                 
AATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTAAAT  >  NC_000913/2151492‑2151924

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