Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 290,012 C→T 100% D167N (GAT→AAT)  yagM ← CP4‑6 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273290,0120CT97.4% 891.1 / 9.4 273D167N (GAT→AAT) yagMCP4‑6 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (5/2);  new base T (134/132);  total (139/134)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 273 total aligned reads displayed.

GACAAGAGTAGAGCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATCGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGTTTTCCAGATTAACGA  >  CP009273/289867‑290156
                                                                                                                                                 |                                                                                                                                                
gACAAGAGTAGAGCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGt                                                                                                                                                <  1:343833/148‑1 (MQ=255)
  cAAGAGTAGAGCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTaa                                                                                                                                              <  1:1068543/148‑1 (MQ=255)
      agTAGAGCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAAtttt                                                                                                                                          >  2:3132452/1‑148 (MQ=255)
         agagCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTccc                                                                                                                                       >  1:3133155/1‑148 (MQ=255)
         agagCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTccc                                                                                                                                       >  2:1050292/1‑148 (MQ=255)
          gagCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCt                                                                                                                                      >  1:494019/1‑148 (MQ=255)
          gagCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAAGGGCTTTTTCCATATTGTAATTTTCCCt                                                                                                                                      <  1:2188994/148‑1 (MQ=255)
             ccTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAg                                                                                                                                   <  1:2199104/148‑1 (MQ=255)
                  ccGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAgg                                                                                                                              <  2:137957/148‑1 (MQ=255)
                  ccGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAgg                                                                                                                              <  2:1938769/148‑1 (MQ=255)
                      cATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAg                                                                                                                          <  1:566832/148‑1 (MQ=255)
                        tCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGcc                                                                                                                        >  2:1293229/1‑148 (MQ=255)
                           aaGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCttt                                                                                                                     >  1:3237296/1‑148 (MQ=255)
                              taataaCCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCAc                                                                                                                  >  2:1069335/1‑148 (MQ=255)
                                ataaCCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGt                                                                                                                <  2:1863356/148‑1 (MQ=255)
                                ataaCCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGt                                                                                                                >  1:1070582/1‑148 (MQ=255)
                                 taaCCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCTCGta                                                                                                               >  1:3319420/1‑148 (MQ=255)
                                 taaCCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGta                                                                                                               <  2:3455366/148‑1 (MQ=255)
                                  aaCCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGtat                                                                                                              >  1:1244770/1‑148 (MQ=255)
                                   aCCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGtata                                                                                                             >  1:606255/1‑148 (MQ=255)
                                    cccTCAAAGGTTGTTCTATGAATAATAATAATGTTAGAATATCCTGAGACACTTTGATGCTTTACATCACAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATaa                                                                                                            <  1:2693216/148‑1 (MQ=255)
                                    cccTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATaa                                                                                                            >  2:3290618/1‑148 (MQ=255)
                                        cAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGt                                                                                                        <  2:1690283/148‑1 (MQ=255)
                                        cAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGt                                                                                                        <  2:903356/148‑1 (MQ=255)
                                        cAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGt                                                                                                        <  1:2476748/148‑1 (MQ=255)
                                         aaaGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATTCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGtt                                                                                                       <  2:1584133/148‑1 (MQ=255)
                                          aaGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGttt                                                                                                      <  2:735732/148‑1 (MQ=255)
                                            ggTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGttttt                                                                                                    <  1:3270760/148‑1 (MQ=255)
                                            ggTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGttttt                                                                                                    <  2:1190263/148‑1 (MQ=255)
                                            ggTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATCGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGttttt                                                                                                    >  2:357932/1‑148 (MQ=255)
                                               tgttCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCa                                                                                                 <  2:3498198/148‑1 (MQ=255)
                                               tgttCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATCGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCa                                                                                                 <  2:503358/148‑1 (MQ=255)
                                                  tCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCt                                                                                              <  2:324415/148‑1 (MQ=255)
                                                  tCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCATTATCGAGGATAGCCTATCACGTATAACTGTTTTTTCAtat                                                                                              >  2:2230352/1‑146 (MQ=255)
                                                   cTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCtt                                                                                             <  1:2486873/148‑1 (MQ=255)
                                                    tATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCttt                                                                                            >  2:583205/1‑148 (MQ=255)
                                                    tATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCttt                                                                                            <  1:2391906/148‑1 (MQ=255)
                                                       gaataagaataaTGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTccc                                                                                         >  2:1441207/1‑148 (MQ=255)
                                                        aataagaataaTGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTcccc                                                                                        >  1:1326027/1‑148 (MQ=255)
                                                        aataagaataaTGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTcccc                                                                                        >  1:2260403/1‑148 (MQ=255)
                                                         ataagaataaTGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATTTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCa                                                                                       <  2:3422344/148‑1 (MQ=255)
                                                          taagaataaTGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCAt                                                                                      >  2:3258758/1‑148 (MQ=255)
                                                          taagaataaTGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCAt                                                                                      >  1:2956058/1‑148 (MQ=255)
                                                            agaataaTGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGa                                                                                    <  1:3385611/148‑1 (MQ=255)
                                                             gaataaTGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGaa                                                                                   >  1:3243764/1‑148 (MQ=255)
                                                              aataatGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAAt                                                                                  >  2:1622260/1‑148 (MQ=255)
                                                              aataatGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAAt                                                                                  >  2:695001/1‑148 (MQ=255)
                                                              aataatGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAAt                                                                                  >  2:696568/1‑148 (MQ=255)
                                                                 aatGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATtca                                                                               >  1:1528948/1‑148 (MQ=255)
                                                                     tttGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTcatcat                                                                           <  1:2343123/148‑1 (MQ=255)
                                                                       tGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTa                                                                         <  1:1293229/148‑1 (MQ=255)
                                                                        gAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAg                                                                        >  1:1613259/1‑148 (MQ=255)
                                                                          atatCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCg                                                                      <  2:635024/148‑1 (MQ=255)
                                                                          atatCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCg                                                                      >  1:1995944/1‑148 (MQ=255)
                                                                           tatCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCACATGAATTCATCATTAGCGt                                                                     >  2:2779283/1‑148 (MQ=255)
                                                                            atCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGtt                                                                    <  2:237530/148‑1 (MQ=255)
                                                                                tGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGt                                                                <  2:164362/148‑1 (MQ=255)
                                                                                  agaCACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTg                                                              <  2:771434/148‑1 (MQ=255)
                                                                                    acacTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGaa                                                            >  2:1043477/1‑148 (MQ=255)
                                                                                    acacTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGaa                                                            <  2:1329547/148‑1 (MQ=255)
                                                                                    acacTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGaa                                                            >  1:1775134/1‑148 (MQ=255)
                                                                                    acacTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGAGTTCAGTTGaa                                                            >  2:2813352/1‑148 (MQ=255)
                                                                                     cacTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGaaa                                                           >  1:291476/1‑148 (MQ=255)
                                                                                     cacTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGaaa                                                           >  1:2089910/1‑148 (MQ=255)
                                                                                        tttGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATCGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTa                                                        >  1:2618732/1‑148 (MQ=255)
                                                                                          tGATGCTTTACCTCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAAt                                                      <  2:1966026/148‑1 (MQ=255)
                                                                                           gATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATc                                                     <  2:2172808/148‑1 (MQ=255)
                                                                                             tGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTg                                                   >  1:849495/1‑148 (MQ=255)
                                                                                             tGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTg                                                   <  2:2183682/148‑1 (MQ=255)
                                                                                               cTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCa                                                 >  1:1887788/1‑148 (MQ=255)
                                                                                                 ttACATCAAAATTTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAAt                                               <  2:1764180/148‑1 (MQ=255)
                                                                                                   aCATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATgg                                             >  1:1283757/1‑148 (MQ=255)
                                                                                                    cATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATggg                                            <  1:1642898/148‑1 (MQ=255)
                                                                                                    cATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATggg                                            <  2:1497495/148‑1 (MQ=255)
                                                                                                        aaaaTGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGAtt                                        >  1:2461373/1‑148 (MQ=255)
                                                                                                          aaTGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGAtttt                                      >  2:1990608/1‑148 (MQ=255)
                                                                                                          aaTGTCTTTGATTGTAAACGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGAtttt                                      >  1:3531180/1‑148 (MQ=255)
                                                                                                           aTGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGAttttt                                     <  1:1945659/148‑1 (MQ=255)
                                                                                                            tGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTa                                    <  1:2734497/148‑1 (MQ=255)
                                                                                                             gTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTaa                                   >  1:1189497/1‑148 (MQ=255)
                                                                                                              tCTTTGATTGTAAATGTAAATGGCTTTTTCCATATTGTAACTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAAt                                  >  1:237031/1‑148 (MQ=255)
                                                                                                                  tGATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGAc                              >  1:2996399/1‑148 (MQ=255)
                                                                                                                   gATTGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACt                             >  1:3359252/1‑148 (MQ=255)
                                                                                                                     ttGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCa                           >  1:3163453/1‑148 (MQ=255)
                                                                                                                     ttGTAAATGTAAATGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCa                           <  1:2858007/148‑1 (MQ=255)
                                                                                                                       gtaaatgtaaatGGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCatg                         <  2:2287134/148‑1 (MQ=255)
                                                                                                                        taaatgtaaatgGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCatga                        >  1:1394023/1‑148 (MQ=255)
                                                                                                                           atgtaaatgGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGt                     <  2:2098200/148‑1 (MQ=255)
                                                                                                                                 atgGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGttt               <  1:773911/148‑1 (MQ=255)
                                                                                                                                  tgGCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGtttt              <  1:3032816/148‑1 (MQ=255)
                                                                                                                                   ggCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGTTTTc             <  2:2878509/148‑1 (MQ=255)
                                                                                                                                   ggCTTTTTCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGTTTTc             <  2:295723/148‑1 (MQ=255)
                                                                                                                                       ttttCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGTTTTCCAGa         <  1:2404900/148‑1 (MQ=255)
                                                                                                                                       ttttCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGTTTTCCAGa         <  1:532733/148‑1 (MQ=255)
                                                                                                                                        tttCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGTTTTCCAGAt        >  2:292531/1‑148 (MQ=255)
                                                                                                                                          tCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGTTTTCCAGATTa      <  1:1710241/148‑1 (MQ=255)
                                                                                                                                          tCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGTTTTCCAGATTa      <  2:1326027/148‑1 (MQ=255)
                                                                                                                                          tCCATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGTTTTCCAGATTa      <  1:745850/148‑1 (MQ=255)
                                                                                                                                           ccATATTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGTTTTCCAGATTaa     <  1:1166544/148‑1 (MQ=255)
                                                                                                                                              tatTGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGTTTTCCAGATTAACGa  <  1:2557279/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GACAAGAGTAGAGCCTGCCCGGCATCAAAGTAATAACCCTCAAAGGTTGTTCTATGAATAAGAATAATGTTTGAATATCCTGAGACACTTTGATGCTTTACATCAAAATGTCTTTGATTGTAAATGTAAATGGCTTTTTCCATATCGTAATTTTCCCTTAGCGAGGATAGCCTTTCACGTATAACTGTTTTTTCATCTTTCCCCATGAATTCATCATTAGCGTTCAGTTGAAAGTAATCTGCAATGGGGATTTTTAATGGACTCATGATGTTCGTTTTCCAGATTAACGA  >  CP009273/289867‑290156

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: