Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,306,787 G→A 16.4% 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.4% 298.8 / 44.0 122intergenic (+19/+34)oppF/yciUoligopeptide ABC transporter ATPase/UPF0263 family protein
Reads supporting (aligned to +/- strand):  ref base G (65/37);  new base A (13/7);  total (78/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-01

Only 100 of 123 total aligned reads displayed.

CCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCA  >  NC_000913/1306688‑1306884
                                                                                                   |                                                                                                 
ccAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAg                                                                                                   <  2:2414595/100‑1 (MQ=255)
 cAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGt                                                                                                  >  1:418881/1‑100 (MQ=255)
   gagTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTg                                                                                                <  1:908523/100‑1 (MQ=255)
     gTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTc                                                                                              <  1:2607692/100‑1 (MQ=255)
     gTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTc                                                                                              <  1:2484700/100‑1 (MQ=255)
      tGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCa                                                                                             <  2:899232/100‑1 (MQ=255)
        cgcCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGc                                                                                           >  1:2321000/1‑100 (MQ=255)
        cgcCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGc                                                                                           <  1:945049/100‑1 (MQ=255)
         gcCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAgat                                                                                          >  1:943553/1‑98 (MQ=255)
         gcCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGcc                                                                                          >  2:1116677/1‑100 (MQ=255)
          ccAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGccc                                                                                         >  1:1316221/1‑100 (MQ=255)
          ccAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGccc                                                                                         >  2:63166/1‑100 (MQ=255)
          ccAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGccc                                                                                         >  1:208678/1‑100 (MQ=255)
           cAAAACACGTCCTGTTCTGGAGGGGAGTTTGAGACACTCCGTTTCATGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCt                                                                                        <  2:637785/100‑1 (MQ=255)
           cAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCt                                                                                        <  1:1270671/100‑1 (MQ=255)
           cAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCt                                                                                        <  1:2026964/100‑1 (MQ=255)
           cAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCt                                                                                        <  2:2060953/100‑1 (MQ=255)
           cAAAACACGTCCTGTTCTGGAGAGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCt                                                                                        <  2:2013299/100‑1 (MQ=255)
            aaaaCACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCtt                                                                                       >  2:566773/1‑100 (MQ=255)
                cacGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTc                                                                                   <  1:750837/100‑1 (MQ=255)
                    tCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGc                                                                               >  2:2044947/1‑100 (MQ=255)
                    tCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGc                                                                               >  2:2591454/1‑100 (MQ=255)
                        gTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTaaa                                                                           >  1:2032795/1‑100 (MQ=255)
                        gTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTCTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTaaa                                                                           <  2:1522308/100‑1 (MQ=255)
                             ggAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGa                                                                      >  1:1686878/1‑100 (MQ=255)
                              gAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGat                                                                     >  2:2704536/1‑100 (MQ=255)
                                ggggAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAACCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGatta                                                                   <  2:2452415/100‑1 (MQ=255)
                                gggcggTTTCAGACACTCCGTTTCTTCCCTGAAAGTCGATCCGCGTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGat                                                                     <  2:2524995/93‑1 (MQ=255)
                                  ggAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGattatt                                                                 <  2:83558/100‑1 (MQ=255)
                                  ggAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGattatt                                                                 <  1:712706/100‑1 (MQ=255)
                                     gTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGAttatttat                                                              <  2:1476604/100‑1 (MQ=255)
                                     gTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGAttatttat                                                              >  1:378694/1‑100 (MQ=255)
                                        tCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAATTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTcc                                                           <  1:511856/100‑1 (MQ=255)
                                            acacTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGcc                                                       <  1:1244739/100‑1 (MQ=255)
                                              acTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAg                                                     <  2:2704981/100‑1 (MQ=255)
                                               cTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAga                                                    >  2:1414424/1‑100 (MQ=255)
                                                tCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAgat                                                   >  1:1847729/1‑100 (MQ=255)
                                                  cGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAgatga                                                 >  1:2591002/1‑100 (MQ=255)
                                                     ttCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGatat                                              <  2:418880/100‑1 (MQ=255)
                                                       cTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATgg                                            >  2:2266549/1‑100 (MQ=255)
                                                       cTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATgg                                            <  1:1223854/100‑1 (MQ=255)
                                                        ttGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGc                                           <  2:715715/100‑1 (MQ=255)
                                                          gacTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTGGGCAGCCATTTTTCACGCAAAAAGCGATTACTTATTCCGGCCATATGatattgaaa                                         >  2:2362477/3‑95 (MQ=255)
                                                          gCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCaa                                         >  1:2719050/1‑100 (MQ=255)
                                                          gCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCaa                                         >  2:2361157/1‑100 (MQ=255)
                                                          gCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCaa                                         <  2:577196/100‑1 (MQ=255)
                                                           ccTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCaaa                                        >  1:601590/1‑100 (MQ=255)
                                                              gAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGtt                                     <  1:88170/100‑1 (MQ=255)
                                                              gAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGtt                                     <  2:140949/100‑1 (MQ=255)
                                                                aaGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTtg                                   >  1:1608377/1‑100 (MQ=255)
                                                                 aGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTtgt                                  <  1:111936/100‑1 (MQ=255)
                                                                 aGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTtgt                                  >  1:78387/1‑100 (MQ=255)
                                                                 aGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTtg                                   <  1:2467878/99‑1 (MQ=255)
                                                                  gTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTtgtg                                 <  1:2003383/100‑1 (MQ=255)
                                                                  gTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTtgtg                                 <  2:2234556/100‑1 (MQ=255)
                                                                  gTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTtgtg                                 >  2:2659461/1‑100 (MQ=255)
                                                                  gTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCAAAAAGCGATTATTTATTACCGCCAGATGATCTGGCACCGttt                                    >  1:227525/1‑97 (MQ=255)
                                                                   tCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGa                                >  1:582019/1‑100 (MQ=255)
                                                                     gATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATc                              >  1:2212281/1‑100 (MQ=255)
                                                                     gATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATc                              >  2:2647826/1‑100 (MQ=255)
                                                                      aTCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCt                             >  2:2557772/1‑100 (MQ=255)
                                                                      aTCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCt                             >  1:790769/1‑100 (MQ=255)
                                                                       tCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCtt                            >  1:2228877/1‑100 (MQ=255)
                                                                       tCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCtt                            >  2:443361/1‑100 (MQ=255)
                                                                         ctgtttAAAACCCGTAAGGGCTGACAGGTGTCGGCCCTTTGTCACGCTAAAAGCGATTATTTGTTCCCGCCAGATGATATGGCAAAGTTTGTGATc                              <  2:2154785/93‑1 (MQ=255)
                                                                         cGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCtttt                          >  1:2212552/1‑100 (MQ=255)
                                                                         cGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCtttt                          >  1:2107439/1‑100 (MQ=255)
                                                                           cTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTc                        >  1:331451/1‑100 (MQ=255)
                                                                            tttAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTcg                       >  2:2569036/1‑100 (MQ=255)
                                                                              tAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTcgcg                     >  2:1559786/1‑100 (MQ=255)
                                                                              tAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTcgcg                     >  1:302543/1‑100 (MQ=255)
                                                                              tAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTcgcg                     <  2:2266096/100‑1 (MQ=255)
                                                                               aaaaaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCgg                    >  1:1518343/1‑100 (MQ=255)
                                                                                aaaaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGc                   >  2:1381035/1‑100 (MQ=255)
                                                                                aaaaCAATAAGGGCTGACAGTTGTCAGCCCATTCTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCGTTTTCGCGGc                   >  2:1803125/1‑100 (MQ=255)
                                                                                 cgtgaaaaaGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCa                  >  1:1342541/8‑100 (MQ=38)
                                                                                 aaaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCa                  >  2:2047129/1‑100 (MQ=255)
                                                                                 aaaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCa                  >  1:2309323/1‑100 (MQ=255)
                                                                                   aCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATa                <  2:928163/100‑1 (MQ=255)
                                                                                    cAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATaa               <  1:2386155/100‑1 (MQ=255)
                                                                                     aaTAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCATATGATATGGCAAATTTTGTGATTTTTTTCGCGGCATAAc              <  2:2186920/100‑1 (MQ=255)
                                                                                      aaaaGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCAt                 >  2:150879/3‑96 (MQ=255)
                                                                                      aaaaGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACa             >  1:383721/3‑100 (MQ=255)
                                                                                       aaaGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAg            >  1:1360768/2‑100 (MQ=255)
                                                                                        aaGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGa           >  2:2106863/1‑100 (MQ=255)
                                                                                         aGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGaa          >  1:649746/1‑100 (MQ=255)
                                                                                         aGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGaa          >  1:1433104/1‑100 (MQ=255)
                                                                                         aGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGaa          >  1:2087722/1‑100 (MQ=255)
                                                                                         aGGGCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTCTGCGATCTTTTTctcagtata                >  2:215005/1‑86 (MQ=255)
                                                                                           ggCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATg        >  2:2227021/1‑100 (MQ=255)
                                                                                            gCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGc       >  1:1065959/1‑100 (MQ=255)
                                                                                            gCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGc       >  2:48757/1‑100 (MQ=255)
                                                                                             cTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCg      >  1:1063626/1‑100 (MQ=255)
                                                                                             cTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCg      >  2:479015/1‑100 (MQ=255)
                                                                                             cTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAAGGCg      <  2:2610877/100‑1 (MQ=255)
                                                                                               gACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCggg    <  1:2070545/100‑1 (MQ=255)
                                                                                                 cAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATTGCAAAGTTTCTGATCTTTTTCGCGGCATAACAGAATGCTGGCa  >  2:1730798/1‑100 (MQ=255)
                                                                                                 cAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCa  >  1:168233/1‑100 (MQ=255)
                                                                                                 cAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCa  >  1:2549612/1‑100 (MQ=255)
                                                                                                 cAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCa  >  2:289780/1‑100 (MQ=255)
                                                                                                   |                                                                                                 
CCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCA  >  NC_000913/1306688‑1306884

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: