Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006584 2,251,239 G→A P28P (CCC→CCT LY180_10740 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0065842,251,2390GA96.8% 61.1 / ‑6.1 32P28P (CCC→CCTLY180_10740hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (18/12);  total (18/14)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

AGTTCGGTCAGGCTCATGGGATACAGTTCTGATGGCGGCCAGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTGGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTT  >  CP006584/2250947‑2251515
                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                     
aGTTCGGTCAGGCTCATGGGATACAGTTCTGATGGCGGCCAGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGttt                                                                                                                                                                                                                                                                                  >  1:222721/1‑298 (MQ=21)
aGTTCGGTCAGGCTCATGGGATACAGTTCTGATGGCGGCCAGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGttt                                                                                                                                                                                                                                                                                  <  1:41295/298‑1 (MQ=21)
    cGGTCAGGCTGATGGGATACAGTTCAGGTGGCTGCCAGTGAAAGATCACTGCCACCTGGCCCAGAAGGCCATCGACCGAGAGATTTTTCGGGAACGTCACTGCGCGGGGTTCGGCGACAAACAACCGACCACCTTCCCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTtcatc                                                                                                                                                                                                                                                                              <  2:136938/298‑1 (MQ=11)
      gTCAGGCTCATGGGATACAGTTCTGATGGCGGCCAGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGCCCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCCCGGCAGCACCTTCATAAGTGCATCGACGTCGGGGTTTGCGCCCGCCGCCAGACCGTCACCCCCGCGCCTCCCGGCATTAGGttctctcga                                                                                                                                                                                                                                                                            >  2:143248/1‑291 (MQ=11)
            ctcctggattACAGTTCGGATGCGGGCCGTTAGAGTAGCGCTGCCACATCCGCAATCAGGTCAGCGACCGAGAGATTTTTCGGGAACGGCACTGCGCCGAGTTCGGCGACAAAAAACCGGCCACCTGACCGGCCAGCGCCACAAGGTCAGTCAGTTCCAGTGCGGCGCCCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTGATCAGTGCATCGGCTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGAc                                                                                                                                                                                                                                                                      <  2:222721/289‑1 (MQ=11)
              cATGGGATACAGTTCTGAGGGCGGCCAGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCt                                                                                                                                                                                                                                                                    <  1:7725/298‑1 (MQ=21)
                                                 tCACTGCCGCATCCGCCATCAGGTAGTCGACAGAAAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGGCCACCTTACCGGCCAGCGCCACCAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTg                                                                                                                                                                                                                                 <  2:400520/298‑1 (MQ=255)
                                                 tCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGCCACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCCTGACCTGTTCGATAACCTGCGCCCCACGCTTCACCggcctt                                                                                                                                                                                                                                 >  2:35783/1‑294 (MQ=255)
                                                               gCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATcac                                                                                                                                                                                                                   >  1:424619/1‑298 (MQ=21)
                                                                          aTCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGt                                                                                                                                                                                                        >  1:145669/1‑298 (MQ=21)
                                                                               cggagataTTTTGCGGGCACGTCACGGCACCGAGGTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATg                                                                                                                                                                                                   <  2:204248/296‑1 (MQ=16)
                                                                                     gaTTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTtc                                                                                                                                                                                             <  1:284319/298‑1 (MQ=18)
                                                                                                  cGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCCGTTCGATACCCTGCTCACACCGCTTGCCCGGGTTGTCCAGGGTATACCCATTTTTTTTTTTTATGGTTTTCACTCTTCTGAATc                                                                                                                                                                                >  2:72123/1‑298 (MQ=11)
                                                                                                            ccGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGCGTTTGCGACCGCCGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGCTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAAcc                                                                                                                                                                     >  2:415696/1‑298 (MQ=18)
                                                                                                                     gCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCa                                                                                                                                                            >  1:217630/1‑298 (MQ=18)
                                                                                                                        aCAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGc                                                                                                                                                         >  1:24774/1‑298 (MQ=18)
                                                                                                                                  gccccccTTACCGGCCAGCGCCACAAGGGCAGGCAGTTCCAGGGCGGCGACCCCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATg                                                                                                                                               <  2:217630/293‑1 (MQ=12)
                                                                                                                                      accTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGaaca                                                                                                                                           >  1:405553/1‑296 (MQ=18)
                                                                                                                                              ggCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACCGACACCGCGCCGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCGCACCCCGCTTGACAGGATTGTCCAGGGTATTCACATTTTCTTTGTTCAGGGTTTTCTGACTTCTGAATCGGGGTCACACGTGTCGCCAGCGCGCCCGGACGAAAATCACCGGc                                                                                                                                   >  2:209468/1‑298 (MQ=255)
                                                                                                                                              ggCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGCCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTACTCAGTGCATCGACTTCGGAGTTTGCGACCTCAGCCAGACCGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTCACCTGTTCGATAACCTGCTCACCACGCTTGACCGGCTTGTCCATGGTAATCACATTTTCTTTCTCCATGGTTTTCTACCTTCTGACTCGGGGTTCACCGCTCAGCCCGGCTTACCGGATCACACTCACAgga                                                                                                                                   >  2:376123/1‑297 (MQ=11)
                                                                                                                                                     gcCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTAATCAGCCTGACCTGTTCGATAACCTGCTCACCACGCCTCACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACCTCTTAATCGGGGTTAACCGGTCCGCCCAGCTACCCGGACGAAAAGCACAGGCAGATCtt                                                                                                                            >  2:142656/1‑298 (MQ=11)
                                                                                                                                                                       tCCAGTGCGGCGACCCCTTGCTCGGTCAGCAGGCGTGGGGGGAGTGGCGGCGGCACCTTAATCAGGGCATCGACGTCGGAGGTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGTCCTGTTCGATAACCTGCTCACCTCGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAg                                                                                                          <  2:293555/298‑1 (MQ=11)
                                                                                                                                                                                  gACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGc                                                                                               <  1:94485/298‑1 (MQ=17)
                                                                                                                                                                                   aCCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGcc                                                                                              >  1:123744/1‑298 (MQ=17)
                                                                                                                                                                                           cTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCt                                                                                      >  1:237754/1‑298 (MQ=17)
                                                                                                                                                                                                                           gcaCCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGAcc                                                      >  1:123423/1‑298 (MQ=14)
                                                                                                                                                                                                                            caCCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCATCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCa                                                     <  2:123744/298‑1 (MQ=14)
                                                                                                                                                                                                                               cTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGTCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCt                                                  <  1:72123/298‑1 (MQ=14)
                                                                                                                                                                                                                                    tCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCAGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTa                                             >  1:65923/1‑298 (MQ=14)
                                                                                                                                                                                                                                      aGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTAGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTAcc                                           >  1:195318/1‑298 (MQ=14)
                                                                                                                                                                                                                                                             gacccgaGCCCGACTGACGCCGTGCAGCGTCCCGTCATTTGGTTTGAGCAGCGTGACCTGTTCGATAACCTGCTGACCGCGCTTGACCTGATGGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACTGGTCAGCCCGGCTGACCGGATGCAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAgg                    <  2:393235/292‑1 (MQ=1)
                                                                                                                                                                                                                                                                               ccggggcggCGTGCCGGTATTGGGTTTCATCGTCGTGACCTGTGCGGTAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTATTGGTTTTCTCACTTCTGGATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGAttt  <  2:368292/290‑1 (MQ=1)
                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                     
AGTTCGGTCAGGCTCATGGGATACAGTTCTGATGGCGGCCAGTGAAATATCACTGCCACATCCGCCATCAGGTCATCGACCGAGAGATTTTTCGGGAACGTCACTGCACCGAGTTCGGCGACAAAAAACCGACCACCTTACCGGCCAGCGCCACAAGGTCAGGCAGTTCCAGTGCGGCGACCTCCTGCTCGGTCAGCATCGGTGCCGTCATGCGCGGCAGCACCTTAATCAGTGCATCGACTTCGGAGTTTGCGACCGCAGCCAGACTGACACCGCGCAGCGTCCCGGCATTGGGTTTCATCAGCGTGACCTGTTCGATAACCTGCTCACCACGCTTGACCGGATTGTCCAGGGTAATCACATTTTCTTTGTTCATGGTTTTCTCACTTCTGAATCGGGGTTAACCGGTCAGCCAGGCTGACCGGATGAAAATCACAGGCCGATATTGCGGCGGTGTTGCTCCAGCCGGTCGACGCCGTTCACCTTCTCAATCATGTTGATGGTGTCGATTTCGACCAGCTCCTTACCGTCCATCGTCAGCCGGAAATAGGTGCAGACCACGGAGATTT  >  CP006584/2250947‑2251515

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: