Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,417,551 A→C 100% 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,5510AC96.5% 996.0 / 10.5 256pseudogene (41/402 nt)ytfApseudogene, related to transcriptional regulators
Reads supporting (aligned to +/- strand):  ref base A (2/7);  new base C (119/128);  total (121/135)
Fisher's exact test for biased strand distribution p-value = 1.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 264 total aligned reads displayed.

GGGGAGAAAATGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGA  >  CP009273/4417404‑4417697
                                                                                                                                                   |                                                                                                                                                  
ggggAGAAAATGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGctt                                                                                                                                                  <  1:3728584/150‑1 (MQ=255)
 gggAGAAAATGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGctta                                                                                                                                                 <  1:1320638/150‑1 (MQ=255)
   gagaAAATGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGcttatt                                                                                                                                               <  2:1771457/150‑1 (MQ=255)
    agaAAATGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCttatta                                                                                                                                              <  1:4462565/150‑1 (MQ=255)
     gaAAATGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTtattat                                                                                                                                             >  2:2390156/1‑150 (MQ=255)
     gaAAATGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTtattat                                                                                                                                             <  2:1632191/150‑1 (MQ=255)
     gaAAATGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTtattat                                                                                                                                             >  1:2468574/1‑150 (MQ=255)
        aaTGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATAtt                                                                                                                                          >  2:2489413/1‑150 (MQ=255)
           gACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCt                                                                                                                                       <  1:3698278/150‑1 (MQ=255)
            aCTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTa                                                                                                                                      >  2:4168340/1‑150 (MQ=255)
              ttCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAAt                                                                                                                                    <  2:4016690/150‑1 (MQ=255)
                 cACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCt                                                                                                                                 <  2:3222276/150‑1 (MQ=255)
                  aCGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCtt                                                                                                                                <  2:4283008/150‑1 (MQ=255)
                   cGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTg                                                                                                                               <  2:644437/150‑1 (MQ=255)
                    gCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCGACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTgg                                                                                                                              <  2:3548795/150‑1 (MQ=255)
                     ccTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGt                                                                                                                             >  1:4030207/1‑150 (MQ=255)
                         acacCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGttttt                                                                                                                         <  2:3574218/150‑1 (MQ=255)
                          cacCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTa                                                                                                                        <  1:2675451/150‑1 (MQ=255)
                           acCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTaa                                                                                                                       >  1:3447409/1‑150 (MQ=255)
                            ccGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTaaa                                                                                                                      >  1:1866337/1‑150 (MQ=255)
                                  aTTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATtgt                                                                                                                <  2:736602/150‑1 (MQ=255)
                                   ttCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATtgtg                                                                                                               >  2:530657/1‑150 (MQ=255)
                                   ttCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATtgtg                                                                                                               <  1:2229733/150‑1 (MQ=255)
                                    tCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGc                                                                                                              <  2:87833/150‑1 (MQ=255)
                                        gttgtGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAg                                                                                                          <  1:2599012/150‑1 (MQ=255)
                                          tgtgCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTg                                                                                                        >  1:2272906/1‑150 (MQ=255)
                                                gTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATc                                                                                                  <  1:372435/150‑1 (MQ=255)
                                                gTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATc                                                                                                  <  1:446633/150‑1 (MQ=255)
                                                  aTGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCaa                                                                                                >  2:2790681/1‑150 (MQ=255)
                                                  aTGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCaa                                                                                                >  2:2900519/1‑150 (MQ=255)
                                                    gTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAg                                                                                              <  2:317579/150‑1 (MQ=255)
                                                     tAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTNTTAAAAATTGTGCATAGTGAAAATCAAAGt                                                                                             <  1:4285287/150‑1 (MQ=255)
                                                        ttattCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTaaa                                                                                          >  1:4545945/1‑150 (MQ=255)
                                                        ttattCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTaaa                                                                                          <  2:1306654/150‑1 (MQ=255)
                                                         tattCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAg                                                                                         <  2:2816702/150‑1 (MQ=255)
                                                          attCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAgg                                                                                        <  2:3079372/150‑1 (MQ=255)
                                                           ttCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGgt                                                                                       <  1:2132580/150‑1 (MQ=255)
                                                            tCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGgtg                                                                                      >  2:1023650/1‑150 (MQ=255)
                                                            tCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGgtg                                                                                      >  2:3295987/1‑150 (MQ=255)
                                                            tCTACGCTGTCGGTATCGGTGTATTAGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTTAAATTTGTGCATAGTGAAAATCAAAGTAAAGgtg                                                                                      <  2:1374299/150‑1 (MQ=255)
                                                               aCGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGc                                                                                   >  2:483113/1‑150 (MQ=255)
                                                                cGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCa                                                                                  >  2:3812061/1‑150 (MQ=255)
                                                                 gCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCAt                                                                                 <  1:1581550/150‑1 (MQ=255)
                                                                 gCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCAt                                                                                 >  2:2993029/1‑150 (MQ=255)
                                                                 gCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCAt                                                                                 <  2:2035800/150‑1 (MQ=255)
                                                                  cTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCAtt                                                                                <  2:628398/150‑1 (MQ=255)
                                                                  cTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCAtt                                                                                <  2:4570015/150‑1 (MQ=255)
                                                                    gTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTAt                                                                              <  1:3374514/150‑1 (MQ=255)
                                                                     tCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATg                                                                             <  1:2481351/150‑1 (MQ=255)
                                                                      cGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGc                                                                            <  1:183846/150‑1 (MQ=255)
                                                                      cGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGc                                                                            >  1:225206/1‑150 (MQ=255)
                                                                       ggTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGTATTATGCa                                                                           >  1:4210393/1‑150 (MQ=255)
                                                                       ggTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCa                                                                           >  1:1225872/1‑150 (MQ=255)
                                                                       ggTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCa                                                                           >  1:1727912/1‑150 (MQ=255)
                                                                        gTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAATATTATTGATTTTACATGTGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGTAAATCAAAGTAAAGGTGTGCATTATGCaa                                                                          <  1:3982961/150‑1 (MQ=255)
                                                                        gTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGACTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCaa                                                                          >  1:84408/1‑150 (MQ=255)
                                                                         tATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCaaa                                                                         >  2:2239771/1‑150 (MQ=255)
                                                                              gtgtATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATcac                                                                    <  1:574842/150‑1 (MQ=255)
                                                                               tgtATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATcaca                                                                   >  2:4103451/1‑150 (MQ=255)
                                                                               tgtATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATcaca                                                                   <  1:3614540/150‑1 (MQ=255)
                                                                                      gTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGtctc                                                            >  2:743697/1‑150 (MQ=255)
                                                                                        tAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCtt                                                          <  2:3621545/150‑1 (MQ=255)
                                                                                        tAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCtt                                                          >  1:1054453/1‑150 (MQ=255)
                                                                                        tAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCtt                                                          >  1:3969217/1‑150 (MQ=255)
                                                                                               cagGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAAtt                                                   <  1:90623/150‑1 (MQ=255)
                                                                                                     agccagccACAAAAAGACATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTACAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCAt                                             <  2:1857014/150‑1 (MQ=255)
                                                                                                      gccagccaCAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAATGTAAAGGTGTGCAGTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATg                                            >  1:3088829/1‑150 (MQ=255)
                                                                                                        cagccaCAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATgcg                                          >  1:3505437/1‑150 (MQ=255)
                                                                                                         agccaCAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATgcgc                                         >  2:1478627/1‑150 (MQ=255)
                                                                                                            cacaAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACt                                      >  1:1776580/1‑150 (MQ=255)
                                                                                                             acaAAAAGCCATTTAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACtt                                     <  2:3615014/150‑1 (MQ=255)
                                                                                                             acaAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACtt                                     <  2:3858270/150‑1 (MQ=255)
                                                                                                             acaAAAAGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGGAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACtt                                     >  2:368235/1‑150 (MQ=255)
                                                                                                                aaaaGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCa                                  >  2:412881/1‑150 (MQ=255)
                                                                                                                 aaaGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAg                                 <  2:619186/150‑1 (MQ=255)
                                                                                                                  aaGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGt                                >  1:2330277/1‑150 (MQ=255)
                                                                                                                   aGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGtt                               >  1:1013485/1‑150 (MQ=255)
                                                                                                                   aGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGtt                               <  1:49928/150‑1 (MQ=255)
                                                                                                                   aGCCATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGtt                               >  2:3365120/1‑150 (MQ=255)
                                                                                                                     ccATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCAAATCCTTGGTTTTTAAAAATTGTGCATAGTTAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAACGTCTCTTAGAAATTAACCATTCGCACTTGCAGTTaa                             >  2:1532763/1‑150 (MQ=255)
                                                                                                                      cATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAAt                            >  2:1679585/1‑150 (MQ=255)
                                                                                                                      cATTGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTAGCAGTTAAt                            >  2:1649910/1‑150 (MQ=255)
                                                                                                                         tGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGa                         <  2:2468574/150‑1 (MQ=255)
                                                                                                                         tGAGAAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGa                         >  1:3615045/1‑150 (MQ=255)
                                                                                                                           agaAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAAt                       >  2:2684486/1‑150 (MQ=255)
                                                                                                                            gaAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAGCCATGCGCACTTGCAGTTAATGGAATc                      <  1:726563/150‑1 (MQ=255)
                                                                                                                            gaAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATc                      <  2:3019754/150‑1 (MQ=255)
                                                                                                                            gaAAATTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATc                      <  2:234792/150‑1 (MQ=255)
                                                                                                                              aaaTTATTGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCt                    <  1:4246594/150‑1 (MQ=255)
                                                                                                                                attattGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTa                  >  2:173014/1‑150 (MQ=255)
                                                                                                                                  tattGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTAtt                <  1:2685876/150‑1 (MQ=255)
                                                                                                                                   attGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGTGCACTTGCAGTTAATGGAATCCTTATTg               <  2:1918736/150‑1 (MQ=255)
                                                                                                                                    ttGATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTgg              >  1:1928633/1‑150 (MQ=255)
                                                                                                                                      gATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGAt            <  2:2981747/150‑1 (MQ=255)
                                                                                                                                      gATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGAt            >  1:3113252/1‑150 (MQ=255)
                                                                                                                                      gATTTTACATGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGAGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAAAGGAATCCTTATTTGAt            >  1:3231939/1‑150 (MQ=255)
                                                                                                                                              aTGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGgaa    <  2:439584/150‑1 (MQ=255)
                                                                                                                                              aTGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGgaa    >  1:2082608/1‑150 (MQ=255)
                                                                                                                                               tGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGgaag   >  2:1532919/1‑150 (MQ=255)
                                                                                                                                               tGGGCTTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGgaag   <  2:4467012/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GGGGAGAAAATGACTTCCACGCCTTACACCGTACATTCTTGTTGTGCGGTATGTAATTATTCTACGCTGTCGGTATCGGTGTATTCGTTAATCAGCAGGTTAGCCAGCCACAAAAAGCCATTGAGAAAATTATTGATTTTACATGGGATTATTATATTGCTAATCCTTGGTTTTTAAAAATTGTGCATAGTGAAAATCAAAGTAAAGGTGTGCATTATGCAAAATCACAGCGTCTCTTAGAAATTAACCATGCGCACTTGCAGTTAATGGAATCCTTATTGGATGAAGGGAAGA  >  CP009273/4417404‑4417697

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: