Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 565,475 C→T 25.5% A373V (GCT→GTT)  ybcK → DLP12 prophage; putative recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273565,4750CT25.5% 180.0 / 67.7 98A373V (GCT→GTT) ybcKDLP12 prophage; putative recombinase
Reads supporting (aligned to +/- strand):  ref base C (31/42);  new base T (11/14);  total (42/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

GTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGTTACCGCTATAGCAGAGAAAATAAGACTATTAGATAAGGAATTACGAAGGGCTTCGGTATCAT  >  CP009273/565328‑565621
                                                                                                                                                   |                                                                                                                                                  
gTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTa                                                                                                                                                  >  1:819568/1‑150 (MQ=255)
gTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTa                                                                                                                                                  >  1:3248634/1‑150 (MQ=255)
 tttGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTaa                                                                                                                                                 >  1:4441483/1‑150 (MQ=255)
 tttGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTaa                                                                                                                                                 <  2:2676643/150‑1 (MQ=255)
 tttGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTaa                                                                                                                                                 >  2:1204738/1‑150 (MQ=255)
 tttGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTaa                                                                                                                                                 <  2:2676651/150‑1 (MQ=255)
  ttGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGTTAAt                                                                                                                                                <  1:41456/150‑1 (MQ=255)
  ttGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGTTAAt                                                                                                                                                <  1:42861/150‑1 (MQ=255)
     cATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGaa                                                                                                                                             >  2:3039214/1‑150 (MQ=255)
     cATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGaa                                                                                                                                             <  1:1375280/150‑1 (MQ=255)
     cATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGaa                                                                                                                                             >  2:3039544/1‑150 (MQ=255)
      aTGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGaaa                                                                                                                                            <  2:1871671/150‑1 (MQ=255)
      aTGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGaaa                                                                                                                                            <  2:3652381/150‑1 (MQ=255)
          cTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACttt                                                                                                                                        >  2:2433216/1‑150 (MQ=255)
          cTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACttt                                                                                                                                        <  1:4441528/150‑1 (MQ=255)
            attatGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGTTAATGAAACTTTag                                                                                                                                      >  2:2773694/1‑150 (MQ=255)
            attatGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGTTAATGAAACTTTag                                                                                                                                      >  2:675670/1‑150 (MQ=255)
            attatGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTat                                                                                                                                      >  2:3233127/1‑149 (MQ=255)
            attatGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTag                                                                                                                                      <  1:228680/150‑1 (MQ=255)
            attatGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTag                                                                                                                                      <  1:1613088/150‑1 (MQ=255)
            attatGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTag                                                                                                                                      >  2:3485279/1‑150 (MQ=255)
            attatGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTag                                                                                                                                      >  2:3485297/1‑150 (MQ=255)
             ttatGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTaga                                                                                                                                     <  2:1611495/150‑1 (MQ=255)
             ttatGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTaga                                                                                                                                     <  1:1502935/150‑1 (MQ=255)
             ttatGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTaga                                                                                                                                     <  1:1502923/150‑1 (MQ=255)
              tatGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTagag                                                                                                                                    >  1:2885063/1‑150 (MQ=255)
               atGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTAGAGt                                                                                                                                   <  1:1615820/150‑1 (MQ=255)
                tGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGTTAATGAAACTTTAGAGtt                                                                                                                                  <  2:1908695/150‑1 (MQ=255)
                tGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTAAGAGtt                                                                                                                                  >  2:1135291/1‑150 (MQ=255)
                     gTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGTTAATGAAACTTTAGAGTTaaaaa                                                                                                                             >  1:2179416/1‑150 (MQ=255)
                      tCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGTTAATGAAACTTTAGAGTTAAAAat                                                                                                                            <  2:4344861/150‑1 (MQ=255)
                      tCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAat                                                                                                                            <  2:2712404/150‑1 (MQ=255)
                      tCCGATGAGAAGAGTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGTTAATGAAACTTTAGAGTTAAAAat                                                                                                                            <  2:4344112/150‑1 (MQ=255)
                             agaagaTTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGTTAATGAAACTTTAGAGTTAAAAATTATTGAg                                                                                                                     <  1:1062348/150‑1 (MQ=255)
                             agaagaTTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGTTAATGAAACTTTAGAGTTAAAAATTATTGAg                                                                                                                     <  1:1063038/150‑1 (MQ=255)
                                agaTTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGTTAATGAAACTTTAGAGTTAAAAATTATTGAGCtt                                                                                                                  <  2:2376210/150‑1 (MQ=255)
                                agaTTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGTTAATGAAACTTTAGAGTTAAAAATTATTGAGCtt                                                                                                                  <  2:4700076/150‑1 (MQ=255)
                                  aTTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCa                                                                                                                <  1:4128277/150‑1 (MQ=255)
                                  aTTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCa                                                                                                                <  2:2311816/150‑1 (MQ=255)
                                  aTTACATCGATGTGACAGGCCATAAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCa                                                                                                                <  2:2852208/150‑1 (MQ=255)
                                      cATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATg                                                                                                            <  1:388198/150‑1 (MQ=255)
                                      cATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATg                                                                                                            <  1:387421/150‑1 (MQ=255)
                                            tgtgACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTGTTGAGCTTCAGATGAAAAtt                                                                                                      <  1:1300286/150‑1 (MQ=255)
                                            tgtgACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAAtt                                                                                                      <  1:4682100/150‑1 (MQ=255)
                                            tgtgACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAAtt                                                                                                      <  1:1299945/150‑1 (MQ=255)
                                              tgACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATtaa                                                                                                    >  1:1053364/1‑150 (MQ=255)
                                                aCAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATtaata                                                                                                  >  2:2274806/1‑150 (MQ=255)
                                                          aaTAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGt                                                                                        >  2:1226410/1‑150 (MQ=255)
                                                             aaaaaGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGTTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGc                                                                                     >  2:2087317/1‑150 (MQ=255)
                                                               aaaGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCAt                                                                                   <  2:2896208/150‑1 (MQ=255)
                                                                aaGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCAtt                                                                                  <  2:913739/150‑1 (MQ=255)
                                                                    gaTTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGTTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATtgtct                                                                              <  1:1314873/150‑1 (MQ=255)
                                                                    gaTTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGTTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATtgtct                                                                              <  1:1314736/150‑1 (MQ=255)
                                                                    gaTTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGTTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATtgtct                                                                              <  1:1314696/150‑1 (MQ=255)
                                                                          gTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGTTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTcgcg                                                                        >  2:3603456/1‑150 (MQ=255)
                                                                              aTTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTg                                                                    <  1:138698/150‑1 (MQ=255)
                                                                                tataATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGaa                                                                  >  1:3777885/1‑150 (MQ=255)
                                                                                tataATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGaa                                                                  <  1:72199/150‑1 (MQ=255)
                                                                                  taataTCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGTTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGt                                                                >  2:662043/1‑150 (MQ=255)
                                                                                  taataTCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGTTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGt                                                                >  2:662224/1‑150 (MQ=255)
                                                                                  taataTCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGTTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGt                                                                >  2:662105/1‑150 (MQ=255)
                                                                                       tCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGTTACCg                                                           <  1:4088109/150‑1 (MQ=255)
                                                                                        cATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGTTACCGc                                                          <  2:1510201/150‑1 (MQ=255)
                                                                                          ttAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGTTACCGCta                                                        >  2:729142/1‑150 (MQ=255)
                                                                                           tAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGTTACCGCtat                                                       >  2:2067045/1‑150 (MQ=255)
                                                                                           tAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGTTACCGCtat                                                       >  2:403806/1‑150 (MQ=255)
                                                                                           tAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGTTACCGCtat                                                       <  2:3880167/150‑1 (MQ=255)
                                                                                               gacttGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGTTACCGCTATAGCa                                                   <  2:4344150/147‑1 (MQ=255)
                                                                                                 aTTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGTTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGTTACCGCTATAGCaga                                                 <  1:3603456/150‑1 (MQ=255)
                                                                                                     cTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGTTACCGCTATAGCAGAGaaa                                             <  1:2647350/150‑1 (MQ=255)
                                                                                                       ttttAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGTTACCGCTATAGCAGAGAAAAt                                           >  1:1939398/1‑150 (MQ=255)
                                                                                                        tttAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGTTACCGCTATAGCAGAGAAAATa                                          <  1:39790/150‑1 (MQ=255)
                                                                                                        tttAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGTTACCGCTATAGCAGAGAAAATa                                          <  1:40887/150‑1 (MQ=255)
                                                                                                         ttAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGTTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGTTACCGCTATAGCAGAGAAAATaa                                         <  1:4579423/150‑1 (MQ=255)
                                                                                                         ttAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGTTACCGCTATAGCAGAGAAAATaa                                         <  2:1391500/150‑1 (MQ=255)
                                                                                                         ttAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGTTACCGCTATAGCAGAGAAAATaa                                         <  1:4748330/150‑1 (MQ=255)
                                                                                                                        ttCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTAAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGTTACCGCTATAGCAGAGAAAATAAGACTATTAGATAAgg                          <  1:1135291/150‑1 (MQ=255)
                                                                                                                                 atcaaaaCAAGAAATATTCTAATTAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGTTACCGCTATAGCAGAGAAAATAAGACTATTAGATAAGGAATTACGaa                 <  2:4165108/147‑1 (MQ=255)
                                                                                                                                  tGAAAACAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGTTACCGCTATAGCAGAGAAAATAAGACTATTAGATAAGGAATTACGAAg                >  2:632917/1‑150 (MQ=255)
                                                                                                                                   gAAAACAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGTTACCGCTATAGCAGAGAAAATAAGACTATTAGATAAGGAATTACGAAgt               >  2:2579359/1‑149 (MQ=255)
                                                                                                                                     aaaCAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGTTACCGCTATAGCAGAGAAAATAAGACTATTAGATAAGGAATTACGAAGGGc             >  2:3609775/1‑150 (MQ=255)
                                                                                                                                     aaaCAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGTTACCGCTATAGCAGAGAAAATAAGACTATTAGATAAGGAATTACGAAGGGc             <  2:517711/150‑1 (MQ=255)
                                                                                                                                     aaaCAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGTTACCGCTATAGCAGAGAAAATAAGACTATTAGATAAGGAATTACGAAGGGc             >  2:3149278/1‑150 (MQ=255)
                                                                                                                                     aaaCAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGTTACCGCTATAGCAGAGAAAATAAGACTATTAGATAAGGAATTACGAAGGGc             <  2:517693/150‑1 (MQ=255)
                                                                                                                                       aCAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGTTACCGCTATAGCAGAGAAAATAAGACTATTAGATAAGGAATTACGAAGGGCtt           >  1:2200156/1‑150 (MQ=255)
                                                                                                                                       aCAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGTTACCGCTATAGCAGAGAAAATAAGACTATTAGATAAGGAATTACGAAGGGCtt           >  1:2200163/1‑150 (MQ=255)
                                                                                                                                       aCAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGTTACCGCTATAGCAGAGAAAATAAGACTATTAGATAAGGAATTACGAAGGGCtt           >  1:3098909/1‑150 (MQ=255)
                                                                                                                                        cAAGAAAGATGTTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGTTACCGCTATAGCAGAGAAAATAAGACTATTAGATAAGGAATTACGAAGGGCTTc          >  2:4561315/1‑150 (MQ=255)
                                                                                                                                        cAAGAAAGATGTTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGTTACCGCTATAGCAGAGAAAATAAGACTATTAGATAAGGAATTACGAAGGGCTTc          >  2:4561927/1‑150 (MQ=255)
                                                                                                                                        cAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGTTACCGCTATAGCAGAGAAAATAAGACTATTAGATAAGGAATTACGAAGGGCTTc          <  1:1069721/150‑1 (MQ=255)
                                                                                                                                         aagaaagaTGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGAGCCTGAAGTTACCACTATAGCAGATAAAATAATACTATTAGATAAGGAATTACTAAGGTCTTCg         <  2:4318691/150‑1 (MQ=255)
                                                                                                                                          agaaagaTGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGTTACCGCTATAGCAGAGAAAATAAGACTATTAGATAAGGAATTACGAAGGGCTTCgg        >  2:4454930/1‑150 (MQ=255)
                                                                                                                                           gaaagaTGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGTTACCGCTATAGCAGAGAAAATAAGACTATTAGATAAGGAATTACGAAGGGCTTCGGt       <  1:1283478/150‑1 (MQ=255)
                                                                                                                                           gaaagaTGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGTTACCGCTATAGCAGAGAAAATAAGACTATTAGATAAGGAATTACGAAGGGCTTCGGt       <  1:3748760/150‑1 (MQ=255)
                                                                                                                                             aagaTGTTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGTTACCGCTATAGCAGAGAAAATAAGACTATTAGATAAGGAATTACGAAGGGCTTCGGTAt     >  2:2726602/1‑150 (MQ=255)
                                                                                                                                             aagaTGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGTTACCGCTATAGCAGAGAAAATAAGACTAGTAGATAAGGAATTACGAAGGGCTTCGGTAt     >  1:1164052/1‑150 (MQ=255)
                                                                                                                                                aTGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGTTACCGCTATAGCAGAGAAAATAAGACTATTAGATAAGGAATTACGAAGGGCTTCGGTATCAt  >  2:2300580/1‑150 (MQ=255)
                                                                                                                                                aTGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGTTACCGCTATAGCAGAGAAAATAAGACTATTAGATAAGGAATTACGAAGGGCTTCGGTATCAt  >  2:2299918/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAATTCAACCAGTTGAAAACAAGAAAGATGCTAATGAAACTTTAGAGTTAAAAATTATTGAGCTTCAGATGAAAATTAATAATTTAATCGTTGCATTGTCTGTCGCGCCTGAAGTTACCGCTATAGCAGAGAAAATAAGACTATTAGATAAGGAATTACGAAGGGCTTCGGTATCAT  >  CP009273/565328‑565621

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: