Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,800,129 C→T 90.8% W71* (TGG→TAG)  ydiY ← acid‑inducible putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,800,1290CT90.8% 671.0 / 48.8 218W71* (TGG→TAG) ydiYacid‑inducible putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (6/14);  new base T (107/91);  total (113/105)
Fisher's exact test for biased strand distribution p-value = 5.86e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

Only 100 of 218 total aligned reads displayed.

GCTTGCCTGACCAAATAAATAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTG  >  CP009273/1799984‑1800273
                                                                                                                                                 |                                                                                                                                                
gcTTGCCTGACCAAATAAATAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAg                                                                                                                                                <  2:34091/148‑1 (MQ=255)
 cTTGCCTGACCAAATAAATAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAgg                                                                                                                                               >  1:2717771/1‑148 (MQ=255)
 cTTGCCTGACCAAATAAATAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAgg                                                                                                                                               >  2:752620/1‑148 (MQ=255)
      cTGACCAAATAAATAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCAtg                                                                                                                                          >  2:2481473/1‑148 (MQ=255)
          ccAAATAAATAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCAtggtgg                                                                                                                                      <  2:2361183/148‑1 (MQ=255)
          ccAAATAAATAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCAtggtgg                                                                                                                                      <  1:2113501/148‑1 (MQ=255)
           cAAATAAATAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATggtggt                                                                                                                                     >  2:3256669/1‑148 (MQ=255)
            aaataaatAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGgtggtg                                                                                                                                    >  2:1818960/1‑148 (MQ=255)
              ataaataGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTc                                                                                                                                  >  2:1570914/1‑148 (MQ=255)
               taaataGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCa                                                                                                                                 >  1:2384444/1‑148 (MQ=255)
                aaataGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAg                                                                                                                                >  1:2512880/1‑148 (MQ=255)
                     gTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGCTGTCAGCCGTa                                                                                                                           >  1:2892776/1‑148 (MQ=255)
                      tCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCCCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTaa                                                                                                                          >  2:1374181/1‑148 (MQ=255)
                        aTAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGt                                                                                                                        >  1:3215055/1‑148 (MQ=255)
                        aTAATCAGTTAAATAAACACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGt                                                                                                                        <  1:574649/148‑1 (MQ=255)
                         tAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATTGTGGTGTCAGCCGTAAGTg                                                                                                                       >  1:48478/1‑148 (MQ=255)
                         tAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATTGTGGTGTCAGCCGTAAGTg                                                                                                                       >  1:48570/1‑148 (MQ=255)
                            tCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTgagg                                                                                                                    <  1:1198998/148‑1 (MQ=255)
                              aGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTgaggag                                                                                                                  <  2:2675306/148‑1 (MQ=255)
                              aGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTgaggag                                                                                                                  <  2:1506090/148‑1 (MQ=255)
                                 taaattaaaACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCtt                                                                                                               <  1:3242647/148‑1 (MQ=255)
                                 taaattaaaACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCtt                                                                                                               <  2:640822/148‑1 (MQ=255)
                                 taaattaaaACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCtt                                                                                                               >  2:3196264/1‑148 (MQ=255)
                                   aattaaaACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCtttt                                                                                                             <  1:3256671/148‑1 (MQ=255)
                                      taacACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGt                                                                                                          <  2:2930457/148‑1 (MQ=255)
                                      taaaACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGt                                                                                                          >  1:1318308/1‑148 (MQ=255)
                                          aCGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTg                                                                                                      >  1:1335333/1‑148 (MQ=255)
                                             gCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCt                                                                                                   >  2:3017596/1‑148 (MQ=255)
                                              cTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTg                                                                                                  >  2:2707377/1‑148 (MQ=255)
                                              cTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTg                                                                                                  >  2:1300382/1‑148 (MQ=255)
                                               tACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGTAGCTTTTCGTATTGCCTGa                                                                                                 <  2:2869256/148‑1 (MQ=255)
                                               tACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGa                                                                                                 <  2:1827308/148‑1 (MQ=255)
                                                  gcccgcccgccGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATtg                                                                                              <  2:1434553/148‑1 (MQ=255)
                                                   cccgcccgccGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATtgt                                                                                             >  2:2991576/1‑148 (MQ=255)
                                                    ccgcccgccGCATATTTTTCTGAAGAACGCTCATAATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATtgtg                                                                                            >  2:1636208/1‑148 (MQ=255)
                                                        ccgccgCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAg                                                                                        <  2:1643575/148‑1 (MQ=255)
                                                        ccgccgCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAg                                                                                        <  1:63807/148‑1 (MQ=255)
                                                         cgccgcATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGa                                                                                       <  1:2558189/148‑1 (MQ=255)
                                                         cgccgcATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGa                                                                                       <  1:1409515/148‑1 (MQ=255)
                                                           ccgcATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGata                                                                                     >  2:1805232/1‑148 (MQ=255)
                                                           ccgcATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGata                                                                                     >  2:2621823/1‑148 (MQ=255)
                                                            cgcATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGatat                                                                                    <  1:1635520/148‑1 (MQ=255)
                                                            cgcATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGatat                                                                                    <  1:2489357/148‑1 (MQ=255)
                                                             gcATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATc                                                                                   >  1:1673630/1‑148 (MQ=255)
                                                               atatTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCt                                                                                 >  1:1129514/1‑148 (MQ=255)
                                                                tatTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTg                                                                                >  2:670459/1‑148 (MQ=255)
                                                                  tttttCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCg                                                                              >  1:2087476/1‑148 (MQ=255)
                                                                   ttttCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGt                                                                             >  1:884819/1‑148 (MQ=255)
                                                                   ttttCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGt                                                                             >  1:2761148/1‑148 (MQ=255)
                                                                   ttttCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGt                                                                             >  2:3073591/1‑148 (MQ=255)
                                                                    tttCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGtt                                                                            <  2:84852/148‑1 (MQ=255)
                                                                       cTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAg                                                                         >  1:2781803/1‑148 (MQ=255)
                                                                            gaaCGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTAc                                                                    <  1:59230/148‑1 (MQ=255)
                                                                              aCGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCt                                                                  <  1:2794460/148‑1 (MQ=255)
                                                                              aCGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCt                                                                  <  1:1327472/148‑1 (MQ=255)
                                                                               cGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCtt                                                                 <  2:2468989/148‑1 (MQ=255)
                                                                               cGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCtt                                                                 <  2:1152534/148‑1 (MQ=255)
                                                                               cGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCtt                                                                 <  1:3029290/148‑1 (MQ=255)
                                                                                   catcatTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGaa                                                             >  1:3174230/1‑148 (MQ=255)
                                                                                    atcatTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGaaa                                                            <  2:2781803/148‑1 (MQ=255)
                                                                                    atcatTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGaaa                                                            <  1:1514243/148‑1 (MQ=255)
                                                                                     tcatTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAg                                                           <  2:1242568/148‑1 (MQ=255)
                                                                                     tcatTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAg                                                           >  1:893811/1‑148 (MQ=255)
                                                                                     tcatTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAg                                                           >  2:2148177/1‑148 (MQ=255)
                                                                                      catTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAgg                                                          <  1:2071128/148‑1 (MQ=255)
                                                                                        ttAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGtt                                                        >  1:1577205/1‑148 (MQ=255)
                                                                                        ttAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGtt                                                        <  1:855787/148‑1 (MQ=255)
                                                                                         tAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGttt                                                       >  1:806992/1‑148 (MQ=255)
                                                                                           gaagaaGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGttttt                                                     >  2:229418/1‑148 (MQ=255)
                                                                                           gaagaaGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGttttt                                                     <  1:2129528/148‑1 (MQ=255)
                                                                                           gaagaaGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGttttt                                                     >  1:1728054/1‑148 (MQ=255)
                                                                                           gaagaaGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGttttt                                                     >  1:339725/1‑148 (MQ=255)
                                                                                           gaagaaGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGttttt                                                     >  2:3299884/1‑148 (MQ=255)
                                                                                              gaagTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAg                                                  <  2:2557960/148‑1 (MQ=255)
                                                                                              gaagTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAg                                                  >  2:20791/1‑148 (MQ=255)
                                                                                                 gTATTACTGGCGTTCCCCCCCAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAgcgc                                               <  2:2528160/148‑1 (MQ=255)
                                                                                                 gTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAgcgc                                               <  1:1601201/148‑1 (MQ=255)
                                                                                                   aTTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCCGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAATGGTTTTTTATCGCTg                                             >  1:1458127/1‑148 (MQ=255)
                                                                                                   aTTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTg                                             >  1:1908014/1‑148 (MQ=255)
                                                                                                      aCTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCa                                          <  2:1064091/148‑1 (MQ=255)
                                                                                                      aCTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCa                                          >  1:1624638/1‑148 (MQ=255)
                                                                                                        tGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTCAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGAGTACCTTCGAAAGGTTTTTTAGCGCTGGCAgg                                        >  2:2344586/1‑148 (MQ=255)
                                                                                                        tGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAgg                                        <  2:2220868/148‑1 (MQ=255)
                                                                                                              tCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGctcatc                                  >  2:2269499/1‑148 (MQ=255)
                                                                                                               cccccACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATcc                                 <  1:1130543/148‑1 (MQ=255)
                                                                                                                 cccACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCAt                               <  2:908033/148‑1 (MQ=255)
                                                                                                                  ccACAGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATg                              >  1:664886/1‑148 (MQ=255)
                                                                                                                     caGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACa                           >  1:2377831/1‑148 (MQ=255)
                                                                                                                     caGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACa                           <  2:2382818/148‑1 (MQ=255)
                                                                                                                     caGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACa                           <  2:1624638/148‑1 (MQ=255)
                                                                                                                     caGCGACCAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACa                           >  2:1553027/1‑148 (MQ=255)
                                                                                                                            cAGGCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTaaaaa                    >  1:1000368/1‑148 (MQ=255)
                                                                                                                              ggCAGTGGTGTGTCCATACTAGGTCATTGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACg                  <  2:48549/148‑1 (MQ=255)
                                                                                                                              ggCAGTGGTGTGTCCATACTAGGTCATTGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACg                  <  2:48602/148‑1 (MQ=255)
                                                                                                                              ggCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACg                  <  1:1818960/148‑1 (MQ=255)
                                                                                                                               gCAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACgg                 >  1:2989943/1‑148 (MQ=255)
                                                                                                                                cAGTGGTGTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGa                <  2:2693605/148‑1 (MQ=255)
                                                                                                                                    ggtgTGTCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAatca            <  2:2717771/148‑1 (MQ=255)
                                                                                                                                     gtgtgtCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAatcat           >  2:1479741/1‑148 (MQ=255)
                                                                                                                                      tgtgtCCATACTAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAatcatc          >  2:472977/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GCTTGCCTGACCAAATAAATAGTCATAATCAGTTAAATTAAAACGGCTACGCCCGCCCGCCGCATATTTTTCTGAAGAACGCTCATCATTAGAAGAAGTATTACTGGCGTTCCCCCACAGCGACCAGGCAGTGGTGTGTCCATACCAGGTCATGGTGGTGTCAGCCGTAAGTGAGGAGCTTTTCGTATTGCCTGATTGTGCAAGATATCCTGCGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTG  >  CP009273/1799984‑1800273

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: