Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,407,969 C→A 100% N249K (AAC→AAA acrE → cytoplasmic membrane lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,407,9690CA97.7% 803.9 / 8.4 263N249K (AAC→AAAacrEcytoplasmic membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (4/2);  new base A (145/111);  total (150/113)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 263 total aligned reads displayed.

AACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTGTGCGTGCACGGATTGA  >  CP009273/3407824‑3408112
                                                                                                                                                 |                                                                                                                                               
     gCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCaaa                                                                                                                                          <  2:2976630/148‑1 (MQ=255)
        aCTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCt                                                                                                                                       >  2:2729048/1‑148 (MQ=255)
        aCTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCt                                                                                                                                       >  2:894640/1‑148 (MQ=255)
          tGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTAt                                                                                                                                     <  1:2908825/148‑1 (MQ=255)
           gTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATc                                                                                                                                    >  1:2732266/1‑148 (MQ=255)
             ccAGCAGCTCTATCCTATCTACGTTTATGTCCCCCATACCAGAAACTATTTTATGTGGCTGACGCAATCGGTAGAGCAAGTAATTTTGAAAAAGTAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATccc                                                                                                                                  <  1:634429/148‑1 (MQ=255)
             ccAGCAGCTCGTCCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCATGGAAATTGGAATAAGCAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATccc                                                                                                                                  <  1:3376116/148‑1 (MQ=255)
              cagcagCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATcccc                                                                                                                                 <  2:1199280/148‑1 (MQ=255)
                gcagcTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACTATTTTATGAGGCTGAAGCAATCCGTATAGCAAGGAAATTTGCATAAGTAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTg                                                                                                                               <  1:1349967/148‑1 (MQ=255)
                 cagcTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGa                                                                                                                              >  1:233404/1‑148 (MQ=255)
                   gcTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGaaa                                                                                                                            >  1:312155/1‑148 (MQ=255)
                    cTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAg                                                                                                                           >  1:1884845/1‑148 (MQ=255)
                    cTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAg                                                                                                                           >  2:235369/1‑148 (MQ=255)
                      cGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGt                                                                                                                         <  2:1846465/148‑1 (MQ=255)
                      cGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGt                                                                                                                         <  1:2295504/148‑1 (MQ=255)
                       gATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTa                                                                                                                        <  2:3434319/148‑1 (MQ=255)
                             aTCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGc                                                                                                                  >  1:440034/1‑148 (MQ=255)
                             aTCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGc                                                                                                                  >  2:2673088/1‑148 (MQ=255)
                                 aCGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAAtt                                                                                                              <  1:1702513/148‑1 (MQ=255)
                                    ttgttgTGACCCAATCCCGAAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGTAAATTTGAATAAGTAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATtctc                                                                                                           <  1:604119/148‑1 (MQ=255)
                                     tGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTcc                                                                                                          <  1:1860541/148‑1 (MQ=255)
                                       aTGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGa                                                                                                        <  2:3564158/148‑1 (MQ=255)
                                       aTGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGa                                                                                                        >  1:1721425/1‑148 (MQ=255)
                                         gtgACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGTAATTTGCATAACGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGAtg                                                                                                      <  1:1268058/148‑1 (MQ=255)
                                          tgACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGAtgt                                                                                                     >  2:202643/1‑148 (MQ=255)
                                          tgACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGAtgt                                                                                                     >  2:3152067/1‑148 (MQ=255)
                                            aCCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGa                                                                                                   >  2:737244/1‑148 (MQ=255)
                                              ccAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTCCAATTCTCCGATGTGAcc                                                                                                 >  2:1933392/1‑148 (MQ=255)
                                                    cAGCACCGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTtgat                                                                                           <  1:3014362/148‑1 (MQ=255)
                                                    cAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTtgat                                                                                           <  2:3237332/148‑1 (MQ=255)
                                                      gCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTtgatga                                                                                         >  1:3058727/1‑148 (MQ=255)
                                                       cAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGaa                                                                                        >  1:1052671/1‑148 (MQ=255)
                                                         aCGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCACCGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAg                                                                                      >  2:1066846/1‑148 (MQ=255)
                                                         aCGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAg                                                                                      >  1:2592284/1‑148 (MQ=255)
                                                         aCGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAg                                                                                      >  1:805339/1‑148 (MQ=255)
                                                         aCGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAg                                                                                      >  1:1102274/1‑148 (MQ=255)
                                                          cGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGc                                                                                     <  1:737244/148‑1 (MQ=255)
                                                                    aGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCa                                                                           >  1:1245418/1‑148 (MQ=255)
                                                                      gCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATa                                                                         >  2:2446932/1‑148 (MQ=255)
                                                                           aGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCt                                                                    <  2:1348583/148‑1 (MQ=255)
                                                                            gCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTa                                                                   >  1:2618727/1‑148 (MQ=255)
                                                                            gCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTa                                                                   >  2:1735356/1‑148 (MQ=255)
                                                                             cAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTAc                                                                  >  2:777494/1‑148 (MQ=255)
                                                                              aaTCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACg                                                                 <  2:901480/148‑1 (MQ=255)
                                                                              aaTCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACg                                                                 >  2:1419909/1‑148 (MQ=255)
                                                                               aTCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGt                                                                >  1:2446913/1‑148 (MQ=255)
                                                                                tCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTg                                                               >  1:585160/1‑148 (MQ=255)
                                                                                tCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTg                                                               <  1:1753364/148‑1 (MQ=255)
                                                                                  cGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCt                                                             >  2:3375044/1‑148 (MQ=255)
                                                                                   gTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTg                                                            <  2:1884818/148‑1 (MQ=255)
                                                                                    tAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGt                                                           >  1:1022051/1‑148 (MQ=255)
                                                                                    tAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGt                                                           <  2:876125/148‑1 (MQ=255)
                                                                                    tAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGt                                                           >  1:2382440/1‑148 (MQ=255)
                                                                                       agCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCtt                                                        <  1:148115/148‑1 (MQ=255)
                                                                                       agCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCtt                                                        <  2:2824217/148‑1 (MQ=255)
                                                                                       agCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCtt                                                        >  1:3107412/1‑148 (MQ=255)
                                                                                       agCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCtt                                                        >  2:148115/1‑148 (MQ=255)
                                                                                         cAAGGAAATTTGAATAAGGAACACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTcc                                                      <  1:954245/148‑1 (MQ=255)
                                                                                                aTTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAAccc                                               <  1:387869/148‑1 (MQ=255)
                                                                                                 tttGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCg                                              >  1:2261127/1‑148 (MQ=255)
                                                                                                  ttGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGc                                             >  2:1804935/1‑148 (MQ=255)
                                                                                                     cATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAAc                                          >  1:2021428/1‑148 (MQ=255)
                                                                                                      aTAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACa                                         <  2:1587499/148‑1 (MQ=255)
                                                                                                       tAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACAt                                        <  1:2180676/148‑1 (MQ=255)
                                                                                                        aaGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATa                                       >  1:3212053/1‑148 (MQ=255)
                                                                                                        aaGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATa                                       <  1:1155402/148‑1 (MQ=255)
                                                                                                        aaGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATa                                       <  2:2048827/148‑1 (MQ=255)
                                                                                                         aGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATAc                                      >  1:987221/1‑148 (MQ=255)
                                                                                                           gAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGc                                    <  2:2097412/148‑1 (MQ=255)
                                                                                                             ataCGACACCAGAATCGTAGAGTTTGTCATTGATAACGGTCAAACCTTTCCCGTGAAAGGTTCGGTGCAATTCTCCGATGTGACGGTTGATGAAAGCACGGGCTCCATAACACTACGTGCTGTCTTCCATAACCCGCTTCATAAGCtt                                  <  2:2895817/146‑1 (MQ=255)
                                                                                                              aaCGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAATGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCttt                                 <  2:197351/148‑1 (MQ=255)
                                                                                                              aaCGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAATGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCttt                                 <  2:233398/148‑1 (MQ=255)
                                                                                                              aaCGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGc                                    <  1:2090959/145‑1 (MQ=255)
                                                                                                              aaCGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCttt                                 <  2:1258000/148‑1 (MQ=255)
                                                                                                                cGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTg                               <  2:710559/148‑1 (MQ=255)
                                                                                                                cGCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTg                               >  2:1338251/1‑148 (MQ=255)
                                                                                                                 gCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGc                              <  1:70824/148‑1 (MQ=255)
                                                                                                                 gCCACCAGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGc                              <  1:2805064/148‑1 (MQ=255)
                                                                                                                     ccaGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCggg                          <  2:436192/148‑1 (MQ=255)
                                                                                                                     ccaGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCggg                          >  1:2749680/1‑148 (MQ=255)
                                                                                                                     ccaGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCggg                          <  2:2951277/148‑1 (MQ=255)
                                                                                                                     ccaGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCggg                          >  1:551123/1‑148 (MQ=255)
                                                                                                                     ccaGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCACCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCggg                          >  2:590569/1‑148 (MQ=255)
                                                                                                                      caGCAACGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGt                         >  1:870791/1‑148 (MQ=255)
                                                                                                                           aCGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGtt                    >  2:1109082/1‑148 (MQ=255)
                                                                                                                           aCGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGtt                    >  1:1924535/1‑148 (MQ=255)
                                                                                                                           aCGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGtt                    >  1:788867/1‑148 (MQ=255)
                                                                                                                            cGTAGAGTTGGTCATGGAAAAAGGTTAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGttt                   >  2:2527535/1‑148 (MQ=255)
                                                                                                                            cGTAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGttt                   >  2:621933/1‑148 (MQ=255)
                                                                                                                              tAGAGTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTtgt                 >  2:211340/1‑148 (MQ=255)
                                                                                                                               agagTTGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTtgtg                <  1:2885121/148‑1 (MQ=255)
                                                                                                                                   ttGGTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTgtgcgtg            >  2:3172702/1‑148 (MQ=255)
                                                                                                                                     ggTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTGTGCGTGCa          <  2:805339/148‑1 (MQ=255)
                                                                                                                                     ggTCATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTGTGCGTGCa          <  2:2592284/148‑1 (MQ=255)
                                                                                                                                        cATGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTGTGCGTGCACgg       >  1:2063791/1‑148 (MQ=255)
                                                                                                                                         aTGGAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTGTGCGTGCACGGa      <  1:891199/148‑1 (MQ=255)
                                                                                                                                           ggAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTGTGCGTGCACGGAtt    >  1:1645589/1‑148 (MQ=255)
                                                                                                                                           ggAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTGTGCGTGCACGGAtt    >  1:3428384/1‑148 (MQ=255)
                                                                                                                                            gAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTGTGCGTGCACGGATtg   >  2:426068/1‑148 (MQ=255)
                                                                                                                                            gAAAAAGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTGTGCGTGCACGGATtg   >  2:1814507/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
AACTGGCGACTGTCCAGCAGCTCGATCCTATCTACGTTGATGTGACCCAATCCAGCAACGATTTTATGAGGCTGAAGCAATCCGTAGAGCAAGGAAATTTGCATAAGGAAAACGCCACCAGCAACGTAGAGTTGGTCATGGAAAACGGTCAAACCTATCCCCTGAAAGGTACGCTGCAATTCTCCGATGTGACCGTTGATGAAAGCACCGGCTCCATAACCCTACGTGCTGTCTTCCCTAACCCGCAACATACGCTTTTGCCGGGTATGTTTGTGCGTGCACGGATTGA  >  CP009273/3407824‑3408112

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: