Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,417,454 C→A 5.3% pseudogene (497/500 nt) ydaW → Rac prophage; pseudogene, DNA‑binding protein family;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,417,4540CA 5.3% 729.9 / 11.9 195pseudogene (497/500 nt)ydaWRac prophage; pseudogene, DNA‑binding protein family;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base C (102/81);  new base A (5/5);  total (107/88)
Fisher's exact test for biased strand distribution p-value = 7.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

Only 100 of 199 total aligned reads displayed.

TACATTGTTTGGTACCACATCACGCAAAGTGGCTTCAACGCTGGCAATGGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAG  >  CP009273/1417308‑1417601
                                                                                                                                                  |                                                                                                                                                   
tACATTGTTTGGTACCACATCACGCAAAGTGGCTTCAACGCTGGCAATGGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCg                                                                                                                                                    >  2:619952/1‑148 (MQ=255)
  cATTGTTTGGTACCACATCACGCAAAGTGGCTTCAACGCTGGCAATGGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGt                                                                                                                                                  <  2:213984/148‑1 (MQ=255)
        ttGGTACCACATCACGCAAAGTGGCTTCAACGCTGGCAATGGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTcc                                                                                                                                            >  2:587619/1‑148 (MQ=255)
               cacaTCACGCAAAGTGGCTTCAACGCTGGCAATGGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTgg                                                                                                                                     >  1:322750/1‑148 (MQ=255)
                   tCACGCAAAGTGGCTTCAACGCTGGCAATGGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTAt                                                                                                                                 >  1:1482907/1‑148 (MQ=255)
                   tCACGCAAAGTGGCTTCAACGCTGGCAATGGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTAt                                                                                                                                 <  1:2145561/148‑1 (MQ=255)
                   tCACGCAAAGTGGCTTCAACGCTGGCAATGGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTAt                                                                                                                                 >  1:2041746/1‑148 (MQ=255)
                    cACGCAAAGTGGCTTCAACGCTGGCAATGGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATc                                                                                                                                >  2:1600363/1‑148 (MQ=255)
                     aCGCAAAGTGGCTTCAACGCTGGCAATGGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCt                                                                                                                               >  2:80987/1‑148 (MQ=255)
                       gCAAAGTGGCTTCAACGCTGGCAATGGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGc                                                                                                                             >  2:137155/1‑148 (MQ=255)
                        cAAAGTGGCTTCAACGCTGGCAATGGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATAGGTATCTCCGCTCTGGTTATCTGCa                                                                                                                            <  2:1721810/148‑1 (MQ=255)
                         aaaGTGGCTTCAACGCTGGCAATGGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCAt                                                                                                                           >  1:1669351/1‑148 (MQ=255)
                           aGTGGCTTCAACGCTGGCAATGGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTAACAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCAtcg                                                                                                                         >  1:302489/1‑148 (MQ=255)
                               gCTTCAACGCTGGCAATGGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATcgtcgt                                                                                                                     >  2:2402139/1‑148 (MQ=255)
                               gCTTCAACGCTGGCAATGGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATcgtcgt                                                                                                                     >  1:597944/1‑148 (MQ=255)
                               gCTTCAACGCTGGCAATGGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATcgtcgt                                                                                                                     <  1:440920/148‑1 (MQ=255)
                                cTTCAACGCTGGCAATGGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATAGGTATCTCCGCTCTGGTTATCTGCATCgtcgtc                                                                                                                    >  2:1932365/1‑148 (MQ=255)
                                cTTCAACGCTGGCAATGGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATAGGTATCTCCGCTCTGGTTATCTGCATCgtcgtc                                                                                                                    >  2:2946517/1‑148 (MQ=255)
                                   cAACGCTGGCAATGGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGc                                                                                                                 <  2:431446/148‑1 (MQ=255)
                                      cGCTGGCAATGGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTg                                                                                                              <  1:3038698/148‑1 (MQ=255)
                                      cGCTGGCAATGGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTg                                                                                                              >  1:1336025/1‑148 (MQ=255)
                                      cGCTGGCAATGGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTg                                                                                                              >  1:939553/1‑148 (MQ=255)
                                        cTGGCAATGGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATAGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTc                                                                                                            <  2:3226490/148‑1 (MQ=255)
                                           gcaatggcaatCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATg                                                                                                         >  1:2731626/1‑148 (MQ=255)
                                            caatggcaatCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATgg                                                                                                        <  1:496786/148‑1 (MQ=255)
                                              atggcaatCAGCATAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAATCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATgggg                                                                                                      <  1:1535870/148‑1 (MQ=255)
                                              atggcaatCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATgggg                                                                                                      <  1:250908/148‑1 (MQ=255)
                                              atggcaatCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATgggg                                                                                                      >  2:2811287/1‑148 (MQ=255)
                                              atggcaatCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATgggg                                                                                                      <  1:1630259/148‑1 (MQ=255)
                                                ggcaatCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGgtg                                                                                                    >  2:1850129/1‑148 (MQ=255)
                                                 gcaatCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGgtgt                                                                                                   >  2:1106800/1‑148 (MQ=255)
                                                   aatCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTa                                                                                                 <  1:612510/148‑1 (MQ=255)
                                                   aatCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTa                                                                                                 >  2:2079112/1‑148 (MQ=255)
                                                    atCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTaa                                                                                                >  1:2633612/1‑148 (MQ=255)
                                                      cAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATc                                                                                              <  2:228780/148‑1 (MQ=255)
                                                       aGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCa                                                                                             <  2:616840/148‑1 (MQ=255)
                                                        gCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCAt                                                                                            <  2:148891/148‑1 (MQ=255)
                                                          aaaGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTa                                                                                          <  2:1995776/148‑1 (MQ=255)
                                                            aGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTAcc                                                                                        >  1:2403918/1‑148 (MQ=255)
                                                            aGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTAcc                                                                                        >  1:6800/1‑148 (MQ=255)
                                                             ggTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCg                                                                                       >  1:2556949/1‑148 (MQ=255)
                                                              gtcgtcTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGt                                                                                      <  2:932200/148‑1 (MQ=255)
                                                               tcgtcTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTg                                                                                     >  2:3282752/1‑148 (MQ=255)
                                                               tcgtcTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTg                                                                                     <  1:2256689/148‑1 (MQ=255)
                                                                  tcTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATa                                                                                  >  1:1277805/1‑148 (MQ=255)
                                                                  tcTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATa                                                                                  >  1:1277569/1‑148 (MQ=255)
                                                                  tcTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATa                                                                                  <  1:1187303/148‑1 (MQ=255)
                                                                  tcTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATa                                                                                  <  2:1471915/148‑1 (MQ=255)
                                                                   cTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATGGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATaa                                                                                 <  2:2448098/148‑1 (MQ=255)
                                                                   cTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATaa                                                                                 <  2:718271/148‑1 (MQ=255)
                                                                   cTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATaa                                                                                 >  2:1188708/1‑148 (MQ=255)
                                                                   cTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCAGTACCGTGATaa                                                                                 <  2:3046616/148‑1 (MQ=255)
                                                                     gATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACg                                                                               >  2:810680/1‑148 (MQ=255)
                                                                     gATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATAGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACg                                                                               <  1:3119152/148‑1 (MQ=255)
                                                                       ttCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGcc                                                                             >  2:1097768/1‑148 (MQ=255)
                                                                        tCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCa                                                                            <  1:1429073/148‑1 (MQ=255)
                                                                        tCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCa                                                                            >  1:2202057/1‑148 (MQ=255)
                                                                           cgTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATAGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGCTAATCATTACCGTGATAACGCCATCg                                                                         >  2:2561013/1‑148 (MQ=255)
                                                                              aaaTCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCt                                                                      <  2:13381/148‑1 (MQ=255)
                                                                                 tCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACa                                                                   >  2:128995/1‑148 (MQ=255)
                                                                                  cAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACaa                                                                  >  2:1437569/1‑148 (MQ=255)
                                                                                  cAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACaa                                                                  <  1:2079112/148‑1 (MQ=255)
                                                                                    gggCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAg                                                                >  1:1548482/1‑148 (MQ=255)
                                                                                    gggCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAg                                                                <  2:855425/148‑1 (MQ=255)
                                                                                     ggCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGa                                                               >  1:1351583/1‑148 (MQ=255)
                                                                                     ggCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGa                                                               <  2:1669487/148‑1 (MQ=255)
                                                                                         gTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAg                                                           <  2:719678/148‑1 (MQ=255)
                                                                                           aaaTTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAgcg                                                         >  2:2524310/1‑148 (MQ=255)
                                                                                            aaTTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAgcgc                                                        >  1:1163617/1‑148 (MQ=255)
                                                                                             aTTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGcgcg                                                       >  2:3266544/1‑148 (MQ=255)
                                                                                               ttCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGAc                                                     >  1:1006720/1‑148 (MQ=255)
                                                                                                   ttACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACaaaa                                                 <  2:2887133/148‑1 (MQ=255)
                                                                                                   ttACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACaaaa                                                 <  2:342536/148‑1 (MQ=255)
                                                                                                    tACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACaaaaa                                                <  1:137155/148‑1 (MQ=255)
                                                                                                      cTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATg                                              <  1:1803044/148‑1 (MQ=255)
                                                                                                        gCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGcc                                            >  2:2652614/1‑148 (MQ=255)
                                                                                                         cATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCa                                           <  1:3183822/148‑1 (MQ=255)
                                                                                                            aCCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCagag                                        >  2:1614457/1‑148 (MQ=255)
                                                                                                            aCCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCagag                                        >  2:3217056/1‑148 (MQ=255)
                                                                                                             ccGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAgaga                                       >  1:2463915/1‑148 (MQ=255)
                                                                                                                  cGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGa                                  <  2:2309121/148‑1 (MQ=255)
                                                                                                                    aTAATTTACCAGCAGTGCCGAAAGCAGCATCGGTATCTCCGCTTTGGTTATCTGCATCGTCGTCTGGCTGTCATGGGGTGTTTATGATTACCGTGATAACGCCATCTCCTACAAAGAACAGCGGTACAAAAATGCCAGAGAACTGatg                                <  2:741742/148‑3 (MQ=25)
                                                                                                                     tAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCTTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGc                               <  2:379142/148‑1 (MQ=255)
                                                                                                                       aTTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTg                             >  1:388345/1‑148 (MQ=255)
                                                                                                                       aTTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTg                             >  1:2362737/1‑148 (MQ=255)
                                                                                                                                cagAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAAcgcg                    <  2:1405618/148‑1 (MQ=255)
                                                                                                                                  gagCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGc                  >  1:491746/1‑148 (MQ=255)
                                                                                                                                  gagCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGc                  >  2:204321/1‑148 (MQ=255)
                                                                                                                                   atcagAAAGCAGCATCTGTATCTCCGCTCTGTTTATCTTAATAGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAAAAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCa                 <  1:293291/144‑1 (MQ=18)
                                                                                                                                   agCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCa                 <  2:266597/148‑1 (MQ=255)
                                                                                                                                    gCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGGTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCCGAGAACTGAAGCTGGCGAACGCGGCac                >  1:629245/1‑147 (MQ=255)
                                                                                                                                      cGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAAtt              <  1:2689050/148‑1 (MQ=255)
                                                                                                                                      cGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAAtt              <  2:1163617/148‑1 (MQ=255)
                                                                                                                                         aaGCAGCATAGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACt           >  1:1764075/1‑148 (MQ=255)
                                                                                                                                            cagcaTCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTgag        <  2:2693826/148‑1 (MQ=255)
                                                                                                                                              gcaTCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGagat      <  1:1671867/148‑1 (MQ=255)
                                                                                                                                               caTCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGagatg     >  2:641105/1‑148 (MQ=255)
                                                                                                                                                 tCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGagatgca   <  1:1128004/148‑1 (MQ=255)
                                                                                                                                                  cGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGagatgcag  <  2:3187680/148‑1 (MQ=255)
                                                                                                                                                  cGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGagatgcag  >  1:2898975/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TACATTGTTTGGTACCACATCACGCAAAGTGGCTTCAACGCTGGCAATGGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAG  >  CP009273/1417308‑1417601

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: