Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,417,755 C→A 6.2% pseudogene (245/402 nt) ytfA → pseudogene, related to transcriptional regulators

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,417,7550CA 6.2% 1194.6 / 31.6 306pseudogene (245/402 nt)ytfApseudogene, related to transcriptional regulators
Reads supporting (aligned to +/- strand):  ref base C (145/141);  new base A (15/5);  total (160/146)
Fisher's exact test for biased strand distribution p-value = 3.90e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 326 total aligned reads displayed.

TAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTACCGCGC  >  CP009273/4417606‑4417904
                                                                                                                                                     |                                                                                                                                                     
tAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATtgctgc                                                                                                                                                       <  2:1043046/150‑1 (MQ=255)
 aaaGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTgcg                                                                                                                                                      >  2:2685545/1‑150 (MQ=255)
   aGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCt                                                                                                                                                    <  2:2609918/150‑1 (MQ=255)
    ggTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCtt                                                                                                                                                   <  2:4595311/150‑1 (MQ=255)
       gtgCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGc                                                                                                                                                >  1:4545187/1‑150 (MQ=255)
         gCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCgg                                                                                                                                              >  1:3988395/1‑150 (MQ=255)
          cATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGa                                                                                                                                             >  1:1523427/1‑150 (MQ=255)
          cATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGa                                                                                                                                             >  2:4378704/1‑150 (MQ=255)
            ttATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAGTATTAATATTGCTGCGCTTGGCGGata                                                                                                                                           >  2:3370441/1‑150 (MQ=255)
            ttATGCAAAATCACAGCGTCTCGTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATAGTGCTGCGCGTGGCGGata                                                                                                                                           >  2:2160818/1‑150 (MQ=255)
               tGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGAtatta                                                                                                                                        <  1:2101098/150‑1 (MQ=255)
               tGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGAtatta                                                                                                                                        <  1:3474028/150‑1 (MQ=255)
                  aaaaTCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTAttt                                                                                                                                     >  2:817361/1‑150 (MQ=255)
                  aaaaTCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTAttt                                                                                                                                     >  1:838087/1‑150 (MQ=255)
                     aTCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGAt                                                                                                                                  >  1:136136/1‑150 (MQ=255)
                       cacaGCGTCTCATAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGAGCTTGGCGGATATTATTTGATCa                                                                                                                                <  2:1332443/150‑1 (MQ=255)
                          aGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAAcc                                                                                                                             <  1:4404225/150‑1 (MQ=255)
                          aGCGTCTCTTAGAAATTAACAATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAAcc                                                                                                                             >  1:1008566/1‑150 (MQ=255)
                           gCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCa                                                                                                                            <  1:655680/150‑1 (MQ=255)
                            cGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAg                                                                                                                           >  2:1203865/1‑150 (MQ=255)
                             gTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGc                                                                                                                          <  2:3633684/150‑1 (MQ=255)
                             gTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGc                                                                                                                          <  1:3220755/150‑1 (MQ=255)
                              tctcTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCa                                                                                                                         >  2:4108723/1‑150 (MQ=255)
                                tctTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATa                                                                                                                       <  2:3926662/150‑1 (MQ=255)
                                tctTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATa                                                                                                                       <  1:3364551/150‑1 (MQ=255)
                                 ctTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATAc                                                                                                                      >  1:1181740/1‑150 (MQ=255)
                                   tAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGc                                                                                                                    >  1:1914428/1‑150 (MQ=255)
                                     gAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCtt                                                                                                                  <  2:1857006/150‑1 (MQ=255)
                                       aaTTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTgg                                                                                                                <  2:4224394/150‑1 (MQ=255)
                                       aaTTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTgg                                                                                                                >  1:3024101/1‑150 (MQ=255)
                                        aTTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGc                                                                                                               >  1:975763/1‑150 (MQ=255)
                                          tAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCt                                                                                                             <  2:745832/150‑1 (MQ=255)
                                           aaCCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTg                                                                                                            <  1:2673417/150‑1 (MQ=255)
                                           aaCCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGAGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTg                                                                                                            >  2:3826323/1‑150 (MQ=255)
                                             ccATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGt                                                                                                          >  2:4874953/1‑150 (MQ=255)
                                             ccATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGt                                                                                                          >  2:2803600/1‑150 (MQ=255)
                                              cATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGtt                                                                                                         <  2:3552684/150‑1 (MQ=255)
                                                  cgcACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATc                                                                                                     >  1:4016617/1‑150 (MQ=255)
                                                    cACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATcac                                                                                                   >  2:4482663/1‑150 (MQ=255)
                                                     aCTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATcaca                                                                                                  >  1:3917340/1‑150 (MQ=255)
                                                      cTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACAt                                                                                                 <  2:1917857/150‑1 (MQ=255)
                                                       ttGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATc                                                                                                >  2:3305772/1‑150 (MQ=255)
                                                        tGCAGTTAATGGAATCCTGATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCa                                                                                               <  1:414225/150‑1 (MQ=255)
                                                         gCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAg                                                                                              >  1:216859/1‑150 (MQ=255)
                                                              tAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATgg                                                                                         >  2:4282556/1‑150 (MQ=255)
                                                              tAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATgg                                                                                         <  1:2028860/150‑1 (MQ=255)
                                                                   gAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCa                                                                                    >  2:4721481/1‑150 (MQ=255)
                                                                    aaTCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCAc                                                                                   <  1:4891566/150‑1 (MQ=255)
                                                                     aTCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCAcc                                                                                  >  2:3731000/1‑150 (MQ=255)
                                                                     aTCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCAcc                                                                                  >  1:234928/1‑150 (MQ=255)
                                                                        cTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCa                                                                               <  1:2942691/150‑1 (MQ=255)
                                                                         ttATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAg                                                                              >  1:3070832/1‑150 (MQ=255)
                                                                          tATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAgg                                                                             >  2:4658966/1‑150 (MQ=255)
                                                                           aTTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACACCAGGc                                                                            >  2:1540603/1‑150 (MQ=255)
                                                                              ggATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACt                                                                         <  2:2174720/150‑1 (MQ=255)
                                                                              ggATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACt                                                                         >  1:1699406/1‑150 (MQ=255)
                                                                               gATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTg                                                                        <  1:2397380/150‑1 (MQ=255)
                                                                                aTGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTgg                                                                       >  1:628653/1‑150 (MQ=255)
                                                                                  gAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGaa                                                                     >  2:1939420/1‑150 (MQ=255)
                                                                                   aaGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAg                                                                    >  1:476962/1‑150 (MQ=255)
                                                                                   aaGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAg                                                                    >  2:2989265/1‑150 (MQ=255)
                                                                                    aGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGc                                                                   >  2:536950/1‑150 (MQ=255)
                                                                                      ggAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCa                                                                 <  2:3879531/150‑1 (MQ=255)
                                                                                      ggAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGAGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCa                                                                 >  2:4978194/1‑150 (MQ=255)
                                                                                        aataaaCACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGa                                                               >  2:5113066/4‑150 (MQ=255)
                                                                                          gaaACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACg                                                             >  2:1294491/1‑150 (MQ=255)
                                                                                          gaaACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACg                                                             <  1:3080934/150‑1 (MQ=255)
                                                                                              cacaATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAg                                                         >  2:1601733/1‑150 (MQ=255)
                                                                                                caATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGAGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGt                                                       >  1:3511261/1‑150 (MQ=255)
                                                                                                 aaTATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGtc                                                      >  1:1319322/1‑150 (MQ=255)
                                                                                                         tAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAaga                                              <  2:4545187/150‑1 (MQ=255)
                                                                                                          aaGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAagag                                             >  1:3573562/1‑150 (MQ=255)
                                                                                                           aGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAgaga                                            >  1:2403343/1‑150 (MQ=255)
                                                                                                           aGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAgaga                                            >  1:370514/1‑150 (MQ=255)
                                                                                                               aGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAAt                                        >  2:1229937/1‑150 (MQ=255)
                                                                                                                gATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATc                                       <  2:4059981/150‑1 (MQ=255)
                                                                                                                 atatCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATcc                                      <  1:5069456/150‑1 (MQ=255)
                                                                                                                        cccATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTg                               <  1:1741351/150‑1 (MQ=255)
                                                                                                                             tGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCttt                          <  2:2218523/150‑1 (MQ=255)
                                                                                                                               cAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGTTTTAACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTAAGCTGGCTTTTg                        >  2:2496377/1‑150 (MQ=255)
                                                                                                                               cAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTg                        >  1:3276853/1‑150 (MQ=255)
                                                                                                                                aGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTgg                       >  1:748283/1‑150 (MQ=255)
                                                                                                                                aGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTgg                       <  2:788441/150‑1 (MQ=255)
                                                                                                                                aGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTgg                       <  1:2504089/150‑1 (MQ=255)
                                                                                                                                aGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTgg                       <  1:1611362/150‑1 (MQ=255)
                                                                                                                                aGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTgg                       <  1:2215750/150‑1 (MQ=255)
                                                                                                                                 ggTGAATATTAATATTGCTGAGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGt                      <  1:2298820/150‑1 (MQ=255)
                                                                                                                                  gTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGtt                     >  1:1216883/1‑150 (MQ=255)
                                                                                                                                  gTGAATATTAATATTGCTGAGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGtt                     >  1:722846/1‑150 (MQ=255)
                                                                                                                                   tGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTg                    <  2:140808/150‑1 (MQ=255)
                                                                                                                                    gAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGa                   <  2:4862727/150‑1 (MQ=255)
                                                                                                                                      atattaatattGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGAcc                 <  2:1711336/150‑1 (MQ=255)
                                                                                                                                             tattGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATc          >  2:2011694/1‑150 (MQ=255)
                                                                                                                                             tattGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATc          <  2:953390/150‑1 (MQ=255)
                                                                                                                                              attGCTGAGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCt         >  2:1778361/1‑150 (MQ=255)
                                                                                                                                                 gctgcGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTAcc      <  2:928179/150‑1 (MQ=255)
                                                                                                                                                  ctgcGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTACcg     <  2:238121/150‑1 (MQ=255)
                                                                                                                                                   tgcGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTACcgc    >  2:1127095/1‑150 (MQ=255)
                                                                                                                                                    gcgcTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTACcgcg   >  1:5125044/1‑150 (MQ=255)
                                                                                                                                                     cgcTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTACCgcgc  >  1:5022259/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAACACAATATTTTTAAGCCAGATATCGACCCATTGCAGGTGAATATTAATATTGCTGCGCTTGGCGGATATTATTTGATCAACCAGCATACGCTTGGCCTGGTTTATCACATCAGTATGGTTTCACCCCAGGCACTGGAAGCCAGACGTAAGGTCATCAAAGAGACAATCCTTAGCTGGCTTTTGGTTGACCCTTCATCTACCGCGC  >  CP009273/4417606‑4417904

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: