Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,417,551 A→C 66.2% pseudogene (41/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,5510AC66.2% 217.9 / 154.8 154pseudogene (41/402 nt)ytfApseudogene, related to transcriptional regulators
Reads supporting (aligned to +/- strand):  ref base A (25/27);  new base C (51/51);  total (76/78)
Fisher's exact test for biased strand distribution p-value = 8.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.20e-01

Only 100 of 158 total aligned reads displayed.

CCGGGGAGAAAATGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAA  >  CP009273/4417402‑4417698
                                                                                                                                                     |                                                                                                                                                   
ccGGGGAGAAAATGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGa                                                                                                                                                     >  1:4373145/1‑150 (MQ=255)
 cGGGGAGAAAATGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGat                                                                                                                                                    <  1:2976483/150‑1 (MQ=255)
      agaAAATGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGAttatta                                                                                                                                               <  1:2406854/150‑1 (MQ=255)
      agaAAATGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGAttatta                                                                                                                                               <  1:4854363/150‑1 (MQ=255)
       gaAAATGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTtattat                                                                                                                                              <  1:3703602/150‑1 (MQ=255)
       gaAAATGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATtattat                                                                                                                                              >  1:1965045/1‑150 (MQ=255)
         aaaTGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTatat                                                                                                                                            >  2:2774453/1‑150 (MQ=255)
          aaTGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATAtt                                                                                                                                           >  2:747635/1‑150 (MQ=255)
          aaTGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATAtt                                                                                                                                           >  2:1950492/1‑150 (MQ=255)
           aTGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAACAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTg                                                                                                                                          <  1:52235/150‑1 (MQ=255)
             gACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCt                                                                                                                                        <  2:1141801/150‑1 (MQ=255)
             gACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCt                                                                                                                                        >  2:4631965/1‑150 (MQ=255)
              aCTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTa                                                                                                                                       >  1:2537706/1‑150 (MQ=255)
              aCTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTa                                                                                                                                       >  1:108370/1‑150 (MQ=255)
                    aCGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCtt                                                                                                                                 >  2:1908585/1‑150 (MQ=255)
                    aCGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCtt                                                                                                                                 >  2:4362763/1‑150 (MQ=255)
                    aCGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCtt                                                                                                                                 >  2:3193503/1‑150 (MQ=255)
                    aCGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCtt                                                                                                                                 <  2:4515775/150‑1 (MQ=255)
                      gCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTgg                                                                                                                               <  1:4257247/150‑1 (MQ=255)
                      gCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTgg                                                                                                                               <  1:4100749/150‑1 (MQ=255)
                         ttACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGttt                                                                                                                            <  1:2174998/150‑1 (MQ=255)
                           acacCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGttttt                                                                                                                          >  2:2646505/1‑150 (MQ=255)
                           acacCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGttttt                                                                                                                          >  2:2157361/1‑150 (MQ=255)
                            cacCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTa                                                                                                                         <  2:3481292/150‑1 (MQ=255)
                             acCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTaa                                                                                                                        >  2:1130218/1‑150 (MQ=255)
                             acCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTaa                                                                                                                        <  2:3862811/150‑1 (MQ=255)
                              ccGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAAAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTaaa                                                                                                                       <  2:712068/150‑1 (MQ=255)
                                 tACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAAt                                                                                                                    <  2:3426485/150‑1 (MQ=255)
                                  aCATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAAtt                                                                                                                   <  2:3964460/150‑1 (MQ=255)
                                    aTTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATtgt                                                                                                                 >  1:4342481/1‑150 (MQ=255)
                                     ttCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATtgtg                                                                                                                <  1:4420376/150‑1 (MQ=255)
                                      tCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGc                                                                                                               >  1:4690986/1‑150 (MQ=255)
                                           ttgtGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTAAATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGt                                                                                                          >  1:4272831/1‑150 (MQ=255)
                                              tgCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGaa                                                                                                       <  2:2597425/150‑1 (MQ=255)
                                              tgCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGaa                                                                                                       <  2:3091588/150‑1 (MQ=255)
                                              tgCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGaa                                                                                                       <  2:3091650/150‑1 (MQ=255)
                                                    aTGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCaa                                                                                                 >  1:248718/1‑150 (MQ=255)
                                                    aTGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCaa                                                                                                 >  1:1515511/1‑150 (MQ=255)
                                                     tGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCaaa                                                                                                <  1:1549288/150‑1 (MQ=255)
                                                      gTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAg                                                                                               >  2:3725172/1‑150 (MQ=255)
                                                      gTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAg                                                                                               >  2:55434/1‑150 (MQ=255)
                                                      gTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAg                                                                                               >  2:55424/1‑150 (MQ=255)
                                                       tAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGt                                                                                              >  1:3690754/1‑150 (MQ=255)
                                                       tAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGt                                                                                              >  1:1312581/1‑150 (MQ=255)
                                                              tCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGgtg                                                                                       >  1:2993836/1‑150 (MQ=255)
                                                                tACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGtgtg                                                                                     >  1:2644115/1‑150 (MQ=255)
                                                                tACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAAAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGtgtg                                                                                     <  2:449178/150‑1 (MQ=255)
                                                                 aCGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGc                                                                                    <  1:3200138/150‑1 (MQ=255)
                                                                 aCGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGc                                                                                    <  1:2280996/150‑1 (MQ=255)
                                                                   gCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCAt                                                                                  <  1:3215193/150‑1 (MQ=255)
                                                                    cTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCAtt                                                                                 >  2:2587659/1‑150 (MQ=255)
                                                                    cTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCAtt                                                                                 <  1:4126980/150‑1 (MQ=255)
                                                                    cTGTCGGTATCGGTGTATTAGTTTATCTGCAGGTTAGACAGCCACAAAATTCCATTGAGAATATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTTGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCAtt                                                                                 <  1:801437/150‑1 (MQ=255)
                                                                      gTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTAt                                                                               >  2:788204/1‑150 (MQ=255)
                                                                       tCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATg                                                                              <  1:2617942/150‑1 (MQ=255)
                                                                       tCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGAGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATg                                                                              >  1:3926337/1‑150 (MQ=255)
                                                                         ggTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCa                                                                            >  1:3554524/1‑150 (MQ=255)
                                                                         ggTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCa                                                                            >  1:1100438/1‑150 (MQ=255)
                                                                           tATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCaaa                                                                          <  2:4042956/150‑1 (MQ=255)
                                                                            aTCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCaaaa                                                                         >  1:2002105/1‑150 (MQ=255)
                                                                            aTCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCaaaa                                                                         >  1:2004285/1‑150 (MQ=255)
                                                                                gtgtATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATcac                                                                     >  1:1941221/1‑150 (MQ=255)
                                                                                gtgtATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATcac                                                                     <  2:1246466/150‑1 (MQ=255)
                                                                                   tATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGc                                                                  >  2:2656441/1‑150 (MQ=255)
                                                                                    aTTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCg                                                                 <  1:3001232/150‑1 (MQ=255)
                                                                                      tCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGtc                                                               <  2:493312/150‑1 (MQ=255)
                                                                                       cGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGtct                                                              <  1:3471352/150‑1 (MQ=255)
                                                                                        gTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGtctc                                                             <  1:2593750/150‑1 (MQ=255)
                                                                                        gTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGtctc                                                             >  1:807889/1‑150 (MQ=255)
                                                                                         ttAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTctct                                                            <  1:295855/150‑1 (MQ=255)
                                                                                             tCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGa                                                        >  1:2321373/1‑150 (MQ=255)
                                                                                                gcagGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAAt                                                     <  1:3337559/150‑1 (MQ=255)
                                                                                                gcagGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAAt                                                     <  1:3337582/150‑1 (MQ=255)
                                                                                                  agGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTa                                                   <  1:1174107/150‑1 (MQ=255)
                                                                                                  agGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTa                                                   <  1:3800381/150‑1 (MQ=255)
                                                                                                      tAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCa                                               <  2:2716742/150‑1 (MQ=255)
                                                                                                         ccagccaCAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATgc                                            >  2:4420609/1‑150 (MQ=255)
                                                                                                           agccaCAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATgcgc                                          >  2:2658254/1‑150 (MQ=255)
                                                                                                                    aaGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGt                                 <  2:4200170/150‑1 (MQ=255)
                                                                                                                    aaGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGt                                 <  2:4342481/150‑1 (MQ=255)
                                                                                                                      gCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTa                               <  2:3383262/150‑1 (MQ=255)
                                                                                                                      gCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTa                               <  2:3383209/150‑1 (MQ=255)
                                                                                                                          ttGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATgg                           >  1:3188327/1‑150 (MQ=255)
                                                                                                                           tGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGa                          >  2:3657648/1‑150 (MQ=255)
                                                                                                                            gagaAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGaa                         >  2:3844765/1‑150 (MQ=255)
                                                                                                                              gaAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATc                       >  2:1690800/1‑150 (MQ=255)
                                                                                                                                  attattTATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTa                   <  1:4749091/150‑1 (MQ=255)
                                                                                                                                  attattGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTa                   <  2:2464178/150‑1 (MQ=255)
                                                                                                                                  attattGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTa                   <  2:2462887/150‑1 (MQ=255)
                                                                                                                                    tattGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTAtt                 <  1:1888068/150‑1 (MQ=255)
                                                                                                                                     attGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTg                >  1:3534351/1‑150 (MQ=255)
                                                                                                                                     attGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTg                <  1:2656441/150‑1 (MQ=255)
                                                                                                                                      ttGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTgg               >  2:1465101/1‑150 (MQ=255)
                                                                                                                                       tGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGa              <  1:2356924/150‑1 (MQ=255)
                                                                                                                                          ttttACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGa           <  2:2004285/150‑1 (MQ=255)
                                                                                                                                          ttttACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGa           <  2:2002105/150‑1 (MQ=255)
                                                                                                                                            ttACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAg         >  1:3920837/1‑150 (MQ=255)
                                                                                                                                               cATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGga      >  1:2624770/1‑150 (MQ=255)
                                                                                                                                                aTGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGGCTCTTAGACATTAACCATGCGCACTTGCTGTTAATGGCATCCTTAGTGGATGAAGGgaa     >  2:3918484/1‑150 (MQ=255)
                                                                                                                                                   ggATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGgaagaa  >  1:4592785/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCGGGGAGAAAATGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAA  >  CP009273/4417402‑4417698

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: