Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,306,787 G→A 16.9% intergenic (+19/+34) oppF → / ← yciU oligopeptide ABC transporter ATPase/UPF0263 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,306,7870GA16.9% 330.3 / 50.8 136intergenic (+19/+34)oppF/yciUoligopeptide ABC transporter ATPase/UPF0263 family protein
Reads supporting (aligned to +/- strand):  ref base G (61/52);  new base A (9/14);  total (70/66)
Fisher's exact test for biased strand distribution p-value = 2.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

Only 100 of 140 total aligned reads displayed.

CCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAA  >  NC_000913/1306688‑1306886
                                                                                                   |                                                                                                   
ccAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAg                                                                                                     >  1:2085957/1‑100 (MQ=255)
ccAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAg                                                                                                     <  2:2408176/100‑1 (MQ=255)
 cAGAGTGCGCCAAAACACGTTCTGTTCTGGAGGGGAGATTCAGACACTACGTTTCTTGCCTGAAAGTCGATCAGCTTTAAAATCAATAAGGGCTGACAg                                                                                                     >  1:78837/1‑99 (MQ=255)
        cgcCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGc                                                                                             >  2:2214331/1‑100 (MQ=255)
         gcCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGcc                                                                                            <  1:2712930/100‑1 (MQ=255)
          ccAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGccc                                                                                           >  1:3137416/1‑100 (MQ=255)
           cAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCt                                                                                          <  1:1598619/100‑1 (MQ=255)
           cAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCt                                                                                          <  1:1077892/100‑1 (MQ=255)
           cAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCt                                                                                          <  1:1927027/100‑1 (MQ=255)
            aaaaCACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCtt                                                                                         <  2:2597779/100‑1 (MQ=255)
            aaaaCACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCtt                                                                                         <  1:118863/100‑1 (MQ=255)
            aaaaCACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCtn                                                                                         <  2:2722206/100‑2 (MQ=255)
            aaaaCACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCATTAAAAACAATAAGGGCTGACAGTTGTCAGCCCtt                                                                                         >  1:439623/1‑100 (MQ=255)
             aaaCACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCttt                                                                                        >  1:2072235/1‑100 (MQ=255)
              aaCACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCttat                                                                                       <  2:260310/100‑3 (MQ=255)
              aaCACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCttat                                                                                       <  2:166447/100‑3 (MQ=255)
                cacGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCttattg                                                                                     <  2:1243200/100‑2 (MQ=39)
                 acGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTTAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCa                                                                                    <  1:2474002/100‑1 (MQ=255)
                 acGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCa                                                                                    <  1:2562550/100‑1 (MQ=255)
                   tctcTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACg                                                                                  <  1:1035585/97‑1 (MQ=255)
                   gTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACg                                                                                  >  2:286487/1‑100 (MQ=255)
                    tCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGc                                                                                 >  2:645551/1‑100 (MQ=255)
                      cTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACgaga                                                                               >  2:1035592/1‑97 (MQ=255)
                       tGTTCTTGAATGAAGTTTCTGACCCTCCGTTTCTTGCCTGAAATTCGATCCGCTTTAAAAACATTAACGGCTGACTACTGTCAGCCCt                                                                                          <  1:2478698/88‑1 (MQ=255)
                        gTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTaaa                                                                             <  2:740875/100‑1 (MQ=255)
                        gTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTaaa                                                                             >  1:2091/1‑100 (MQ=255)
                            tGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCg                                                                         >  2:2049801/1‑100 (MQ=255)
                            tGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCg                                                                         <  2:1392099/100‑1 (MQ=255)
                              gAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGat                                                                       <  2:2799651/100‑1 (MQ=255)
                              gAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGat                                                                       <  2:2426499/100‑1 (MQ=255)
                              gAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGat                                                                       <  2:159954/100‑1 (MQ=255)
                               aGGGGAGTTTCGGACACTCCGCTTTTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGat                                                                       <  1:1782107/99‑1 (MQ=255)
                                ggggAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGatta                                                                     <  1:1090848/100‑1 (MQ=255)
                                  ggAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGattatt                                                                   >  2:1569859/1‑100 (MQ=255)
                                   gAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGAttattt                                                                  <  2:40791/100‑1 (MQ=255)
                                    aGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGAttattta                                                                 <  1:115788/100‑1 (MQ=255)
                                     gTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGAttatttat                                                                >  2:1768238/1‑100 (MQ=255)
                                     gTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGAttatttat                                                                >  1:1924226/1‑100 (MQ=255)
                                     gTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGAttatttat                                                                >  1:2239176/1‑100 (MQ=255)
                                      tttCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATtatttatt                                                               >  2:1773858/1‑100 (MQ=255)
                                         cAGACACTTCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGATAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTc                                                              <  1:3214665/98‑1 (MQ=255)
                                          aGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTccc                                                            >  1:2928563/1‑99 (MQ=255)
                                          aGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCg                                                           >  2:302524/1‑100 (MQ=255)
                                           gACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGc                                                          <  2:1909085/100‑1 (MQ=255)
                                           gACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGc                                                          >  2:3404325/1‑100 (MQ=255)
                                            acacTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGcc                                                         <  2:561924/100‑1 (MQ=255)
                                               cTCCGTTTGTTGCCTGAAAGTCGATCCGCTTTAACAACAATAAGGGGTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGcc                                                         <  1:983533/97‑1 (MQ=255)
                                               cTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAga                                                      >  1:2775080/1‑100 (MQ=255)
                                                tCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAgat                                                     <  1:3151888/100‑1 (MQ=255)
                                                tCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAgat                                                     >  2:3583628/1‑100 (MQ=255)
                                                tCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAgat                                                     <  2:1328244/100‑1 (MQ=255)
                                                    tttCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGata                                                 <  1:3449801/100‑1 (MQ=255)
                                                    tttCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGata                                                 <  2:691784/100‑1 (MQ=255)
                                                       cTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATgg                                              >  1:1765138/1‑100 (MQ=255)
                                                        ttGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGc                                             <  2:2091/100‑1 (MQ=255)
                                                         tGCCTTAAAGTCTATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCa                                            <  1:3417018/100‑1 (MQ=255)
                                                         tGCCTGAGAGTCGATCCGCTTTACCAACAATAAGGGCTGACAGCTGCAAGACATTTTTCAAGCTCAACGCGATTATTTACTCCCGCCAGATGATATGGCa                                            >  1:2667008/1‑100 (MQ=255)
                                                           ccTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCaaa                                          >  2:3497553/1‑100 (MQ=255)
                                                           ccTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCAAGATGATATGGCaaa                                          >  1:1811580/1‑100 (MQ=255)
                                                            cTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAg                                         >  2:1528243/1‑100 (MQ=255)
                                                                 aGTCGATCCGCTTTAAAAACAATAAGGGCTGAGAGTTGTAAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTtgt                                    <  1:2679355/100‑1 (MQ=255)
                                                                 aGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTtgt                                    <  2:2239162/100‑1 (MQ=255)
                                                                  gTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGTAAAGTTtgtg                                   >  1:2201180/1‑100 (MQ=255)
                                                                  gTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTtgtg                                   <  2:676044/100‑1 (MQ=255)
                                                                   tCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGa                                  <  1:148368/100‑1 (MQ=255)
                                                                     gATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATc                                <  1:2448512/100‑1 (MQ=255)
                                                                     gATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATc                                >  2:1174592/1‑100 (MQ=255)
                                                                     gATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATc                                >  2:3258432/1‑100 (MQ=255)
                                                                     gATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATc                                >  1:2979552/1‑100 (MQ=255)
                                                                      aaaaaaaatAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTtgtg                                   <  1:2600905/88‑1 (MQ=255)
                                                                         cGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGCTGATATGGCAAAGTTTGTGATCtttt                            >  2:1372941/1‑100 (MQ=255)
                                                                         cGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCtttt                            <  1:2990313/100‑1 (MQ=255)
                                                                         cGCTTTAAAAACAATAAGGGCTGACAGTTGTCAACCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCtttt                            >  1:3304152/1‑100 (MQ=255)
                                                                          gCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCttttt                           >  1:441166/1‑100 (MQ=255)
                                                                          gCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCttttt                           >  2:1982190/1‑100 (MQ=255)
                                                                             ttAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTcgc                        >  2:1886733/1‑100 (MQ=255)
                                                                              tAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTcgcg                       >  1:1398387/1‑100 (MQ=255)
                                                                               aaaaaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCgg                      <  1:1291517/100‑1 (MQ=255)
                                                                                aaaaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGc                     <  2:3111130/100‑1 (MQ=255)
                                                                                aaaaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGc                     <  2:3112343/100‑1 (MQ=255)
                                                                                 aaaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCa                    >  2:1772216/1‑100 (MQ=255)
                                                                                 aaaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCa                    <  2:671541/100‑1 (MQ=255)
                                                                                  gtgaaaaaGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCAt                   >  2:623057/7‑100 (MQ=38)
                                                                                   aCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATa                  <  1:3258414/100‑1 (MQ=255)
                                                                                    cAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATaa                 >  1:2986361/1‑100 (MQ=255)
                                                                                      aaaaGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACa               >  2:527136/3‑100 (MQ=255)
                                                                                      aTAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACa               <  1:1593055/100‑1 (MQ=255)
                                                                                       tAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAg              <  2:154436/100‑1 (MQ=255)
                                                                                       tAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAg              <  2:1280049/100‑1 (MQ=255)
                                                                                        aaGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGa             >  2:3408558/1‑100 (MQ=255)
                                                                                         aGGGCTGACAGTTGTCAGCCCTTTTTTTTTTTTTTATCGATTATTTATTGCCGCCAGATGATTTGTCAAAGTTTGTGATCTTTTTCGCGGCATAACAGaa            >  1:548860/1‑100 (MQ=255)
                                                                                         aGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGaa            >  2:1564178/1‑100 (MQ=255)
                                                                                         aGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGaa            >  2:3670328/1‑100 (MQ=255)
                                                                                         aGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGaa            >  2:1527970/1‑100 (MQ=255)
                                                                                           ggCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATg          >  2:2430842/1‑100 (MQ=255)
                                                                                             ntGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGCTGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCg        >  2:453676/2‑100 (MQ=255)
                                                                                              tGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCgg       <  2:723571/100‑1 (MQ=255)
                                                                                               gACAGTTGTCAGCCCTTTTTTTTTTTTTTTTCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGc         >  1:710855/1‑97 (MQ=255)
                                                                                               gACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCggg      <  2:439629/100‑1 (MQ=255)
                                                                                                   gTTGTCAACCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCaaa  >  2:2727043/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAA  >  NC_000913/1306688‑1306886

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