Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,417,551 A→C 29.9% 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,5510AC29.9% 264.9 / 176.1 189pseudogene (41/402 nt)ytfApseudogene, related to transcriptional regulators
Reads supporting (aligned to +/- strand):  ref base A (63/69);  new base C (28/29);  total (91/98)
Fisher's exact test for biased strand distribution p-value = 8.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 194 total aligned reads displayed.

GGGGAGAAAATGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAAC  >  CP009273/4417404‑4417700
                                                                                                                                                   |                                                                                                                                                     
ggggAGAAAATGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGctt                                                                                                                                                     >  1:4419379/1‑150 (MQ=255)
ggggAGAAAATGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGctt                                                                                                                                                     <  2:2333812/150‑1 (MQ=255)
    agaAAATGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGAttatta                                                                                                                                                 >  2:3904976/1‑150 (MQ=255)
     gaAAATGACTTCCACGCCTTACACCTTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTAGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATtattat                                                                                                                                                <  1:1072956/150‑1 (MQ=255)
     gaAAATGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTtattat                                                                                                                                                <  2:411478/150‑1 (MQ=255)
           gACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCt                                                                                                                                          <  2:1716753/150‑1 (MQ=255)
            acgtCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTa                                                                                                                                         <  2:2984828/147‑1 (MQ=255)
             cTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTaa                                                                                                                                        <  2:581983/150‑1 (MQ=255)
                ccACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATcc                                                                                                                                     <  2:84972/150‑1 (MQ=255)
                ccACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATcc                                                                                                                                     >  2:1138132/1‑150 (MQ=255)
                 cACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCt                                                                                                                                    <  1:3557756/150‑1 (MQ=255)
                 cACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCt                                                                                                                                    >  1:3962977/1‑150 (MQ=255)
                  aCGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCtt                                                                                                                                   <  2:1484247/150‑1 (MQ=255)
                    gCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTgg                                                                                                                                 >  2:86913/1‑150 (MQ=255)
                    gCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTgg                                                                                                                                 <  2:3664220/150‑1 (MQ=255)
                       ttACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGttt                                                                                                                              >  1:3521421/1‑150 (MQ=255)
                       ttACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGttt                                                                                                                              <  1:1426423/150‑1 (MQ=255)
                        tACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGtttt                                                                                                                             <  2:5051196/150‑1 (MQ=255)
                         acacCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGttttt                                                                                                                            <  2:4942665/150‑1 (MQ=255)
                         acacCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGttttt                                                                                                                            >  1:1055613/1‑150 (MQ=255)
                         acacCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGttttt                                                                                                                            >  1:1055707/1‑150 (MQ=255)
                          cacCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTa                                                                                                                           <  1:919692/150‑1 (MQ=255)
                          cacCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTa                                                                                                                           <  1:2747998/150‑1 (MQ=255)
                           acCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTaa                                                                                                                          >  1:768918/1‑150 (MQ=255)
                           acCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTaa                                                                                                                          <  1:493830/150‑1 (MQ=255)
                              gTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTaaaaa                                                                                                                       >  2:338125/1‑150 (MQ=255)
                                 cATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATtg                                                                                                                    >  2:631319/1‑150 (MQ=255)
                                   ttCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATtgtg                                                                                                                  >  2:1030854/1‑150 (MQ=255)
                                     cTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCa                                                                                                                >  1:1303335/1‑150 (MQ=255)
                                     cTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCa                                                                                                                >  1:1850019/1‑150 (MQ=255)
                                      ttgttgTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCAt                                                                                                               >  1:2503343/1‑150 (MQ=255)
                                       tgttgtGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATa                                                                                                              <  2:4141340/150‑1 (MQ=255)
                                          tgtgCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTg                                                                                                           >  1:345817/1‑150 (MQ=255)
                                               ggTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTTCAGAGTGAAAAt                                                                                                      >  2:4208914/1‑150 (MQ=255)
                                                gTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATc                                                                                                     >  2:229179/1‑150 (MQ=255)
                                                 tATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCa                                                                                                    >  1:5061922/1‑150 (MQ=255)
                                                   tGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCaaa                                                                                                  <  2:1300905/150‑1 (MQ=255)
                                                    gTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAg                                                                                                 >  2:2520663/1‑150 (MQ=255)
                                                    gTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAg                                                                                                 >  2:5098792/1‑150 (MQ=255)
                                                    gTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAg                                                                                                 >  2:1162977/1‑150 (MQ=255)
                                                         tattCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAg                                                                                            <  1:3662587/150‑1 (MQ=255)
                                                           ttCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGgt                                                                                          <  1:1942715/150‑1 (MQ=255)
                                                             cTACGCTGTCGGTATCGGTGTATTAGTTAATCTGCATGTTAGCCAGCAACAAAAAGCCATTGAGATAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGTTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGgtgt                                                                                        <  1:522693/150‑1 (MQ=255)
                                                              tACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGtgtg                                                                                       >  1:4126709/1‑150 (MQ=255)
                                                               aCGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTANCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGc                                                                                      >  1:4788689/1‑150 (MQ=255)
                                                                cGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCa                                                                                     <  2:20776/150‑1 (MQ=255)
                                                                cGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCa                                                                                     <  1:2116551/150‑1 (MQ=255)
                                                                cGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCa                                                                                     <  1:4169403/150‑1 (MQ=255)
                                                                 gCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCAt                                                                                    <  2:3603642/150‑1 (MQ=255)
                                                                 gCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCAt                                                                                    <  1:3634611/150‑1 (MQ=255)
                                                                      cGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGc                                                                               <  2:1333022/150‑1 (MQ=255)
                                                                       ggTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCa                                                                              >  2:3510282/1‑150 (MQ=255)
                                                                         tATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCaaa                                                                            >  1:3880184/1‑150 (MQ=255)
                                                                                      gTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGtctc                                                               >  2:4750672/1‑150 (MQ=255)
                                                                                      gTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGtctc                                                               <  1:1548281/150‑1 (MQ=255)
                                                                                          aTCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAg                                                           >  1:158101/1‑150 (MQ=255)
                                                                                            cagcagGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAATGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGaa                                                         >  2:2696014/1‑150 (MQ=255)
                                                                                            cagcagGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGaa                                                         <  2:3883466/150‑1 (MQ=255)
                                                                                            cagcagGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGaa                                                         <  2:345612/150‑1 (MQ=255)
                                                                                             agcagGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGaaa                                                        >  1:850062/1‑150 (MQ=255)
                                                                                               cagGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAAtt                                                      >  1:2539563/1‑150 (MQ=255)
                                                                                               cagGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAAtt                                                      <  1:1204126/150‑1 (MQ=255)
                                                                                                agGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAATGTGTGCATTATGCAAAATCACAGCGTCTCTGAGAAATTa                                                     >  1:3290862/1‑150 (MQ=255)
                                                                                                 ggTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTaa                                                    <  1:1221404/150‑1 (MQ=255)
                                                                                                   ttAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAAcc                                                  <  1:2518945/150‑1 (MQ=255)
                                                                                                    tAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCa                                                 >  1:33041/1‑150 (MQ=255)
                                                                                                    tAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCa                                                 >  1:1745472/1‑150 (MQ=255)
                                                                                                     agccagccACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCAt                                                <  1:724738/150‑1 (MQ=255)
                                                                                                        cagccaCAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATgcg                                             >  2:1880662/1‑150 (MQ=255)
                                                                                                          gccaCAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCa                                           >  2:3740092/1‑150 (MQ=255)
                                                                                                           ccaCAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCAc                                          <  1:3726930/150‑1 (MQ=255)
                                                                                                           ccaCAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCAc                                          >  1:4291205/1‑150 (MQ=255)
                                                                                                           ccaCAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCAc                                          <  1:3750345/150‑1 (MQ=255)
                                                                                                           ccaCAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCAc                                          <  1:3750200/150‑1 (MQ=255)
                                                                                                           ccaCAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCAc                                          >  2:477599/1‑150 (MQ=255)
                                                                                                           ccaCAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCAc                                          <  2:2742931/150‑1 (MQ=255)
                                                                                                               aaaaaGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGc                                      >  2:4359152/1‑150 (MQ=255)
                                                                                                                 aaaGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTAGCAg                                    >  1:2710431/1‑150 (MQ=255)
                                                                                                                  aaGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGt                                   >  1:3540705/1‑150 (MQ=255)
                                                                                                                   aGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGtt                                  <  2:4531362/150‑1 (MQ=255)
                                                                                                                    gCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTa                                 >  2:899650/1‑150 (MQ=255)
                                                                                                                      cATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAAt                               >  1:1376814/1‑150 (MQ=255)
                                                                                                                          gagaAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGaa                           <  2:419726/150‑1 (MQ=255)
                                                                                                                          gagaAAATTATTGATATTAAATGGGATTATTATATTGCTAATCATTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGGGTGCATAATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGaa                           >  2:4010957/1‑150 (MQ=255)
                                                                                                                            gaATATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATc                         <  1:522078/150‑1 (MQ=255)
                                                                                                                             aaaaTTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCAGAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATcc                        <  2:2450485/150‑1 (MQ=255)
                                                                                                                              aaaTTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCt                       <  2:1580300/150‑1 (MQ=255)
                                                                                                                                 ttattGATGTTACATGGGATTAGTATCTGGCTTATCCTTTGTTTTTATTATTTGTGCATAGTGTTAATCTAAGTAAAGGTGTGCATTATGCAAATTCTCAGCGTCTCTTAGTATTTAACCTTTCGCACTTGAAGTTAATGGAATCCTTAt                    <  2:4918754/150‑1 (MQ=255)
                                                                                                                                  tattGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTAtt                   >  2:468461/1‑150 (MQ=255)
                                                                                                                                   attGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTg                  <  2:2393543/150‑1 (MQ=255)
                                                                                                                                       aTTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATg              <  1:305228/150‑1 (MQ=255)
                                                                                                                                         tttACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGaa            <  1:2208534/150‑1 (MQ=255)
                                                                                                                                         tttACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGaa            <  2:3201013/150‑1 (MQ=255)
                                                                                                                                          ttACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAg           <  1:4719573/150‑1 (MQ=255)
                                                                                                                                          ttACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAg           <  2:876759/150‑1 (MQ=255)
                                                                                                                                           tACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAgg          >  1:2481373/1‑150 (MQ=255)
                                                                                                                                            aCATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAggg         <  1:2339046/150‑1 (MQ=255)
                                                                                                                                            aCATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAggg         <  2:4570590/150‑1 (MQ=255)
                                                                                                                                               tGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGgaag      >  2:2607990/1‑150 (MQ=255)
                                                                                                                                                 ggATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGgaagaa    <  1:3655680/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GGGGAGAAAATGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAAC  >  CP009273/4417404‑4417700

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: