Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,567,278 G→T 88.4% I166I (ATC→ATA)  mcrC ← 5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,567,2780GT88.4% 611.1 / 77.7 248I166I (ATC→ATA) mcrC5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrC
Reads supporting (aligned to +/- strand):  ref base G (13/16);  new base T (98/121);  total (111/137)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 249 total aligned reads displayed.

ACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGC  >  CP009273/4567133‑4567422
                                                                                                                                                 |                                                                                                                                                
 cGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTtata                                                                                                                                               <  2:3110790/148‑1 (MQ=255)
  gTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATaa                                                                                                                                              <  1:1756365/148‑1 (MQ=255)
   ttAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAAc                                                                                                                                             >  1:379870/1‑148 (MQ=255)
        tCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAAtt                                                                                                                                        <  1:2377714/148‑1 (MQ=255)
         cACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAAttt                                                                                                                                       <  1:3123233/148‑1 (MQ=255)
         cACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAAttt                                                                                                                                       <  1:1607133/148‑1 (MQ=255)
           cgacgGCAAAATTCATAAAGAAACTTTTGATAGAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTtat                                                                                                                                     <  2:1207134/148‑1 (MQ=255)
            gacgGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTtata                                                                                                                                    <  2:3126627/148‑1 (MQ=255)
                gCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTAtaataa                                                                                                                                <  2:1506293/148‑1 (MQ=255)
                gCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTAtaataa                                                                                                                                <  2:2043265/148‑1 (MQ=255)
                gCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTAtaataa                                                                                                                                <  2:597814/148‑1 (MQ=255)
                  aaaaTTCATAAATAAACTTTTTATAAAGTAATGACATATATTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCATTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAAcg                                                                                                                              <  2:2720268/148‑1 (MQ=255)
                  aaaaTTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAAcg                                                                                                                              >  2:2125582/1‑148 (MQ=255)
                   aaaTTCATAAAGAAACTTTTGATAAAGTAATGCCATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAAcgc                                                                                                                             >  1:1143366/1‑148 (MQ=255)
                   aaaTTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAAcgc                                                                                                                             >  1:1660874/1‑148 (MQ=255)
                    aaTTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAAcgcg                                                                                                                            >  2:3275054/1‑148 (MQ=255)
                    aaTTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAAcgcg                                                                                                                            <  1:2625965/148‑1 (MQ=255)
                    aaTTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAAcgcg                                                                                                                            >  1:3235073/1‑148 (MQ=255)
                      ttCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTa                                                                                                                          >  1:101624/1‑148 (MQ=255)
                       tCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTAt                                                                                                                         <  2:2273411/148‑1 (MQ=255)
                         aTAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGGTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTAttt                                                                                                                       >  1:952368/1‑148 (MQ=255)
                            aaGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTAtttttt                                                                                                                    <  2:482171/148‑1 (MQ=255)
                            aaGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTAttt                                                                                                                       >  1:2421801/1‑145 (MQ=255)
                                  cTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCa                                                                                                              <  1:1232833/148‑1 (MQ=255)
                                    tttGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCAtt                                                                                                            <  1:615713/148‑1 (MQ=255)
                                          aaaGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATa                                                                                                      >  1:633342/1‑148 (MQ=255)
                                          aaaGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGAGTATTTTTTCCGCCATTCAGATa                                                                                                      >  2:1688243/1‑148 (MQ=255)
                                           aaGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAg                                                                                                     <  1:1094737/148‑1 (MQ=255)
                                            aGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGc                                                                                                    >  2:2363475/1‑148 (MQ=255)
                                             gTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCt                                                                                                   >  1:1829598/1‑148 (MQ=255)
                                               aaTGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGa                                                                                                 >  1:250383/1‑148 (MQ=255)
                                                aTGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGaa                                                                                                >  2:1144846/1‑148 (MQ=255)
                                                       ctctTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAAtgctgc                                                                                         >  2:33703/1‑148 (MQ=255)
                                                          tttttCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGa                                                                                      <  2:621422/148‑1 (MQ=255)
                                                           ttttCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAg                                                                                     >  1:811934/1‑148 (MQ=255)
                                                           ttttCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAg                                                                                     >  2:1760566/1‑148 (MQ=255)
                                                               cGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTaa                                                                                 <  1:1213562/148‑1 (MQ=255)
                                                                 tttctttcAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAAt                                                                               >  2:2405539/1‑148 (MQ=255)
                                                                  ttctttcAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATg                                                                              <  2:911101/148‑1 (MQ=255)
                                                                   tctttcAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGa                                                                             <  1:1907101/148‑1 (MQ=255)
                                                                   tctttcAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGa                                                                             <  1:1629699/148‑1 (MQ=255)
                                                                      ttcAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAaga                                                                          <  1:2286091/148‑1 (MQ=255)
                                                                       tcAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAagag                                                                         <  2:661602/148‑1 (MQ=255)
                                                                       tcAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAagag                                                                         <  2:1993112/148‑1 (MQ=255)
                                                                        cAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGt                                                                        <  2:379870/148‑1 (MQ=255)
                                                                          aaaTCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGc                                                                      <  2:1923198/148‑1 (MQ=255)
                                                                           aaTCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCt                                                                     >  2:1796063/1‑148 (MQ=255)
                                                                             tCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTATTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAAGGAAGAGTGCTaa                                                                   >  2:2614580/1‑148 (MQ=255)
                                                                              cATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAAt                                                                  <  1:2048363/148‑1 (MQ=255)
                                                                               aTAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATg                                                                 <  2:2392766/148‑1 (MQ=255)
                                                                                   aaaCGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCg                                                             <  2:2964772/148‑1 (MQ=255)
                                                                                   aaaCGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCg                                                             >  1:1754281/1‑148 (MQ=255)
                                                                                     aCGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGt                                                           >  1:92548/1‑148 (MQ=255)
                                                                                      cGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTa                                                          <  1:3271611/148‑1 (MQ=255)
                                                                                           gtgtCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAAtttt                                                     <  2:2737359/148‑1 (MQ=255)
                                                                                           gtgtCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAAtttt                                                     <  1:1190212/148‑1 (MQ=255)
                                                                                             gtCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCt                                                   <  2:1044358/148‑1 (MQ=255)
                                                                                             gtCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCt                                                   <  2:511371/148‑1 (MQ=255)
                                                                                              tCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAAAGTGCTAATGACCGGTAATTTTCta                                                  >  2:2576081/1‑148 (MQ=255)
                                                                                               ccTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCtat                                                 <  2:2615190/148‑1 (MQ=255)
                                                                                                cTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCtata                                                >  1:1404148/1‑148 (MQ=255)
                                                                                                cTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCtata                                                >  1:1004941/1‑148 (MQ=255)
                                                                                                  ttGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATaaa                                              <  1:1996061/148‑1 (MQ=255)
                                                                                                  ttGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATaaa                                              >  1:2575290/1‑148 (MQ=255)
                                                                                                      tttGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGa                                          >  1:3146096/1‑148 (MQ=255)
                                                                                                       ttGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGaa                                         <  2:689936/148‑1 (MQ=255)
                                                                                                        tGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAAc                                        >  2:848559/1‑148 (MQ=255)
                                                                                                          aCCTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGa                                      <  1:1366869/148‑1 (MQ=255)
                                                                                                            cTGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGc                                    >  2:2031082/1‑148 (MQ=255)
                                                                                                             tGGAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCt                                   >  1:2859256/1‑148 (MQ=255)
                                                                                                              gtaaTAGAATTATTGACTATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAAAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCtt                                  <  1:750701/146‑1 (MQ=255)
                                                                                                              ggAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCtt                                  >  1:1725641/1‑148 (MQ=255)
                                                                                                               gAATAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTc                                 >  1:2756283/1‑148 (MQ=255)
                                                                                                                aaTAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAATAGTGCTAATGCCCGGTAATTTTCTATAAATTGAACGAGCTTCa                                >  2:2336189/1‑148 (MQ=255)
                                                                                                                aaTAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCa                                >  2:2586839/1‑148 (MQ=255)
                                                                                                                aaTAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCa                                <  1:2090502/148‑1 (MQ=255)
                                                                                                                 aTAGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCAt                               <  2:793415/148‑1 (MQ=255)
                                                                                                                   aGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCAtct                             <  2:2508608/148‑1 (MQ=255)
                                                                                                                   aGAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCAtct                             <  2:1780041/148‑1 (MQ=255)
                                                                                                                    gAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCAtctc                            >  2:106956/1‑148 (MQ=255)
                                                                                                                    gAATTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCAtctc                            <  1:2100392/148‑1 (MQ=255)
                                                                                                                     aaTTATTGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCAGCCATACAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACTAGCTTCATctct                           >  1:1722230/1‑148 (MQ=255)
                                                                                                                      attattGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTg                          >  2:2611240/1‑148 (MQ=255)
                                                                                                                      attattGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTg                          <  1:3082412/148‑1 (MQ=255)
                                                                                                                      attattGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTg                          <  2:2926338/148‑1 (MQ=255)
                                                                                                                        tattGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGAt                        <  2:1332711/148‑1 (MQ=255)
                                                                                                                         attGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATa                       <  2:1032754/148‑1 (MQ=255)
                                                                                                                          ttGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAg                      >  2:3282398/1‑148 (MQ=255)
                                                                                                                          ttGACGATGAATATGCAGAAACTTATAACGAATTTATAATAACGCGTATTTTTTCCTCCATTCAGATAGCTCAAATGATGAGGAGTTAAATGAAGAGTGCTAATGCCCGGAAATTTTCTATAAAGTGAACGAGCTTCATATCTGagat                      >  2:444143/1‑145 (MQ=255)
                                                                                                                           tGACGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGt                     <  2:1754343/148‑1 (MQ=255)
                                                                                                                              cGATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGa                  <  2:914851/148‑1 (MQ=255)
                                                                                                                               gATGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGaa                 >  2:2781251/1‑148 (MQ=255)
                                                                                                                                 tGAATTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAAtt               >  1:74989/1‑148 (MQ=255)
                                                                                                                                  gAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAAAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAAttt              >  1:1353875/1‑148 (MQ=255)
                                                                                                                                   aaTTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAAAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTa             <  1:340643/148‑1 (MQ=255)
                                                                                                                                    aTTTGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTaa            >  1:2324550/1‑148 (MQ=255)
                                                                                                                                     tttGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAAc           >  1:2627463/1‑148 (MQ=255)
                                                                                                                                     tttGCAGACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAAc           >  2:941719/1‑148 (MQ=255)
                                                                                                                                           gACACTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCa     <  1:30509/148‑1 (MQ=255)
                                                                                                                                             cacTTATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATg   <  2:250383/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
ACGTTAATTCACGACGGCAAAATTCATAAAGAAACTTTTGATAAAGTAATGACATCTCTTTTTCGTTTCTTTCAAAATCATAGAAACGGTAGTGTCCTTTGTTTTGACCTGGAATAGAATTATTGACGATGAATTTGCAGACACTGATAACGAATTTATAATAACGCGTATTTTTTCCGCCATTCAGATAGCTGAAATGCTGCGGAGTTAAATGAAGAGTGCTAATGCCCGGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGC  >  CP009273/4567133‑4567422

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