Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,864,752 C→A 47.6% S37R (AGC→AGA yeaI → putative membrane‑anchored diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,864,7520CA47.6% 173.1 / inf 275S37R (AGC→AGAyeaIputative membrane‑anchored diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (61/83);  new base A (64/67);  total (125/150)
Fisher's exact test for biased strand distribution p-value = 3.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 277 total aligned reads displayed.

TTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTCCATGGGGTTATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCAC  >  CP009273/1864603‑1864900
                                                                                                                                                     |                                                                                                                                                    
tttACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTCCATGGGATTATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGc                                                                                                                                                      >  1:4235097/1‑150 (MQ=255)
        aTTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTCCATGGGGTTATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACAt                                                                                                                                              <  1:634880/150‑1 (MQ=255)
         tttAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTCCATGGGGTTATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACAtt                                                                                                                                             >  1:690759/1‑150 (MQ=255)
          ttAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTCCATGGGGTTATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGAACTCACATTg                                                                                                                                            >  1:1982808/1‑150 (MQ=255)
           tAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTCCATGGGGTTATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGAACTCACATTGc                                                                                                                                           >  1:1822039/1‑150 (MQ=255)
            aGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTCCATGGGGTTATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGAACTCACATTGct                                                                                                                                          >  2:3935644/1‑150 (MQ=255)
              ttCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTCCATGGGGTTATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGAACTCACATTGctct                                                                                                                                        >  2:1527802/1‑150 (MQ=255)
                  aGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTCCATGGGGTTATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGAACTCACATTGCTCTCCtt                                                                                                                                    >  1:150134/1‑150 (MQ=255)
                        cAATCATTCGGAAAAATGATTCAGTCAACACGTATTTCCATGGGGTTATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGAACTCACATTGCTCTCCTTTTTATg                                                                                                                              >  1:3325365/1‑150 (MQ=255)
                         aaTCATTCGGAAAAATGATTCAGTCAACACGTATTTCCATGGGGTTATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGAACTCACATTGCTCTCCTTTTTATgg                                                                                                                             <  2:980186/150‑1 (MQ=255)
                            cATTCGGAAAAATGATTCAGTCAACACGTATTTCCATGGGGTTATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGttt                                                                                                                          <  2:3076078/150‑1 (MQ=255)
                             aTTCGGAAAAATGATTCAGTCAACACGTATTTCCATGGGGTTATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTc                                                                                                                         >  2:3829562/1‑150 (MQ=255)
                                cGGAAAAATGATTCAGTCAACACGTATTTCCATGGGGTTATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTAt                                                                                                                      <  1:3371790/150‑1 (MQ=255)
                                cGGAAAAATGATTCAGTCAACACGTATTTCCATGGGGTTATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTAt                                                                                                                      <  1:3390118/150‑1 (MQ=255)
                                cGGAAAAATGATTCAGTCAACACGTATTTCCATGGGGTTATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGAACTCACATTGCTCTCCTTTTTATGGTTTCTAt                                                                                                                      <  2:776543/150‑1 (MQ=255)
                                   aaaaaTGATTCAGTCAACACGTATTTCCATGGGGTTATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATggg                                                                                                                   >  1:2760700/1‑150 (MQ=255)
                                    aaaaTGATTCAGTCAACACGTATTTCCATGGGGTTATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGt                                                                                                                  >  1:977553/1‑150 (MQ=255)
                                    aaaaTGATTCAGTCAACACGTATTTCCATGGGGTTATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCATTTTTATGGTTTCTATGGGt                                                                                                                  >  2:3957251/1‑150 (MQ=255)
                                             cAGTCAACACGTATTTCCATGGGGTTATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGAACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACaaaaa                                                                                                         >  1:3843172/1‑150 (MQ=255)
                                              aGTCAACACGTATTTCCATGGGGTTATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAAt                                                                                                        >  2:91956/1‑150 (MQ=255)
                                                tCAACACGTATTTCCATGGGGTTATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTATGCCATCAATAAAATCCAGAACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTa                                                                                                      >  2:1767710/1‑150 (MQ=255)
                                                  aaCACGTATTTCCATGGGGTTATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTaaa                                                                                                    <  1:3544067/150‑1 (MQ=255)
                                                  aaCACGTATTTCCATGGGGTTATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTaaa                                                                                                    <  1:2587096/150‑1 (MQ=255)
                                                   acacGTATTTCCATGGGGTTATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTaaaa                                                                                                   <  2:2115347/150‑1 (MQ=255)
                                                      cGTATTTCCATGGGGTTATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCt                                                                                                <  1:680805/150‑1 (MQ=255)
                                                           ttCCATGGGGTTATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAg                                                                                           <  2:1563507/150‑1 (MQ=255)
                                                              cATGGGGTTATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCtt                                                                                        <  1:535049/150‑1 (MQ=255)
                                                               aTGGGGTTATTCTTTAAAAATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCttt                                                                                       <  2:128965/150‑1 (MQ=255)
                                                                 ggggTTATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCttttt                                                                                     >  2:4197187/1‑150 (MQ=255)
                                                                     ttATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGAACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTAtt                                                                                 >  1:2994017/1‑150 (MQ=255)
                                                                         tCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGAACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCa                                                                             <  1:2924433/150‑1 (MQ=255)
                                                                         tCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGAACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCa                                                                             <  2:2594006/150‑1 (MQ=255)
                                                                          cTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGAACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAg                                                                            <  1:1884828/150‑1 (MQ=255)
                                                                           tttAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGt                                                                           <  2:2170251/150‑1 (MQ=255)
                                                                           tttAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGAACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGt                                                                           <  2:2597556/150‑1 (MQ=255)
                                                                            ttAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGAACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGtt                                                                          <  1:4469921/150‑1 (MQ=255)
                                                                             tAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGAACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGttt                                                                         >  1:1443900/1‑150 (MQ=255)
                                                                              aaaTATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTa                                                                        >  2:2300894/1‑150 (MQ=255)
                                                                               aaTATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTAt                                                                       <  1:907093/150‑1 (MQ=255)
                                                                                atatTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGAACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTAtt                                                                      <  2:3150309/150‑1 (MQ=255)
                                                                                  atTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGc                                                                    <  2:4628023/150‑1 (MQ=255)
                                                                                  atTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGAACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGc                                                                    <  1:849322/150‑1 (MQ=255)
                                                                                     ttttATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACAGTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGa                                                                 <  2:4293862/150‑1 (MQ=255)
                                                                                      tttATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGAc                                                                >  2:4493983/1‑150 (MQ=255)
                                                                                      tttATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGAACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGAc                                                                <  2:3325365/150‑1 (MQ=255)
                                                                                      tttATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGAACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGAc                                                                <  2:1006237/150‑1 (MQ=255)
                                                                                        tATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGt                                                              >  1:1196295/1‑150 (MQ=255)
                                                                                        tATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGt                                                              >  1:3292173/1‑150 (MQ=255)
                                                                                         aTCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGtt                                                             <  2:2067867/150‑1 (MQ=255)
                                                                                         aTCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGAACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGtt                                                             <  2:1255441/150‑1 (MQ=255)
                                                                                          tCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTa                                                            <  1:2052613/150‑1 (MQ=255)
                                                                                           cGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGAACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTaa                                                           >  1:4499885/1‑150 (MQ=255)
                                                                                            gTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGAACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAat                                                          >  1:4465361/1‑150 (MQ=255)
                                                                                               aaCGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAatttt                                                       <  1:1688483/150‑1 (MQ=255)
                                                                                                    aaaTTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAattttatttt                                                  <  2:1370708/150‑1 (MQ=255)
                                                                                                    aaaTTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGAACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAattttatttt                                                  <  2:750579/150‑1 (MQ=255)
                                                                                                    aaaTTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGAACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAattttatttt                                                  >  1:2479217/1‑150 (MQ=255)
                                                                                                     aaTTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTg                                                 <  2:2795342/150‑1 (MQ=255)
                                                                                                     aaTTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGAACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTg                                                 >  2:569473/1‑150 (MQ=255)
                                                                                                     aaTTGATCCTGGTCAAAACTATATATCACTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTg                                                 <  2:3291741/150‑1 (MQ=255)
                                                                                                       ttGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGAACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCa                                               <  2:3757590/150‑1 (MQ=255)
                                                                                                       ttGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGAACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCa                                               >  1:3995207/1‑150 (MQ=255)
                                                                                                        tGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGAACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCAt                                              >  1:417793/1‑150 (MQ=255)
                                                                                                            ccTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATtact                                          >  2:1211578/1‑150 (MQ=255)
                                                                                                            ccTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATtact                                          >  1:3251443/1‑150 (MQ=255)
                                                                                                            ccTGGTCAAAACTATATATCTCTGCCAACAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAAATAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTATGCATtaat                                          >  1:2669242/1‑148 (MQ=255)
                                                                                                              tGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATtactac                                        >  2:192583/1‑150 (MQ=255)
                                                                                                              tGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGAACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATtactac                                        >  2:2371931/1‑150 (MQ=255)
                                                                                                              tGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGAACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATtactac                                        <  2:2455896/150‑1 (MQ=255)
                                                                                                               ggTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTAcc                                       >  2:3811467/1‑150 (MQ=255)
                                                                                                                 tCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCtt                                     <  1:1922668/150‑1 (MQ=255)
                                                                                                                     aaCTATATATCTCTGCCATCAATAAATTCCAGAACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTtatat                                 >  2:1114697/1‑150 (MQ=255)
                                                                                                                      aCTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTtatata                                <  1:1638926/150‑1 (MQ=255)
                                                                                                                      aCTATATATCTCTGCCATCAATAAAATCCAGAACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTtatata                                <  1:4140750/150‑1 (MQ=255)
                                                                                                                        tatataTCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAAc                              <  1:228845/150‑1 (MQ=255)
                                                                                                                         atatatCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACg                             <  2:1730884/150‑1 (MQ=255)
                                                                                                                          tatatCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGa                            >  2:1807382/1‑150 (MQ=255)
                                                                                                                           atatCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGaa                           <  1:2198037/150‑1 (MQ=255)
                                                                                                                             atCTCTGCCATCAATAAAATCCAGAACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGaaaa                         <  1:1017524/150‑1 (MQ=255)
                                                                                                                               ctctGCCATCAATAAAATCCAGAACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAaca                       <  2:584057/150‑1 (MQ=255)
                                                                                                                               ctctGCCATCAATAAAATCCAGAACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAaca                       <  2:2896655/150‑1 (MQ=255)
                                                                                                                                tctGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAacac                      >  1:3626038/1‑150 (MQ=255)
                                                                                                                                tctGCCATCAATAAAATCCAGAACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAacac                      >  2:2285401/1‑150 (MQ=255)
                                                                                                                                tctGCCATCAATAAAATCCAGAACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAacac                      <  2:4254376/150‑1 (MQ=255)
                                                                                                                                  tGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACaa                    >  1:1741946/1‑150 (MQ=255)
                                                                                                                                   gCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACaaa                   >  1:3066748/1‑150 (MQ=255)
                                                                                                                                     cATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATg                 >  2:3517552/1‑150 (MQ=255)
                                                                                                                                        cAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTaa              >  2:2707360/1‑150 (MQ=255)
                                                                                                                                        cAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTaa              >  1:1131018/1‑150 (MQ=255)
                                                                                                                                         aaTAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTaaa             <  1:1000590/150‑1 (MQ=255)
                                                                                                                                          aTAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTaaaa            <  2:3184536/150‑1 (MQ=255)
                                                                                                                                           tAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAAc           <  2:3716495/150‑1 (MQ=255)
                                                                                                                                            aaaaTCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACt          <  1:1330381/150‑1 (MQ=255)
                                                                                                                                            aaaaTCCAGAACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACt          <  2:2657658/150‑1 (MQ=255)
                                                                                                                                            aaaaTCCAGAACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACt          <  2:3612877/150‑1 (MQ=255)
                                                                                                                                              aaTCCAGAACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCa        <  1:2109451/150‑1 (MQ=255)
                                                                                                                                              aaTCCAGAACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCa        <  2:860391/150‑1 (MQ=255)
                                                                                                                                               aTCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCAt       <  1:4592814/150‑1 (MQ=255)
                                                                                                                                                  cAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATccc    >  1:4580264/1‑150 (MQ=255)
                                                                                                                                                    gCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCcac  >  1:2087183/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTCCATGGGGTTATTCTTTAAATATTTTTTATCGTTAACGAAAATTGATCCTGGTCAAAACTATATATCTCTGCCATCAATAAAATCCAGCACTCACATTGCTCTCCTTTTTATGGTTTCTATGGGTACACAAAAATTAAAAGCTCAAAGCTTTTTTATTTTCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCAC  >  CP009273/1864603‑1864900

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: