Predicted mutation
evidence seq id position mutation annotation gene description
RA CP028101 1,054,222 (G)11→10 coding (138/150 nt) C7Y58_05125 → site‑specific DNA‑methyltransferase (Type III DNA modification enzyme)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0281011,054,2220G.93.8% 482.5 / 44.9 163coding (138/150 nt)C7Y58_05125site‑specific DNA‑methyltransferase (Type III DNA modification enzyme)
Reads supporting (aligned to +/- strand):  ref base G (9/1);  new base . (142/10);  total (151/12)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 173 total aligned reads displayed.

GTAGAAAATGGAAAAGTTAGATGGAACTAGTATGAATTTAGTTCAAGAAAATATTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGGGTACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTAAGACCAGCCAAAGAAGAATCAAAAAATTGGGATAATACAGAAAATATCTA  >  CP028101/1054078‑1054370
                                                                                                                                                |                                                                                                                                                      
gTAGAAAATGGAAAAGTTAGATGGAACTAGTATGAATTTAGTTCAAGAAAATATTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAgg                                                                                                                                                  >  1:1316336/1‑150 (MQ=255)
gTAGAAAATGGAAAAGTTAGATGGAACTAGTATGAATTTAGTTCAAGAAAATATTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAgg                                                                                                                                                  >  2:344676/1‑150 (MQ=255)
gTAGAAAATGGAAAAGTTAGATGGAACTAGTATGAATTTAGTTCAAGAAAATATTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAgg                                                                                                                                                  >  2:26494/1‑150 (MQ=255)
  aGAAAATGGAAAAGTTAGATGGAACTAGTATGAATTTAGTTCAAGAAAATATTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGAt                                                                                                                                                >  2:118717/1‑150 (MQ=255)
  aGAAAATGGAAAAGTTAGATGGAACTAGTATGAATTTAGTTCAAGAAAATATTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGAt                                                                                                                                                >  2:118718/1‑150 (MQ=255)
    aaaaTGGAAAAGTTAGATGGAACTAGTATGAATTTAGTTCAAGAAAATATTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGG‑‑‑TACAGGATTtag                                                                                                                                            <  1:1571677/150‑1 (MQ=255)
      aaTGGAAAAGTTAGATGGAACTAGTATGAATTTAGTTCAAGAAAATATTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGG‑‑TACAGGATTtagt                                                                                                                                           <  2:439382/150‑1 (MQ=255)
        tGGAAAAGTTAGATGGAACTAGTATGAATTTAGTTCAAGAAAATATTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTtagt‑a                                                                                                                                         >  2:1497480/1‑150 (MQ=255)
         ggAAAAGTTAGATGGAACTAGTATGAATTTAGTTCAAGAAAATATTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTtagtcg                                                                                                                                         >  2:1743631/1‑148 (MQ=255)
          gAAAAGTTAGATGGAACTAGTATGAATTTAGTTCAAGAAAATANTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGG‑‑TACAGGATTTAGT‑AGAt                                                                                                                                      >  1:1786382/1‑150 (MQ=255)
              aGTTAGATGGAACTAGTATGAATTTAGTTCAAGAAAATATTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑GGATTCt                                                                                                                                   >  1:2206032/1‑150 (MQ=255)
               gTTAGATGGAACTAGTATGAATTTAGTTCAAGAAAATATTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGGGTACAGGATTTAGT‑AGATTCt                                                                                                                                   >  2:259487/1‑150 (MQ=255)
                 tAGATGGAACTAGTATGAATTTAGTTCAAGAAAATATTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAg                                                                                                                                >  2:570528/1‑150 (MQ=255)
                  aGATGGAACTAGTATGAATTTAGTTCAAGAAAATATTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAgg                                                                                                                               >  1:603921/1‑150 (MQ=255)
                   gATGGAACTAGTATGAATTTAGTTCAAGAAAATATTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGG‑‑TACAGGATTTAGT‑AGATTCTAAGGaa                                                                                                                             >  1:1016770/1‑150 (MQ=255)
                            aGTATGAATTTAGTTCAAGAAAATATTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATActc                                                                                                                     >  2:461757/1‑150 (MQ=255)
                            aGTATGAATTTAGTTCAAGAAAATATTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATActc                                                                                                                     >  2:461758/1‑150 (MQ=255)
                            aGTATGAATTTAGTTCAAGAAAATATTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATActc                                                                                                                     >  1:3119655/1‑150 (MQ=255)
                              tATGAATTTAGTTCAAGAAAATATTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCtt                                                                                                                   >  1:2299201/1‑150 (MQ=255)
                              tATGAATTTAGTTCAAGAAAATATTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCtt                                                                                                                   >  1:528240/1‑150 (MQ=255)
                               aTGAATTTAGTTCAAGAAAATATTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCttt                                                                                                                  >  2:1019105/1‑150 (MQ=255)
                                tGAATTTAGTTTGAGAAAATATTAAAAAGTTAAAAGAAATTTGTCCAGAGATATTTATAGAAGGAAGAATAGCCTTTGATGTATTAAAGCAAATTTTAGTTAGGGGGTGGCG‑CACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTa                                                                                                                 <  1:830854/150‑1 (MQ=11)
                                 gAATTTAGTTCAAGAAAATATTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGG‑‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAAc                                                                                                               <  1:645275/150‑1 (MQ=255)
                                  aagttAGTTCAAGAAAAGATGAAAAAGTTAAAAGAAATTTTTCTAGAGGTGTTTATAGAAGGAGGAAGAGATTTTGATTTATTAAAGCCCATTTTAGT‑CCTGGGGTTGCTTACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAAc                                                                                                               <  1:1079738/147‑1 (MQ=2)
                                  aaTTTAGTTCAAGAAAATATTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGGGTACAGGATTTAGT‑AGATTCTAAGGAAAGCTACTCTTTaa                                                                                                                >  2:2437203/1‑150 (MQ=255)
                                      tatttCAAGAAAATATTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTAGAGAAGGAAAAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGA‑‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTtgga                                                                                                          <  2:33039/147‑1 (MQ=11)
                                        gTTCAAGAAAATATTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTACCTtggaa                                                                                                         >  1:1446828/1‑150 (MQ=255)
                                        gTTCAAGAAAATATTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTtggaa                                                                                                         >  1:3322240/1‑150 (MQ=255)
                                        gTTCAAGAAAATATTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGG‑‑TACAGGATTTAGT‑AAATTCTAAGGAAAGATACTCTTTAACTtggaat                                                                                                        >  1:3473717/1‑150 (MQ=255)
                                         ttCAAGAAAATATTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGG‑‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTtggaatg                                                                                                       <  2:958045/150‑1 (MQ=9)
                                          tCAAGAAAATATTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACATGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTtggaatg                                                                                                       >  1:3249314/1‑150 (MQ=255)
                                          tCAAGAAAATATTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGG‑‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTtggaatgg                                                                                                      >  1:544398/1‑150 (MQ=2)
                                           cAAGAAAATATTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTtggaatgg                                                                                                      >  1:1051676/1‑150 (MQ=37)
                                           cAAGAAAATATTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTtggaatgg                                                                                                      >  1:1752069/1‑150 (MQ=37)
                                           cAAGAAAATATTAAAAAGTTAAAAAAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAACAAATTTTAGTT‑CGGGGGGGGGGTACAGGATTTAGT‑AGATTCTAAGGAAAGATGCTCCTTAACCtggaatgg                                                                                                      >  2:3496849/1‑150 (MQ=11)
                                              gAAAATATTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGaaa                                                                                                   >  2:2760065/1‑150 (MQ=37)
                                               aaaaTATTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGACTGGaaaa                                                                                                  >  2:2151289/1‑150 (MQ=255)
                                               aaaaGATTAAAAAGTGAAAAGAAATTTTTCCAGAGATGTTGATAAAAGGAAGAGGAGACTTTGGTTCATGAAAGCAAATTTTAGT‑‑GGGGGGGGTGGTACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAAt                                                                                                 <  2:814058/150‑1 (MQ=11)
                                                aaaTATTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGGGTACAGGATTTAGT‑AGACTCTAAGGAAAGATACTCTTTAACTTGGAATGGaaaa                                                                                                  >  2:1744797/1‑150 (MQ=255)
                                                aaaTATTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAAt                                                                                                 >  1:3397748/1‑150 (MQ=25)
                                                aaaTATTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAAt                                                                                                 >  1:693882/1‑150 (MQ=25)
                                                aaaTACTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGGGTACCGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGaaaa                                                                                                  >  2:1744790/1‑150 (MQ=255)
                                                    atTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAag                                                                                             >  1:1169577/1‑150 (MQ=17)
                                                     ttAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTACGGAAAGATACTCTTTAACTTGGAATGGAAAATCAaga                                                                                            >  1:1299228/1‑150 (MQ=12)
                                                     ttAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAaga                                                                                            >  1:1299233/1‑150 (MQ=17)
                                                       aaaaaGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTCAACTTGGAATGGAAAATCAAGAG‑c                                                                                         >  2:843553/1‑150 (MQ=2)
                                                           aGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑Caaaa                                                                                     >  1:571985/1‑148 (MQ=11)
                                                            gTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGAc                                                                                    >  1:816393/1‑150 (MQ=11)
                                                              tAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACaa                                                                                  >  2:1476030/1‑150 (MQ=11)
                                                               aaaaGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACaaa                                                                                 >  2:852601/1‑150 (MQ=11)
                                                                  aGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAg                                                                              >  2:3505645/1‑150 (MQ=11)
                                                                     aaTTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGGGTACAGGACTTAGTCAG‑TTCTCAGAAAAGATACTCTTTAACTTGGAATGGAAAATCAA‑AGACAAAACAAATAGca                                                                            >  2:3039313/1‑150 (MQ=11)
                                                                            ccAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGaagt                                                                    >  1:619889/1‑150 (MQ=11)
                                                                             cAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAACTCAAGAG‑CAAGACAAATAGCACAAGaagta                                                                   >  2:3343659/1‑150 (MQ=2)
                                                                               gagaTATTTATAGAAGGAAGACTAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGAAAAATCACAAAAGaagtaag                                                                 >  2:2594757/1‑150 (MQ=11)
                                                                                 gaTATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑TGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAagtaagta                                                               >  2:1465078/1‑150 (MQ=2)
                                                                                       tataGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGaa                                                         >  1:2129237/1‑150 (MQ=11)
                                                                                       tataGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGACTCTAAGGAAAGACACTCCTTCACCTGGAATGGAAAACTAAAAA‑CAAAACAAACAACACAAGAAAGAAGGACAGGaa                                                         >  2:1790683/1‑150 (MQ=2)
                                                                                       tataGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TAAAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAACCAAGAG‑CAAAAAAAATAGCACAAGaagtaagcccaggac                                                         >  1:2786432/1‑142 (MQ=11)
                                                                                       tataGAAGGAAGAATAGAATTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑CACAGGATTCAGT‑AAATTCCAAGGAAAGACACCCTTCAACTTGGAATGGAAAACCAAGAG‑CAAGACAAACAGCACAAGAAGCAAGTACAGGaa                                                         >  1:2786437/1‑150 (MQ=11)
                                                                                             aggaagAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTaa                                                   >  2:545146/1‑150 (MQ=11)
                                                                                             aggaagAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTaa                                                   >  1:3524674/1‑150 (MQ=11)
                                                                                             aggaagAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTaa                                                   >  2:675444/1‑150 (MQ=11)
                                                                                             aggaagAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTaa                                                   >  2:887678/1‑150 (MQ=11)
                                                                                                aagaaTAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTAAGAc                                                >  2:579779/1‑150 (MQ=11)
                                                                                                       gACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTAAGACCAGCCaa                                         >  1:3465371/1‑150 (MQ=17)
                                                                                                        aCTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTAAGACCAGCCaaa                                        >  1:335889/1‑150 (MQ=18)
                                                                                                         cTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTAAGACCAGCCAaag                                       >  1:1936689/1‑150 (MQ=18)
                                                                                                           ttGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTAAGACCAGCCAaagaa                                     >  1:2901224/1‑150 (MQ=18)
                                                                                                             gATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTAAGACCAGCCAAagaaga                                   >  2:2288093/1‑150 (MQ=18)
                                                                                                             gATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTAAGACCAGCCAAagaaga                                   >  2:2288094/1‑150 (MQ=18)
                                                                                                             gATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTAAGACCAGCCAAagaaga                                   >  1:1479153/1‑150 (MQ=18)
                                                                                                              aTTTATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTAAGACCAGCCAAAgaagaa                                  >  2:2995899/1‑150 (MQ=21)
                                                                                                               tttATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTAAGACCAGCCAAAGAAGAAt                                 >  1:1653562/1‑150 (MQ=25)
                                                                                                               tttATTAAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTAAGACCAGCCAAAGAAGAAt                                 >  1:2513451/1‑150 (MQ=25)
                                                                                                                ttattaAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTAAGACCAGCCAAAGAAGAATc                                >  2:796603/1‑150 (MQ=25)
                                                                                                                  attaAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTAAGACCAGCCAAAGAAGAATCaa                              >  2:1528419/1‑150 (MQ=37)
                                                                                                                   ttaAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTAAGACCAGCCAAAGAAGAATCaaa                             >  1:142220/1‑150 (MQ=37)
                                                                                                                   ttaAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTAAGACCAGCCAAAGAAGAATCaaa                             >  1:2140866/1‑150 (MQ=37)
                                                                                                                    taAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTAAGACCAGCCAAAGAAGAATCaaaa                            >  2:2702054/1‑150 (MQ=37)
                                                                                                                    taAAGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTAAGACCAGCCAAAGAAGAATCaaaa                            >  2:1682311/1‑150 (MQ=37)
                                                                                                                       aGCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATCCTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTAAGACCAGCCAAAGAAGAATCAAAAAAt                         >  1:1581431/1‑150 (MQ=255)
                                                                                                                        gCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTAAGACCAGCCAAAGAAGAATCAAAAAAtt                        >  1:1527029/1‑150 (MQ=37)
                                                                                                                        gCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTAAGACCAGCCAAAGAAGAATCAAAAAAtt                        >  1:453062/1‑150 (MQ=37)
                                                                                                                        gCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGAAAAATAGCACAAGAAGTAAGTACAGGAACATTAAGACCAGCCAAAGAAGAATCAAAAAAtt                        >  1:2949696/1‑150 (MQ=255)
                                                                                                                        gCAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAAATTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTAAGACCAGCCAAAGAAGAATCAAAAAAtt                        >  2:3409523/1‑150 (MQ=255)
                                                                                                                         cAAATTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTAAGACCAGCCAAAGAAGAATCAAAAAATTg                       >  2:200005/1‑150 (MQ=37)
                                                                                                                          aaaTTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTAAGACCAGCCAAAGAAGAATCAAAAAATTgg                      >  2:411107/1‑150 (MQ=37)
                                                                                                                           aaTTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTAAGGCCAGCCAAAGAAGAATCAAAAAATTggg                     >  1:3030157/1‑150 (MQ=255)
                                                                                                                           aaTTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTAAGACCAGCCAAAGAAGAATCAAAAAATTggg                     >  1:3030158/1‑150 (MQ=37)
                                                                                                                           aaTTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTAAGACCAGCCAAAGAAGAATCAAAAAATTggg                     >  1:240800/1‑150 (MQ=37)
                                                                                                                           aaTTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTAAGACCAGCCAAAGAAGAATCAAAAAATTggg                     >  1:240801/1‑150 (MQ=37)
                                                                                                                            aTTTTAGTTAGGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTAAGACCAGCCAAAGAAGAATCAAAAAATTGGGa                    >  2:1831761/1‑150 (MQ=255)
                                                                                                                               ttttatttagggggggtgTACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTAAGACCAGCCAAAGAAGAATCAAAAAATTGGGata                  <  2:2503647/141‑1 (MQ=255)
                                                                                                                                 agttagGGGGGGGGG‑TACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTACCTTGGACTGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTAAGACCAGCCAAAGAAGAATCAAAAAATTGGGataata               >  1:189948/1‑150 (MQ=255)
                                                                                                                                   gttaggggggggggTACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTAAGACCAGCCAAAGAAGAATCAAAAAATTGGGATAATAc              >  2:3394019/5‑150 (MQ=255)
                                                                                                                                   gttaggggggggggTACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTAAGACCAGCCAAAGAAGAATCAAAAAATTGGGATAATAc              >  2:2925286/5‑150 (MQ=255)
                                                                                                                                   gttaggggggggggTACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTAAGACCAGCCAAAGAAGAATCAAAAAATTGGGATAATAc              >  2:2925282/5‑150 (MQ=255)
                                                                                                                                           ggggggTACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTAAGACCAGCCAAAGAAGAATCAAAAAATTGGGATAATACAGAAAata      >  1:3502656/1‑150 (MQ=39)
                                                                                                                                              cagTACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTA‑AGGAACATTAAGACCAGCCAAAGAAGAATCAAAAAATTGGGATAATACAGAAAATATCta  <  2:1482040/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                      
GTAGAAAATGGAAAAGTTAGATGGAACTAGTATGAATTTAGTTCAAGAAAATATTAAAAAGTTAAAAGAAATTTTTCCAGAGATATTTATAGAAGGAAGAATAGACTTTGATTTATTAAAGCAAATTTTAGTTAGGGGGGGGGGGTACAGGATTTAGT‑AGATTCTAAGGAAAGATACTCTTTAACTTGGAATGGAAAATCAAGAG‑CAAGACAAATAGCACAAGAAGTAAGTACAGGAACATTAAGACCAGCCAAAGAAGAATCAAAAAATTGGGATAATACAGAAAATATCTA  >  CP028101/1054078‑1054370

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