Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,417,551 A→C 93.7% 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,5510AC93.7% 1087.2 / 32.7 286pseudogene (41/402 nt)ytfApseudogene, related to transcriptional regulators
Reads supporting (aligned to +/- strand):  ref base A (12/6);  new base C (120/148);  total (132/154)
Fisher's exact test for biased strand distribution p-value = 8.84e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

Only 100 of 291 total aligned reads displayed.

GGGGAGAAAATGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGAAAC  >  CP009273/4417404‑4417700
                                                                                                                                                   |                                                                                                                                                     
  ggAGAAAATGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGcttat                                                                                                                                                   <  1:47011/150‑1 (MQ=255)
    agaAAATGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCttatta                                                                                                                                                 <  1:2996552/150‑1 (MQ=255)
    agaAAATGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCttatta                                                                                                                                                 <  2:534574/150‑1 (MQ=255)
     gaAAATGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTtattat                                                                                                                                                >  1:1042446/1‑150 (MQ=255)
       aaaTGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTatat                                                                                                                                              <  1:2460859/150‑1 (MQ=255)
         aTGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTg                                                                                                                                            <  2:4026967/150‑1 (MQ=255)
              ttCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAAt                                                                                                                                       <  1:4589592/150‑1 (MQ=255)
                ccACGCCTTACACCGTACCTTCTTGTTGTGCGGTATGTCATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATcc                                                                                                                                     >  2:3196481/1‑150 (MQ=255)
                ccACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACGAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATcc                                                                                                                                     <  1:205048/150‑1 (MQ=255)
                 cACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCt                                                                                                                                    <  2:2561533/150‑1 (MQ=255)
                  aCGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCtt                                                                                                                                   <  2:915735/150‑1 (MQ=255)
                  aCGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCtt                                                                                                                                   <  1:4742386/150‑1 (MQ=255)
                     ccTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGt                                                                                                                                <  2:4710701/150‑1 (MQ=255)
                        tACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGtttt                                                                                                                             >  2:619612/1‑150 (MQ=255)
                        tACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGtttt                                                                                                                             >  1:1064184/1‑150 (MQ=255)
                         acacCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGttttt                                                                                                                            <  1:2343600/150‑1 (MQ=255)
                          cacCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTa                                                                                                                           <  2:35188/150‑1 (MQ=255)
                            ccGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTaaa                                                                                                                         <  2:1855906/150‑1 (MQ=255)
                              gTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTaaaaa                                                                                                                       >  2:2263622/1‑150 (MQ=255)
                              gTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTaaaaa                                                                                                                       >  2:4398311/1‑150 (MQ=255)
                              gTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGATAGCCAGCCACAAAAAGCCATTGAGAAAAATATTGATTTTACATGGGCTTATTATATAGCTAATCCTTGGTTTTTaaaaa                                                                                                                       >  1:1051943/1‑150 (MQ=255)
                                aCATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAAtt                                                                                                                     >  2:3979129/1‑150 (MQ=255)
                                 cATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATtg                                                                                                                    <  2:1988292/150‑1 (MQ=255)
                                  aTTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATtgt                                                                                                                   <  1:1878239/150‑1 (MQ=255)
                                  aTTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATtgt                                                                                                                   >  1:3657132/1‑150 (MQ=255)
                                   gtCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATtgtg                                                                                                                  <  2:2893660/149‑1 (MQ=255)
                                       tgttgtGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATa                                                                                                              >  1:4952514/1‑150 (MQ=255)
                                          tgtgCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTg                                                                                                           <  1:437466/150‑1 (MQ=255)
                                          tgtgCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTg                                                                                                           >  1:5118298/1‑150 (MQ=255)
                                             gCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGaaa                                                                                                        >  2:2397381/1‑150 (MQ=255)
                                                  aTGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCaa                                                                                                   >  1:3472316/1‑150 (MQ=255)
                                                   tGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCaaa                                                                                                  <  2:4350436/150‑1 (MQ=255)
                                                    gTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAg                                                                                                 <  2:803306/150‑1 (MQ=255)
                                                      aaTTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTa                                                                                               <  2:1081558/150‑1 (MQ=255)
                                                       attattCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTaa                                                                                              >  2:2125889/1‑150 (MQ=255)
                                                        ttattCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTaaa                                                                                             >  2:1375865/1‑150 (MQ=255)
                                                          attCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAgg                                                                                           >  2:2298862/1‑150 (MQ=255)
                                                          attCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAgg                                                                                           <  1:3196481/150‑1 (MQ=255)
                                                             cTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAATGTAAAGgtgt                                                                                        <  2:4295816/150‑1 (MQ=255)
                                                                  cTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCAtt                                                                                   >  1:2174637/1‑150 (MQ=255)
                                                                    gTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTAt                                                                                 <  1:1269540/150‑1 (MQ=255)
                                                                      cGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGc                                                                               >  2:4491104/1‑150 (MQ=255)
                                                                        gTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCaa                                                                             <  2:982121/150‑1 (MQ=255)
                                                                            cGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATc                                                                         >  2:2386345/1‑150 (MQ=255)
                                                                             ggTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATca                                                                        >  1:3879531/1‑150 (MQ=255)
                                                                             ggTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATca                                                                        <  1:2698298/150‑1 (MQ=255)
                                                                               tgtATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAACGTAAAGGTGTGCATTATGCAAAATcaca                                                                      <  1:1105257/150‑1 (MQ=255)
                                                                                 tATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGc                                                                    <  1:4071159/150‑1 (MQ=255)
                                                                                  aTTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCg                                                                   >  1:1332443/1‑150 (MQ=255)
                                                                                  aTTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCg                                                                   >  1:3607586/1‑150 (MQ=255)
                                                                                  aTTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCg                                                                   <  2:1341858/150‑1 (MQ=255)
                                                                                    tCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGtc                                                                 <  2:2670581/150‑1 (MQ=255)
                                                                                     cGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGtct                                                                >  1:1300982/1‑150 (MQ=255)
                                                                                      gTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTTAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGtctc                                                               <  2:1680552/150‑1 (MQ=255)
                                                                                      gTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGtctc                                                               <  2:2482278/150‑1 (MQ=255)
                                                                                       ttAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTctct                                                              >  2:655710/1‑150 (MQ=255)
                                                                                        tAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCtt                                                             >  1:3806319/1‑150 (MQ=255)
                                                                                          aTCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAATATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAg                                                           >  2:3238488/1‑150 (MQ=255)
                                                                                           tCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGa                                                          <  2:3021329/150‑1 (MQ=255)
                                                                                           tCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGa                                                          >  2:3449468/1‑150 (MQ=255)
                                                                                            cagcagGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGaa                                                         <  1:1594890/150‑1 (MQ=255)
                                                                                                agGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTa                                                     >  2:2942691/1‑150 (MQ=255)
                                                                                                     agccagccACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCAt                                                >  2:3315378/1‑150 (MQ=255)
                                                                                                     agccagccACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCAt                                                <  1:3521792/150‑1 (MQ=255)
                                                                                                     agccagccACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCAt                                                <  1:4398952/150‑1 (MQ=255)
                                                                                                      gccagccaCAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATg                                               <  1:3436454/150‑1 (MQ=255)
                                                                                                      gccagccaCAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATg                                               <  2:3093779/150‑1 (MQ=255)
                                                                                                       ccagccaCAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATgc                                              <  2:3462534/150‑1 (MQ=255)
                                                                                                        cagccaCAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATgcg                                             <  1:757277/150‑1 (MQ=255)
                                                                                                         agccaCAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATgcgc                                            <  1:3033410/150‑1 (MQ=255)
                                                                                                           ccaCAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTCCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCAc                                          <  2:4658459/150‑1 (MQ=255)
                                                                                                            cacaAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACt                                         >  1:2565949/1‑150 (MQ=255)
                                                                                                            cacaAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACt                                         >  1:4185992/1‑150 (MQ=255)
                                                                                                               aaaaaGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGc                                      >  2:64324/1‑150 (MQ=255)
                                                                                                                 aaaGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAg                                    >  1:788441/1‑150 (MQ=255)
                                                                                                                 aaaGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAg                                    >  1:1951580/1‑150 (MQ=255)
                                                                                                                   aGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGtt                                  <  1:4633037/150‑1 (MQ=255)
                                                                                                                     ccATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTaa                                >  1:812821/1‑150 (MQ=255)
                                                                                                                        ttGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATgg                             >  1:4667557/1‑150 (MQ=255)
                                                                                                                          gagaAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGaa                           >  2:1491132/1‑150 (MQ=255)
                                                                                                                           agaAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAAt                          >  2:3859452/1‑150 (MQ=255)
                                                                                                                            gaAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATc                         <  1:4860408/150‑1 (MQ=255)
                                                                                                                             aaaaTTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATcc                        <  1:2597057/150‑1 (MQ=255)
                                                                                                                             aaaaTTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATcc                        >  1:638158/1‑150 (MQ=255)
                                                                                                                             aaaaTTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATcc                        >  2:729405/1‑150 (MQ=255)
                                                                                                                               aaTTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCtt                      <  1:3635115/150‑1 (MQ=255)
                                                                                                                                 ttattGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTAt                    <  2:422693/150‑1 (MQ=255)
                                                                                                                                  tattGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTAtt                   >  1:4059981/1‑150 (MQ=255)
                                                                                                                                  tattGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTAtt                   <  1:3470487/150‑1 (MQ=255)
                                                                                                                                   attGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTg                  >  1:953390/1‑150 (MQ=255)
                                                                                                                                   attGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTg                  >  2:2200818/1‑150 (MQ=255)
                                                                                                                                    ttGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTgg                 >  1:140808/1‑150 (MQ=255)
                                                                                                                                     tGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGa                >  2:2390684/1‑150 (MQ=255)
                                                                                                                                       aTTTTACATGGGCTTATTATATTGCTAATACTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAACGTTTTGCATTATGCAAAATCACAGCGTATCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATg              <  1:2074196/150‑1 (MQ=255)
                                                                                                                                        ttttACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGa             <  1:1542937/150‑1 (MQ=255)
                                                                                                                                          ttACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAg           <  1:496662/150‑1 (MQ=255)
                                                                                                                                           tACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAgg          >  1:2218523/1‑150 (MQ=255)
                                                                                                                                             cATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGga        >  2:2028860/1‑150 (MQ=255)
                                                                                                                                                gggCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGgaaga     <  2:4892656/150‑1 (MQ=255)
                                                                                                                                                gggCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGgaaga     <  2:4120370/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: