Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 259,494 G→T 100% C75* (TGC→TGA yafW ← CP4‑6 prophage; antitoxin of the YkfI‑YafW toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273259,4940GT96.2% 564.5 / 12.3 184C75* (TGC→TGAyafWCP4‑6 prophage; antitoxin of the YkfI‑YafW toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base G (2/5);  new base T (89/88);  total (91/93)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 184 total aligned reads displayed.

GCCCGCTGAGTTATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACGCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTCAGCCAGATAATG  >  CP009273/259349‑259637
                                                                                                                                                 |                                                                                                                                               
gcccgcTGAGTTATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCa                                                                                                                                               >  1:1050437/1‑148 (MQ=255)
gcccgcTGAGTTATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCa                                                                                                                                               >  1:332866/1‑148 (MQ=255)
gcccgcTGAGTTAGTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCa                                                                                                                                               <  1:2232248/148‑1 (MQ=255)
   cgcTGAGTTATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTg                                                                                                                                            <  2:3092964/148‑1 (MQ=255)
    gcTGAGTTATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTgg                                                                                                                                           <  1:1538044/148‑1 (MQ=255)
     cTGAGTTATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTggg                                                                                                                                          >  2:1413329/1‑148 (MQ=255)
     cTGAGTTATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTggg                                                                                                                                          >  2:1413409/1‑148 (MQ=255)
      tGAGTTATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGc                                                                                                                                         >  1:2367933/1‑148 (MQ=255)
       gAGTTATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCa                                                                                                                                        >  2:2027062/1‑148 (MQ=255)
          ttATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCAtgg                                                                                                                                     >  2:1496939/1‑148 (MQ=255)
          ttATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACGCAGTGGGCAtgg                                                                                                                                     >  1:1401825/1‑148 (MQ=255)
          ttATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACGCAGTGGGCAtgg                                                                                                                                     >  1:1782839/1‑148 (MQ=255)
              tGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTggg                                                                                                                                 >  1:1446342/1‑148 (MQ=255)
                cAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTggggg                                                                                                                               <  1:2975573/148‑1 (MQ=255)
                     aaaGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAg                                                                                                                          >  1:3341978/1‑148 (MQ=255)
                     aaaGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAg                                                                                                                          <  2:1439196/148‑1 (MQ=255)
                     aaaGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAg                                                                                                                          <  2:1439549/148‑1 (MQ=255)
                      aaGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGt                                                                                                                         >  2:1799692/1‑148 (MQ=255)
                      aaGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGt                                                                                                                         <  2:1059585/148‑1 (MQ=255)
                      aaGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGt                                                                                                                         >  2:1799728/1‑148 (MQ=255)
                       aGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGtt                                                                                                                        <  1:2654043/148‑1 (MQ=255)
                        gTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTc                                                                                                                       <  1:1315056/148‑1 (MQ=255)
                        gTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTc                                                                                                                       <  1:2467180/148‑1 (MQ=255)
                          tttCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCAc                                                                                                                     >  2:428289/1‑148 (MQ=255)
                           ttCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCAcc                                                                                                                    >  1:2742159/1‑148 (MQ=255)
                               tGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCggtg                                                                                                                <  2:1568935/148‑1 (MQ=255)
                               tGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCggtg                                                                                                                >  2:3132477/1‑148 (MQ=255)
                               tGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCggtg                                                                                                                <  1:581166/148‑1 (MQ=255)
                                gTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCggtgg                                                                                                               >  1:87471/1‑148 (MQ=255)
                                gTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCggtgg                                                                                                               >  1:3378482/1‑148 (MQ=255)
                                 tttGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCggtggt                                                                                                              >  2:434415/1‑148 (MQ=255)
                                     ctctCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGc                                                                                                          >  2:2854959/1‑148 (MQ=255)
                                          gtgtAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGa                                                                                                     <  1:1678483/148‑1 (MQ=255)
                                             tAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTc                                                                                                  <  2:1069571/148‑1 (MQ=255)
                                             tAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTc                                                                                                  <  1:765119/148‑1 (MQ=255)
                                             tAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTc                                                                                                  <  1:765318/148‑1 (MQ=255)
                                               ggTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCa                                                                                                <  1:2275924/148‑1 (MQ=255)
                                                gTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCAt                                                                                               <  1:840832/148‑1 (MQ=255)
                                                gTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCAt                                                                                               >  2:301733/1‑148 (MQ=255)
                                                   aTTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTg                                                                                            <  1:3132477/148‑1 (MQ=255)
                                                   aTTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTg                                                                                            <  2:3358842/148‑1 (MQ=255)
                                                   aTTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTg                                                                                            <  1:1348578/148‑1 (MQ=255)
                                                   aTTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTg                                                                                            >  2:1162289/1‑148 (MQ=255)
                                                   aTTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTg                                                                                            >  2:2270899/1‑148 (MQ=255)
                                                    ttAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTgg                                                                                           >  1:2993425/1‑148 (MQ=255)
                                                    ttAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTgg                                                                                           <  1:2415642/148‑1 (MQ=255)
                                                       aCGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTggagg                                                                                        <  1:2817423/148‑1 (MQ=255)
                                                       aCGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTggagg                                                                                        <  1:2027110/148‑1 (MQ=255)
                                                        cGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTggagga                                                                                       >  2:1255100/1‑148 (MQ=255)
                                                            gAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACGCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATaaa                                                                                   <  1:1635312/148‑1 (MQ=255)
                                                            gAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACGCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATaaa                                                                                   <  1:1635335/148‑1 (MQ=255)
                                                             aGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAg                                                                                  >  2:68727/1‑148 (MQ=255)
                                                             aGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAg                                                                                  <  2:857599/148‑1 (MQ=255)
                                                             aGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAg                                                                                  <  1:810236/148‑1 (MQ=255)
                                                             aGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAg                                                                                  <  1:1941464/148‑1 (MQ=255)
                                                             aGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAg                                                                                  <  1:3391063/148‑1 (MQ=255)
                                                               tGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTg                                                                                <  2:1000021/148‑1 (MQ=255)
                                                                 gggTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCg                                                                              >  2:1761239/1‑148 (MQ=255)
                                                                  ggTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACGCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCgg                                                                             <  2:1782839/148‑1 (MQ=255)
                                                                  ggTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACGCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCgg                                                                             <  2:1401825/148‑1 (MQ=255)
                                                                   gTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCggg                                                                            >  1:1119073/1‑148 (MQ=255)
                                                                      aaTGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAAt                                                                         >  2:1309047/1‑148 (MQ=255)
                                                                       aTGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATg                                                                        >  2:1121742/1‑148 (MQ=255)
                                                                          gCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGttt                                                                     <  1:946365/148‑1 (MQ=255)
                                                                          gCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGttt                                                                     <  1:947912/148‑1 (MQ=255)
                                                                          gCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGttt                                                                     <  1:3503791/148‑1 (MQ=255)
                                                                               gTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATcc                                                                <  2:2403700/148‑1 (MQ=255)
                                                                               gTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATcc                                                                >  2:2071688/1‑148 (MQ=255)
                                                                                  taCGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCTTGGTGGGGGCTGAGTTCACAGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGc                                                             <  1:1882975/148‑1 (MQ=255)
                                                                                   aCGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCt                                                            >  2:1620785/1‑148 (MQ=255)
                                                                                    cGTAGAAGCAACTTCCAAGGGTGTCGGCTTCTCAGGTAAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATAGATTCAATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCta                                                           >  2:3561121/1‑147 (MQ=255)
                                                                                     gTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTc                                                          <  2:3255721/148‑1 (MQ=255)
                                                                                      tAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTcc                                                         >  2:588148/1‑148 (MQ=255)
                                                                                          cGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCa                                                     <  2:2569212/148‑1 (MQ=255)
                                                                                            cAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCAtt                                                   <  2:2603074/148‑1 (MQ=255)
                                                                                              aCTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCg                                                 <  1:129309/148‑1 (MQ=255)
                                                                                                  ccAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGt                                             <  2:3548231/148‑1 (MQ=255)
                                                                                                    aaGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCa                                           <  1:1531813/148‑1 (MQ=255)
                                                                                                     aGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGATGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCAc                                          >  2:1204271/1‑148 (MQ=255)
                                                                                                        gtgtCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGa                                       >  1:649987/1‑148 (MQ=255)
                                                                                                         tgtCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGaa                                      <  1:434415/148‑1 (MQ=255)
                                                                                                           tCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACt                                    >  1:118679/1‑148 (MQ=255)
                                                                                                           tCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACt                                    >  1:1233273/1‑148 (MQ=255)
                                                                                                             ggCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTg                                  <  1:1604415/148‑1 (MQ=255)
                                                                                                                    cAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGa                           <  1:2467716/148‑1 (MQ=255)
                                                                                                                    cAGGTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGa                           >  1:3269962/1‑148 (MQ=255)
                                                                                                                      ggTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATg                         >  1:1216984/1‑148 (MQ=255)
                                                                                                                       gTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGc                        <  1:955754/148‑1 (MQ=255)
                                                                                                                       gTTAAATCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGc                        <  1:2765180/148‑1 (MQ=255)
                                                                                                                            aTCGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGcc                   >  1:2021343/1‑148 (MQ=255)
                                                                                                                              cGTTGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCg                 >  2:236794/1‑148 (MQ=255)
                                                                                                                                ttGTGGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGt               <  1:2799565/148‑1 (MQ=255)
                                                                                                                                   tgGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTCAg            >  1:2390919/1‑148 (MQ=255)
                                                                                                                                   tgGTACAGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTCAg            >  1:2920068/1‑148 (MQ=255)
                                                                                                                                         aGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTCAGCCAGAt      <  1:289482/148‑1 (MQ=255)
                                                                                                                                         aGGGTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTCAGCCAGAt      <  2:3341978/148‑1 (MQ=255)
                                                                                                                                           ggTAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTCAGCCAGATaa    <  1:3468184/148‑1 (MQ=255)
                                                                                                                                             tAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTCAGCCAGATAATg  <  1:1413329/148‑1 (MQ=255)
                                                                                                                                             tAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTCAGCCAGATAATg  <  1:1413409/148‑1 (MQ=255)
                                                                                                                                             tAACTCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTCAGCCAGATAATg  >  2:1691251/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
GCCCGCTGAGTTATTGCAGGTAAAGTTTTCATGTTTGCTCTCGTGTAGGTAATTAACGCTGAGTGGGGTAAATGGCGATGTATACGTAGCCGCAACTGCCAAGGGTGTCGGCTTCGCAGGTTAAATCGTTGTGGTACAGGGTAACGCAGTGGGCATGGTGGGGGCTGAGTTCACCGGTGGTCAGCATCGATTCCATCTGGAGGATAAAGTGCGGGAATGTTTCATCCAGCTTCCGGCATTCGATGTCACTGAACTTGCCGGTGATGCTGGCCCGGTCAGCCAGATAATG  >  CP009273/259349‑259637

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: