Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,306,787 G→A 5.2% 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,7870GA 5.2% 485.4 / 8.9 151intergenic (+19/+34)oppF/yciUoligopeptide ABC transporter ATPase/UPF0263 family protein
Reads supporting (aligned to +/- strand):  ref base G (72/69);  new base A (5/3);  total (77/74)
Fisher's exact test for biased strand distribution p-value = 7.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 152 total aligned reads displayed.

AACCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAGACCAATCACCACT  >  NC_000913/1306640‑1306935
                                                                                                                                                   |                                                                                                                                                    
  cccGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTg                                                                                                                                                   >  1:63959/1‑149 (MQ=255)
   ccGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGt                                                                                                                                                  >  4:831416/1‑149 (MQ=255)
    cGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTc                                                                                                                                                 >  3:250408/1‑149 (MQ=255)
        ttCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGcc                                                                                                                                             <  4:764339/149‑1 (MQ=255)
         tCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGccc                                                                                                                                            <  4:619305/149‑1 (MQ=255)
            ggTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCtta                                                                                                                                         <  4:9484/149‑2 (MQ=255)
              ttGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGATTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAACCCttttt                                                                                                                                       >  2:1311200/1‑149 (MQ=255)
                 tgtTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAATGGCTGACAGTTGTCAGCCCTTTTTCAc                                                                                                                                    >  1:1095642/1‑149 (MQ=255)
                   ttttCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGc                                                                                                                                  >  2:960270/1‑149 (MQ=255)
                   ttttCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGc                                                                                                                                  >  2:785051/1‑149 (MQ=255)
                   ttttCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGc                                                                                                                                  >  2:534680/1‑149 (MQ=255)
                   ttttCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGc                                                                                                                                  >  1:942557/1‑149 (MQ=255)
                        cGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTaaaa                                                                                                                             <  2:608620/149‑1 (MQ=255)
                          tACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGc                                                                                                                           >  1:1439828/1‑149 (MQ=255)
                             ccGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGat                                                                                                                        <  4:497413/149‑1 (MQ=255)
                             ccGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGat                                                                                                                        <  3:508131/149‑1 (MQ=255)
                             ccGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGat                                                                                                                        <  3:1372160/149‑1 (MQ=255)
                             ccGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGat                                                                                                                        >  4:1327628/1‑149 (MQ=255)
                             ccGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGat                                                                                                                        >  3:1206916/1‑149 (MQ=255)
                             ccGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGat                                                                                                                        >  4:356143/1‑149 (MQ=255)
                                    ccGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGAttatttat                                                                                                                 <  3:1032473/149‑1 (MQ=255)
                                    ccGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTATCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAATTGTCAGCCCTTTTTCACGCTAAAAGCGAttatttat                                                                                                                 >  1:516428/1‑149 (MQ=255)
                                           ccGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGcc                                                                                                          >  2:254324/1‑149 (MQ=255)
                                            cGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCa                                                                                                         >  4:1378801/1‑149 (MQ=255)
                                            cGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCa                                                                                                         <  3:99086/149‑1 (MQ=255)
                                              gTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAga                                                                                                       >  3:880790/1‑149 (MQ=255)
                                              gTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCt                                                                                                                                           >  1:885829/1‑113 (MQ=255)
                                              gTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAACTGTCAGCCCt                                                                                                                                           <  2:885829/113‑1 (MQ=255)
                                               tCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAgat                                                                                                      <  1:550820/149‑1 (MQ=255)
                                                    agTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGatat                                                                                                 >  2:1292390/1‑149 (MQ=255)
                                                      tGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATgg                                                                                               >  2:1319681/1‑149 (MQ=255)
                                                          ccAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCaaa                                                                                           <  3:1293695/149‑1 (MQ=255)
                                                               acacGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTtg                                                                                      <  1:1468792/149‑1 (MQ=255)
                                                                cacGTCCTGTTCTGGAGGGGAGTTTCTGCCTCTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGACCTTTTTCACGCTAACAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTagt                                                                                     <  2:516428/149‑1 (MQ=255)
                                                                cacGTCCTATTCTGGCGGGTAGTTGCAGACACTTCTTTTCTTGCCTCACATTCGATCCGCTTTCGCAACGATAAGGGTTTACCTTTGTCACCCCTTTTCCACGCTACAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTtgt                                                                                     <  2:1210503/149‑1 (MQ=255)
                                                                  cGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGa                                                                                   <  3:347203/149‑1 (MQ=255)
                                                                    tCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATc                                                                                 >  4:1265401/1‑149 (MQ=255)
                                                                      cTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCtt                                                                               >  1:139743/1‑149 (MQ=255)
                                                                      cTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCtt                                                                               <  4:62627/149‑1 (MQ=255)
                                                                       tGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTTGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCttt                                                                              <  2:1444708/149‑1 (MQ=255)
                                                                         ttCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCttttt                                                                            >  2:836187/1‑149 (MQ=255)
                                                                                 gggAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCAt                                                                    >  4:359228/1‑149 (MQ=255)
                                                                                 gggAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCAt                                                                    <  2:1439828/149‑1 (MQ=255)
                                                                                   gAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTtgt                                                                                     >  2:36899/1‑130 (MQ=255)
                                                                                   gAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATaa                                                                  >  2:1279285/1‑149 (MQ=255)
                                                                                   gAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATaa                                                                  >  4:1427388/1‑149 (MQ=255)
                                                                                      tttCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAg                                                               >  4:1162086/1‑149 (MQ=255)
                                                                                      tttCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAg                                                               >  4:284914/1‑149 (MQ=255)
                                                                                        tctgACTCTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGGCAGTCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTAATCTTTTTCGCGGCATAACAGaa                                                             <  2:276701/146‑1 (MQ=255)
                                                                                         cAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAAt                                                            <  1:924374/149‑1 (MQ=255)
                                                                                            acacTCCGTTTCTTGCCTGAAAGTCGATCCGCGTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCg                                                         >  1:145559/1‑149 (MQ=255)
                                                                                             cacTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCgg                                                        <  4:718679/149‑1 (MQ=255)
                                                                                             cacTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCgg                                                        <  1:960270/149‑1 (MQ=255)
                                                                                                tCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAAc                                                                 >  1:1148471/1‑137 (MQ=255)
                                                                                                tCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAAc                                                                 <  2:1148471/137‑1 (MQ=255)
                                                                                                  cGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCaaa                                                   >  4:1132653/1‑149 (MQ=255)
                                                                                                   gTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCaaaa                                                  >  3:1166151/1‑149 (MQ=255)
                                                                                                    tttCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCTATTATTTATTCCCGCCAGATGAGATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCaaaaa                                                 <  2:63959/149‑1 (MQ=255)
                                                                                                        ttGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGGCAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATc                                             >  4:845435/1‑149 (MQ=255)
                                                                                                        ttGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATc                                             >  2:877123/1‑149 (MQ=255)
                                                                                                        ttGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATc                                             >  1:99263/1‑149 (MQ=255)
                                                                                                           ccTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGtt                                          <  1:1292390/149‑1 (MQ=255)
                                                                                                           ccTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGtt                                          <  1:474135/149‑1 (MQ=255)
                                                                                                            cTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTg                                         >  2:437900/1‑149 (MQ=255)
                                                                                                             tGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGAttatttat                                                                                                                 <  1:256068/76‑1 (MQ=255)
                                                                                                              gAAAGTCGATCCGCTTTAAAAGCAATAAGGGCGGACAGTTGTCAGCGCTTTTTCACGCTAAAAGCGATTACTTATACCCGCCAGATGATATGGCAAAGTTTGTGGTCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGAt                                       <  2:293251/149‑1 (MQ=255)
                                                                                                                  gTCGATCCGCTTTAAAAACAATAAGGGCTGACATTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAGCAGAATGCGGGCAAAAACATCGTTGATTTCa                                   <  2:148040/149‑1 (MQ=255)
                                                                                                                   tCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCAc                                  <  3:795016/149‑1 (MQ=255)
                                                                                                                   tCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCAc                                  <  3:831416/149‑1 (MQ=255)
                                                                                                                     gATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCa                                >  4:1382470/1‑149 (MQ=255)
                                                                                                                     gATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCa                                >  3:878885/1‑149 (MQ=255)
                                                                                                                      aTCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGAttattt                                                                                                                   <  3:23837/65‑1 (MQ=255)
                                                                                                                      aTCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGAttattt                                                                                                                   >  4:23837/1‑65 (MQ=255)
                                                                                                                      aTCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCAt                               >  1:276110/1‑149 (MQ=255)
                                                                                                                        ccGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCt                             <  3:1265401/149‑1 (MQ=255)
                                                                                                                        ccGCTTTAAAAACAATAAGGACTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCt                             <  3:78408/149‑1 (MQ=255)
                                                                                                                         cGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCtt                            >  4:765278/1‑149 (MQ=255)
                                                                                                                         cGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCtt                            <  2:942557/149‑1 (MQ=255)
                                                                                                                          gCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTc                           <  3:359143/149‑1 (MQ=255)
                                                                                                                           cTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCa                          <  1:1371396/149‑1 (MQ=255)
                                                                                                                           cTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCa                          <  3:560130/149‑1 (MQ=255)
                                                                                                                           cTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCa                          <  1:300842/149‑1 (MQ=255)
                                                                                                                             ttAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACt                        <  1:887803/149‑1 (MQ=255)
                                                                                                                                aaaaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTc                     >  3:1027063/1‑149 (MQ=255)
                                                                                                                                aaaaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTc                     >  1:777963/1‑149 (MQ=255)
                                                                                                                                 aaaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTcc                    <  3:1162086/149‑1 (MQ=255)
                                                                                                                                 aaaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTcc                    >  2:1203578/1‑149 (MQ=255)
                                                                                                                                  aaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCa                                   >  2:164163/1‑133 (MQ=255)
                                                                                                                                  aaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCg                   >  4:500784/1‑149 (MQ=255)
                                                                                                                                      aTAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAg               >  2:808867/1‑149 (MQ=255)
                                                                                                                                       aaaGGGCTGACAACTGTCAGCCCTTATTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAGa              >  1:1305004/2‑149 (MQ=255)
                                                                                                                                        aaGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAGAc             >  1:1292983/1‑149 (MQ=255)
                                                                                                                                          gggCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAGACCa           >  3:808847/1‑149 (MQ=255)
                                                                                                                                           ggCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAGACCaa          >  3:929924/1‑149 (MQ=255)
                                                                                                                                           ggCTGACAACTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAGACCaa          >  1:171613/1‑149 (MQ=255)
                                                                                                                                             cTGACAGTAGTCAGCCCTTTTTCACGCTAAAACCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCGTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAGACCAATc        >  2:245190/1‑149 (MQ=255)
                                                                                                                                              tGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAGACCAATca       >  4:381347/1‑149 (MQ=255)
                                                                                                                                               gACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAGACCAATcac      <  1:836187/149‑1 (MQ=255)
                                                                                                                                                  aGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAGACCAATcaccac   >  4:57634/1‑149 (MQ=255)
                                                                                                                                                   gTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAGACCAATCACCACt  >  4:1038528/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AACCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAGACCAATCACCACT  >  NC_000913/1306640‑1306935

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