Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,271,128 A→C 92.7% T145P (ACC→CCC)  yhaV → toxin of the SohB(PrlF)‑YhaV toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,271,1280AC92.7% 762.9 / 22.1 205T145P (ACC→CCC) yhaVtoxin of the SohB(PrlF)‑YhaV toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base A (10/5);  new base C (113/77);  total (123/82)
Fisher's exact test for biased strand distribution p-value = 7.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 207 total aligned reads displayed.

CTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCAGCCTTTCGCAGTCCTTCCAGGCTAT  >  CP009273/3270983‑3271270
                                                                                                                                                 |                                                                                                                                              
 tttCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAAccc                                                                                                                                             <  1:545211/148‑1 (MQ=255)
  ttCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGTACATCCTCCTGCCGACTGGGAACCCCt                                                                                                                                            <  2:516834/148‑1 (MQ=255)
  ttCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCt                                                                                                                                            >  2:1751846/1‑148 (MQ=255)
    cGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCa                                                                                                                                          >  1:1271671/1‑148 (MQ=255)
      ttATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCAcc                                                                                                                                        <  2:105279/148‑1 (MQ=255)
       tataGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCAccc                                                                                                                                       >  2:990915/1‑148 (MQ=255)
       tataGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCAccc                                                                                                                                       >  2:3195145/1‑148 (MQ=255)
         taGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCga                                                                                                                                     >  1:877478/1‑148 (MQ=255)
          aGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCgag                                                                                                                                    >  2:952852/1‑148 (MQ=255)
          aGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCgag                                                                                                                                    >  2:2125432/1‑148 (MQ=255)
          aGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCgag                                                                                                                                    >  2:1406261/1‑148 (MQ=255)
           gTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCgaga                                                                                                                                   <  2:1213169/148‑1 (MQ=255)
            tGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGaa                                                                                                                                  <  2:2784037/148‑1 (MQ=255)
            tGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGaa                                                                                                                                  >  2:2495977/1‑148 (MQ=255)
            tGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGaa                                                                                                                                  >  1:1496449/1‑148 (MQ=255)
            tGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGaa                                                                                                                                  >  2:1800477/1‑148 (MQ=255)
             gAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGaaa                                                                                                                                 <  1:1836663/148‑1 (MQ=255)
             gAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGaaa                                                                                                                                 >  2:2148054/1‑148 (MQ=255)
             gAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGaaa                                                                                                                                 >  1:994247/1‑148 (MQ=255)
              aaaaaGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAAc                                                                                                                                <  1:2963514/148‑1 (MQ=255)
                 aaGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACaga                                                                                                                             >  1:2035963/1‑148 (MQ=255)
                  agagAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACagaa                                                                                                                            >  2:2024545/1‑148 (MQ=255)
                  agagAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACagaa                                                                                                                            >  1:2133201/1‑148 (MQ=255)
                    agaAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACagaaga                                                                                                                          <  2:3109811/148‑1 (MQ=255)
                      aaaGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGaaa                                                                                                                        >  1:2564539/1‑148 (MQ=255)
                      aaaGTCATCATTCTGGGAAGGATGAACCATGAAAACACTCTGCGCACCAACGGTAAAAAAACAGATCCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGaaa                                                                                                                        >  1:1324726/1‑148 (MQ=255)
                       aaGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAAc                                                                                                                       >  2:488105/1‑148 (MQ=255)
                          tcatcaTTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCa                                                                                                                    >  1:1126706/1‑148 (MQ=255)
                            atcatTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCAtt                                                                                                                  <  1:485026/148‑1 (MQ=255)
                            atcatTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCAtt                                                                                                                  >  2:2651187/1‑148 (MQ=255)
                             tcatTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTg                                                                                                                 >  2:720857/1‑148 (MQ=255)
                              catTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGa                                                                                                                >  2:3167842/1‑148 (MQ=255)
                                ttCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATg                                                                                                              >  2:2285565/1‑148 (MQ=255)
                                ttCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATg                                                                                                              <  2:1858415/148‑1 (MQ=255)
                                 tCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACAGGGAACCCCTCACCCGAGAAACAGAAGAAACCCAGTGATgg                                                                                                             >  2:3033163/1‑148 (MQ=255)
                                   tGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGgtg                                                                                                           <  1:193304/148‑1 (MQ=255)
                                     ggatggatGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGtt                                                                                                         >  2:1497377/1‑148 (MQ=255)
                                     ggatggatGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGtt                                                                                                         <  2:392717/148‑1 (MQ=255)
                                       atggatgAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTca                                                                                                       >  2:3236541/1‑148 (MQ=255)
                                         ggatgAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTcaca                                                                                                     <  2:2362772/148‑1 (MQ=255)
                                            tgAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATg                                                                                                  <  2:2856607/148‑1 (MQ=255)
                                             gAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGc                                                                                                 <  2:1732759/148‑1 (MQ=255)
                                                  aTGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCccgcc                                                                                            >  1:896503/1‑148 (MQ=255)
                                                    gAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGa                                                                                          >  1:2305196/1‑148 (MQ=255)
                                                      aaaCACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAAc                                                                                        >  2:492907/1‑148 (MQ=255)
                                                      aaaCACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAAc                                                                                        >  1:2679336/1‑148 (MQ=255)
                                                      aaaCACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAAc                                                                                        >  1:2679243/1‑148 (MQ=255)
                                                      aaaCACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAAc                                                                                        >  2:3134207/1‑148 (MQ=255)
                                                        acacTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGc                                                                                      <  1:2193735/148‑1 (MQ=255)
                                                         cacTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCa                                                                                     >  2:2407895/1‑148 (MQ=255)
                                                          acTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCAt                                                                                    <  1:3050470/148‑1 (MQ=255)
                                                                  cACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTg                                                                            <  2:567856/148‑1 (MQ=255)
                                                                   aCCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGa                                                                           >  1:1459214/1‑148 (MQ=255)
                                                                     cTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGAcc                                                                         >  1:650928/1‑148 (MQ=255)
                                                                        cGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGc                                                                      >  1:2765239/1‑148 (MQ=255)
                                                                          gTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCgg                                                                    >  1:1234221/1‑148 (MQ=255)
                                                                          gTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCgg                                                                    >  1:2169716/1‑148 (MQ=255)
                                                                           tAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCggg                                                                   >  1:1707116/1‑148 (MQ=255)
                                                                            aaaaaaaCAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGc                                                                  >  2:1298098/1‑148 (MQ=255)
                                                                            aaaaaaaCAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGc                                                                  >  2:562548/1‑148 (MQ=255)
                                                                             aaaaaaCAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCa                                                                 >  1:2364611/1‑148 (MQ=255)
                                                                             aaaaaaCAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCa                                                                 >  1:113467/1‑148 (MQ=255)
                                                                              aaaaaCAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCAt                                                                >  1:2305895/1‑148 (MQ=255)
                                                                                aaaCAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATcc                                                              <  2:2675192/148‑1 (MQ=255)
                                                                                 aaCAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCt                                                             <  2:930823/148‑1 (MQ=255)
                                                                                   cAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCttt                                                           <  2:2563925/148‑1 (MQ=255)
                                                                                     gATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTtct                                                         >  1:477591/1‑148 (MQ=255)
                                                                                     gATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTtct                                                         >  1:2841246/1‑148 (MQ=255)
                                                                                       tGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCa                                                       >  2:635060/1‑148 (MQ=255)
                                                                                       tGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCACGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCa                                                       >  1:1249067/1‑148 (MQ=255)
                                                                                          cTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTc                                                    <  1:1751846/148‑1 (MQ=255)
                                                                                           tataCCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTcc                                                   <  1:1095453/148‑1 (MQ=255)
                                                                                           tataCCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTcc                                                   >  1:873539/1‑148 (MQ=255)
                                                                                            ataCCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTccc                                                  <  2:1237354/148‑1 (MQ=255)
                                                                                              aCCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCg                                                >  1:1431584/1‑148 (MQ=255)
                                                                                                  tATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCa                                            >  2:2608959/1‑148 (MQ=255)
                                                                                                   aTTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAg                                           >  1:1963989/1‑148 (MQ=255)
                                                                                                   aTTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAg                                           >  1:1781982/1‑148 (MQ=255)
                                                                                                       aGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATc                                       <  2:468594/148‑1 (MQ=255)
                                                                                                        gCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCa                                      >  1:892494/1‑148 (MQ=255)
                                                                                                        gCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCa                                      >  1:210896/1‑148 (MQ=255)
                                                                                                           aaaTGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCTCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACtt                                   >  1:427177/1‑148 (MQ=255)
                                                                                                              tGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCaa                                >  1:696398/1‑148 (MQ=255)
                                                                                                               gTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCCTCCTTTCTCACTCCCCGACCAGAATCACTTCAAc                               >  2:2929574/1‑148 (MQ=255)
                                                                                                               gTTAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAAc                               >  2:2743084/1‑148 (MQ=255)
                                                                                                                ttAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAAcc                              >  1:2532841/1‑148 (MQ=255)
                                                                                                                  aaaaaGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTTGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAAcccc                            <  2:1271671/148‑1 (MQ=255)
                                                                                                                    aaaGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCAg                          >  2:1075561/1‑148 (MQ=255)
                                                                                                                    aaaGAGGACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCAg                          >  1:1898010/1‑148 (MQ=255)
                                                                                                                         ggACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCAGCCttt                     <  2:1126706/148‑1 (MQ=255)
                                                                                                                         ggACATCCTCCTGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCAGCCttt                     <  1:3176198/148‑1 (MQ=255)
                                                                                                                                ctcctGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCCTTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCAGCCTTTCGCAGTc              <  1:173728/148‑1 (MQ=255)
                                                                                                                                 tcctGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCAGCCTTTCGCAGTcc             >  2:586913/1‑148 (MQ=255)
                                                                                                                                  cctGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCAGCCTTTCGCAGTCCt            >  2:2460688/1‑148 (MQ=255)
                                                                                                                                    tGCCGACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCAGCCTTTCGCAGTCCTTc          <  1:3195145/148‑1 (MQ=255)
                                                                                                                                        gACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCAGCCTTTCGCAGTCCTTCCAgg      >  2:3016393/1‑148 (MQ=255)
                                                                                                                                        gACTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCAGCCTTTCGCAGTCCTTCCAgg      >  2:3016235/1‑148 (MQ=255)
                                                                                                                                          cTGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCAGCCTTTCGCAGTCCTTCCAGGCt    <  1:633051/148‑1 (MQ=255)
                                                                                                                                           tGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCAGCCTTTCGCAGTCCTTCCAGGCTa   <  1:941275/148‑1 (MQ=255)
                                                                                                                                           tGGGAACCCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCAGCCTTTCGCAGTCCTTCCAGGCTa   <  1:2595695/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
CTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATTCGTGTTGACCAGCGGGCATCCTTTCTCACTCCCCGACCAGAATCACTTCAACCCCAGCCTTTCGCAGTCCTTCCAGGCTAT  >  CP009273/3270983‑3271270

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: