Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,988,355 C→A 100% pseudogene (97/633 nt) ygeK ← pseudogene; putative 2‑component transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,988,3550CA92.6% 635.4 / 42.9 243pseudogene (97/633 nt)ygeKpseudogene; putative 2‑component transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (11/7);  new base A (112/113);  total (123/120)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

Only 100 of 251 total aligned reads displayed.

ATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCATATATTCATCCTGAGTGTAAATGCTTTAACGCATTCATTTATAATCATC  >  CP009273/2988210‑2988499
                                                                                                                                                 |                                                                                                                                                
atatGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCcaaaaa                                                                                                                                                <  1:2366420/148‑1 (MQ=255)
  atGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAAcc                                                                                                                                              >  2:2990651/1‑148 (MQ=255)
     aGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCa                                                                                                                                           >  1:538314/1‑148 (MQ=255)
      gCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCat                                                                                                                                          >  2:2852612/1‑148 (MQ=255)
       cATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCata                                                                                                                                         >  1:3171672/1‑148 (MQ=255)
         tctatctTATGCGATTTAACCCATTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCAtata                                                                                                                                       >  1:1626459/1‑148 (MQ=255)
          ctatctTATGCGATTTAACCCATTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATaa                                                                                                                                      >  1:1004607/1‑148 (MQ=255)
          ctatctTATGCGATTTAACCCATTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATaa                                                                                                                                      >  2:2887926/1‑148 (MQ=255)
           tatctTATGCGATTTAACCCATTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATaaa                                                                                                                                     <  1:3212956/148‑1 (MQ=255)
            atctTATGCGATTTAACCCATTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAAt                                                                                                                                    <  1:2038331/148‑1 (MQ=255)
            atctTATGCGATTTAACCCATTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAAt                                                                                                                                    >  2:1452341/1‑148 (MQ=255)
             tctTATGCGATTTAACCCATTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATc                                                                                                                                   <  2:2373730/148‑1 (MQ=255)
                tATGCGATTTAACCCATTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAgg                                                                                                                                <  1:17962/148‑1 (MQ=255)
                 aTGCGATTTAACCCATTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGa                                                                                                                               <  2:2764777/148‑1 (MQ=255)
                   gCGATTTAACCCATTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATa                                                                                                                             >  2:2852184/1‑148 (MQ=255)
                      aTTTAACCCATTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATg                                                                                                                          <  2:2803876/148‑1 (MQ=255)
                      aTTTAACCCATTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATg                                                                                                                          <  2:218512/148‑1 (MQ=255)
                      aTTTAACCCATTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATg                                                                                                                          >  2:1285698/1‑148 (MQ=255)
                       tttAACCCATTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATCAGGATAATGt                                                                                                                         >  1:1696067/1‑148 (MQ=255)
                       tttAACCCATTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGt                                                                                                                         >  2:317372/1‑148 (MQ=255)
                       tttAACCCATTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTATTCAGATCTTTATAGCCCCCAACAACCTCATATAAATCAGGATAATGt                                                                                                                         <  1:711883/148‑1 (MQ=255)
                       tttAACCCATTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATGCTCCTAAAATAAATAATTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGt                                                                                                                         >  1:3249470/1‑148 (MQ=255)
                        ttAACCCATTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTc                                                                                                                        <  2:1681513/148‑1 (MQ=255)
                           aCCCATTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCga                                                                                                                     >  2:3235100/1‑148 (MQ=255)
                              cATTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGaa                                                                                                                  >  2:2363412/1‑148 (MQ=255)
                               aTTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGaaa                                                                                                                 <  2:1240152/148‑1 (MQ=255)
                                ttttACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAAt                                                                                                                >  2:1073073/1‑148 (MQ=255)
                                ttttACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAAt                                                                                                                <  1:1509185/148‑1 (MQ=255)
                                 tttACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATc                                                                                                               <  1:919425/148‑1 (MQ=255)
                                 tttACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATc                                                                                                               <  1:1128978/148‑1 (MQ=255)
                                   tACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCa                                                                                                             <  2:2533926/148‑1 (MQ=255)
                                     cGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAAt                                                                                                           <  2:1726362/148‑1 (MQ=255)
                                        gTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTAt                                                                                                        <  1:672712/148‑1 (MQ=255)
                                         ttCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATc                                                                                                       >  2:2707137/1‑148 (MQ=255)
                                                    tcatTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCAt                                                                                            >  2:2863861/1‑148 (MQ=255)
                                                    tcatTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTCGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCAt                                                                                            >  2:534056/1‑148 (MQ=255)
                                                     catTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATa                                                                                           <  2:890815/148‑1 (MQ=255)
                                                        tCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAAc                                                                                        >  1:1273433/1‑148 (MQ=255)
                                                             cTGAAAATTATCATAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGATATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATATTGTCCGATAAATCCAATTATCCAATCAATCATAAACGgctg                                                                                   <  2:1278085/148‑1 (MQ=255)
                                                                cAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGctgctg                                                                                <  2:2587575/148‑1 (MQ=255)
                                                                  aaTTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGAt                                                                              >  2:1683262/1‑148 (MQ=255)
                                                                   aTTCTCATAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCTAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATc                                                                             <  2:3148160/148‑1 (MQ=255)
                                                                    ttatcCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCt                                                                            >  1:388200/4‑148 (MQ=255)
                                                                    ttCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCt                                                                            >  2:2325613/1‑148 (MQ=255)
                                                                        ccTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGaaa                                                                        <  2:2728611/148‑1 (MQ=255)
                                                                           aaaaTAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTa                                                                     <  2:775247/148‑1 (MQ=255)
                                                                            aaataaatATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTAc                                                                    >  2:182271/1‑148 (MQ=255)
                                                                             aataaataTTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACa                                                                   >  1:1702798/1‑148 (MQ=255)
                                                                                  atatTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTATTCAGATCTTTATAGCCCCCAACAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACaatttt                                                              >  1:2468110/1‑148 (MQ=255)
                                                                                   tatTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACaatttta                                                             >  2:1778237/1‑148 (MQ=255)
                                                                                   tatTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACaatttta                                                             >  2:3033131/1‑148 (MQ=255)
                                                                                   tatTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACaatttta                                                             >  1:1886943/1‑148 (MQ=255)
                                                                                    atTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACaattttaa                                                            >  2:1172557/1‑148 (MQ=255)
                                                                                     tttAGGCTGTATATTTGTTACACTCCGCTATAGCTTTCTACAGATCTTTATCGCCCCCAAAAACCACATATAAATCAGTATAATGTCCCATAAATCCAATTATCCCAACAATCATAAACGGCTGCTGATCTGAAACTACaattttaat                                                           <  1:1161524/148‑1 (MQ=255)
                                                                                      ttAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACaattttaatt                                                          >  1:3019867/1‑148 (MQ=255)
                                                                                       tAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACaattttaattt                                                         <  1:2128340/148‑1 (MQ=255)
                                                                                         ggCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAAttttt                                                       >  2:980667/1‑148 (MQ=255)
                                                                                               aGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCata                                                 >  2:727494/1‑148 (MQ=255)
                                                                                                gATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATCGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCatat                                                <  1:172481/148‑1 (MQ=255)
                                                                                                gATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCatat                                                <  1:259178/148‑1 (MQ=255)
                                                                                                 aTTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCatata                                               >  2:2396132/1‑148 (MQ=255)
                                                                                                  tttGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCatatat                                              >  1:1412404/1‑148 (MQ=255)
                                                                                                  tttGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCatatat                                              <  2:444407/148‑1 (MQ=255)
                                                                                                    tGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCATATATTc                                            <  2:2942011/148‑1 (MQ=255)
                                                                                                       tACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGTATAATGTCCGAGAAATCCAATTATCCCATCAATAATAAACGGCTGCTGATCTTAAACTACAATTTTAATTTTTCCCATATATTCATc                                         <  2:1574623/148‑1 (MQ=255)
                                                                                                       tACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCATATATTCATc                                         <  1:1379560/148‑1 (MQ=255)
                                                                                                       tACACTCCGCTATAGATTTCTTCAGATCTTTATAGCCCCCAAAAACCTAATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCATATATTCATc                                         <  1:3196359/148‑1 (MQ=255)
                                                                                                                 tataGATTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCATATATTCATCCTGAGTGTaa                               <  2:1505056/148‑1 (MQ=255)
                                                                                                                  atcgCTTTATTCATATATTTTTAGCCCACAAAACCCTCATATAAAACAGGATAATGTCCGATAAATCCATTTATCCCAACAATAATACACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCATATATTCATCCTGAGTGTaaa                              <  1:582761/145‑1 (MQ=255)
                                                                                                                  ataGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCATATATTCATCCTGAGTGTaaa                              <  2:1542807/148‑1 (MQ=255)
                                                                                                                   taGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCATATATTCATCCTGAGTGTAAAt                             >  2:1586639/1‑148 (MQ=255)
                                                                                                                   taGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCATATATTCATCCTGAGTGTAAAt                             <  1:1924220/148‑1 (MQ=255)
                                                                                                                   taGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCATATATTCATCCTGAGTGTAAAt                             >  2:3120231/1‑148 (MQ=255)
                                                                                                                   taGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCATATATTCATCCTGAGTGTAAAt                             >  1:2561470/1‑148 (MQ=255)
                                                                                                                   taGCTTTATTCAGATCTTTATAGCCCCCAACAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCATATATTCATCCTGAGTGTAAAt                             >  2:187394/1‑148 (MQ=255)
                                                                                                                    aGCTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCATATATTCATCCTGAGTGTAAATg                            >  2:899662/1‑148 (MQ=255)
                                                                                                                      cTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCATATATTCATCCTGAGTGTAAATGCt                          >  1:574273/1‑148 (MQ=255)
                                                                                                                      cTTTCTTCAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCATATATTCATCCTGAGTGTAAATGCt                          >  2:507131/1‑148 (MQ=255)
                                                                                                                        ttcttcAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAAGCATACACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCATATATTCATCCTGAGTGTAAATGCttt                        <  1:326778/148‑1 (MQ=255)
                                                                                                                         tcttcAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCATATATTCATCCTGAGTGTAAATGCTTTa                       <  1:2946349/148‑1 (MQ=255)
                                                                                                                           ttcAGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCATATATTCATCCTGAGTGTAAATGCTTTAAc                     >  2:268208/1‑148 (MQ=255)
                                                                                                                             cAGATCTTTAGAGCCCCCAAACACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCATATATTCATCCTGAGTGTAAATGCTTTAACGc                   <  2:1970355/148‑1 (MQ=255)
                                                                                                                             cAGATATTTATAGCCCCCAAAAACCAAATATCAATAAGGATAATGGCCGATAAATACATTTATCCCATCACTCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCATATATTCATCCTGAGTGTAAATGCTTTAACGc                   <  1:961567/148‑1 (MQ=255)
                                                                                                                              aGATCTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCATATATTCATCCTGAGTGTAAATGCTTTAACGca                  <  1:2348782/148‑1 (MQ=255)
                                                                                                                               gATCTTTATAGCCCCCAACAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCATATATTCATCCTGAGTGTAAATGCTTTAACGcat                 >  1:3161785/1‑148 (MQ=255)
                                                                                                                                 tatttATAGCCCCCAAAAACCTCATATAAATCACGATAATGTCCGAGAAATCCAATTATCCCAACAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCATATATTCATCCTGAGTGTAAATGCTTTAACGcattc               <  1:1690437/146‑1 (MQ=255)
                                                                                                                                  cTTTATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCATATATTCATCCTGAGTGTAAATGCTTTAACGcattca              >  2:1536292/1‑148 (MQ=255)
                                                                                                                                   tttATCCCCCCAAAAAACCTCTTATAATTCAGCATAATGTCAGATAAATCCAATTATCCAATCAACCTTCCACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCATATATTCATCCTGAGTGTAAATGCTTTAACGcattcat             <  1:3125802/148‑1 (MQ=255)
                                                                                                                                    ttATAGCCCCCAAAACCCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCATATATTCATCCTGAGTGTAAATGCTTTAACGcattcatt            <  2:1595834/148‑1 (MQ=255)
                                                                                                                                    ttATAGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCATATATTCATCCTGAGTGTAAATGCTTTAACGcattcatt            <  1:229394/148‑1 (MQ=255)
                                                                                                                                     tataGCCCCCAAAACCCTCATATAAATCAGTATAATGTCCGAGAAATCCAATAATCCCATCAATCATACACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCATATATTCATCCTGAGTGTAAATGCTTTAACGCATTCAttt           <  1:2889461/148‑1 (MQ=255)
                                                                                                                                     tataGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCATATATTCATCCTGAGTGTAAATGCTTTAACGCATTCAttt           <  2:5778/148‑1 (MQ=255)
                                                                                                                                     tataGCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCATATATTCATCCTGAGTGTAAATGCTTTAACGCATTCAttt           <  1:2453180/148‑1 (MQ=255)
                                                                                                                                     tataGCCCCCAAAAACCTAATATAAATCAGGATAATGTACGAGAAATCCAATTATCCCAACAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCATATATTCATCCTGAGTGTAAATGCTTTAACGCATTCAttt           <  2:2840443/148‑1 (MQ=255)
                                                                                                                                         gCCCCCAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCATATATTCATCCTGAGTGTAAATGCTTTAACGCATTCATTTATaa       <  1:2572698/148‑1 (MQ=255)
                                                                                                                                            cccAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCATATATTCATCCTGAGTGTAAATGCTTTAACGCATTCATTTATAatca    <  1:1544769/148‑1 (MQ=255)
                                                                                                                                            cccAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCATATATTCATCCTGAGTGTAAATGCTTTAACGCATTCATTTATAatca    >  2:2037312/1‑148 (MQ=255)
                                                                                                                                            cccAAAAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCATATATTCATCCTGAGTGTAAATGCTTTAACGCATTCATTTATAatca    <  1:2842181/148‑1 (MQ=255)
                                                                                                                                              caaaaaCCTCATATAAATCAGTATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCATATATTCATCCTGAGTGTAAATGCTTTAACGCATTCATTTATAatcatc  <  1:2254735/148‑1 (MQ=255)
                                                                                                                                              caaaaaCCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCATATATTCATCCTGAGTGTAAATGCTTTAACGCATTCATTTATAatcatc  >  2:3137067/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
ATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTATCCCATCAATCATAAACGGCTGCTGATCTGAAACTACAATTTTAATTTTTCCCATATATTCATCCTGAGTGTAAATGCTTTAACGCATTCATTTATAATCATC  >  CP009273/2988210‑2988499

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