Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,860,183 G→T 89.3% R183S (CGC→AGC)  mipA ← scaffolding protein for murein synthesizing machinery

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,860,1830GT89.3% 446.2 / 53.0 187R183S (CGC→AGC) mipAscaffolding protein for murein synthesizing machinery
Reads supporting (aligned to +/- strand):  ref base G (11/9);  new base T (80/87);  total (91/96)
Fisher's exact test for biased strand distribution p-value = 6.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

Only 100 of 188 total aligned reads displayed.

TCGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCGAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAGGCCATATCCCAGACGATGCCGTTGCTG  >  CP009273/1860039‑1860325
                                                                                                                                                |                                                                                                                                              
tCGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGc                                                                                                                                             >  1:3019833/1‑148 (MQ=255)
tCGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGc                                                                                                                                             >  1:2737795/1‑148 (MQ=255)
 cGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCg                                                                                                                                            >  2:1984854/1‑148 (MQ=255)
   ggCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGa                                                                                                                                          >  1:510455/1‑148 (MQ=255)
    gCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGACGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCGAGCGGAc                                                                                                                                         >  2:319351/1‑148 (MQ=255)
      tGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGActc                                                                                                                                       >  1:2438613/1‑148 (MQ=255)
       gTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGActct                                                                                                                                      <  1:1320600/148‑1 (MQ=255)
       gTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGActct                                                                                                                                      >  1:2377511/1‑148 (MQ=255)
       gTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCGAGCGGActct                                                                                                                                      <  1:1285214/148‑1 (MQ=255)
       gTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCGAGCGGActct                                                                                                                                      <  1:2413330/148‑1 (MQ=255)
       gTCAGTAACTTAATCAGACAGACGTGTGTAGCGTGCTGTACCGTAAACACTCCAGTCGCCGATGAAGTTGTAGCTGGCGCTCAGCTCCATGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGActct                                                                                                                                      <  1:1124969/148‑1 (MQ=255)
        tCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCtt                                                                                                                                     <  1:410760/148‑1 (MQ=255)
         cAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCttt                                                                                                                                    <  2:3230278/148‑1 (MQ=255)
         cAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCttt                                                                                                                                    <  2:2929902/148‑1 (MQ=255)
          aGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTg                                                                                                                                   <  2:1148031/148‑1 (MQ=255)
            tAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTCGCGGACTCTTTgcg                                                                                                                                 >  2:2596933/1‑148 (MQ=255)
            tAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTgcg                                                                                                                                 <  1:950368/148‑1 (MQ=255)
            tAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTgcg                                                                                                                                 >  1:2415525/1‑148 (MQ=255)
             aaCTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTgcgc                                                                                                                                >  1:1772153/1‑148 (MQ=255)
               cTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGa                                                                                                                              <  1:2136837/148‑1 (MQ=255)
                     cagacagaCGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGcc                                                                                                                        <  2:863280/148‑1 (MQ=255)
                     cagacagaCGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGcc                                                                                                                        >  1:2896435/1‑148 (MQ=255)
                     cagacagaCGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGcc                                                                                                                        >  1:2673314/1‑148 (MQ=255)
                      agacagacGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCa                                                                                                                       >  2:1842957/1‑148 (MQ=255)
                        acagacGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCata                                                                                                                     <  1:1650763/148‑1 (MQ=255)
                         cagacGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCataa                                                                                                                    <  2:119043/148‑1 (MQ=255)
                           gacGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCataata                                                                                                                  <  2:419901/148‑1 (MQ=255)
                           gacGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCataata                                                                                                                  <  1:2797046/148‑1 (MQ=255)
                            acGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAg                                                                                                                 <  2:1545884/148‑1 (MQ=255)
                             cGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGt                                                                                                                <  2:2781500/148‑1 (MQ=255)
                              gggTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTa                                                                                                               <  1:2528286/148‑1 (MQ=255)
                                gtgtAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTAtt                                                                                                             <  2:1930546/148‑1 (MQ=255)
                                    aGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGtt                                                                                                         <  1:1006013/148‑1 (MQ=255)
                                     gCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTc                                                                                                        >  2:123314/1‑148 (MQ=255)
                                     gCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTc                                                                                                        <  1:2896075/148‑1 (MQ=255)
                                       gTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTg                                                                                                      <  1:217264/148‑1 (MQ=255)
                                        tGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTgg                                                                                                     >  1:3220259/1‑148 (MQ=255)
                                        tGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCGAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTgg                                                                                                     >  1:986845/1‑148 (MQ=255)
                                        tGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCGAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTgg                                                                                                     >  1:2640909/1‑148 (MQ=255)
                                         gCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGt                                                                                                    <  2:1801494/148‑1 (MQ=255)
                                         gCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGt                                                                                                    <  1:2814318/148‑1 (MQ=255)
                                         gCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGt                                                                                                    <  1:2814308/148‑1 (MQ=255)
                                                     acacTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCGAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTcc                                                                                        >  1:2617958/1‑148 (MQ=255)
                                                        cTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCGAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACt                                                                                     >  2:774613/1‑148 (MQ=255)
                                                         tCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATATCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTg                                                                                    >  2:1647184/1‑148 (MQ=255)
                                                          ccAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGc                                                                                   <  2:2429321/148‑1 (MQ=255)
                                                            aGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGcac                                                                                 <  2:3104448/148‑1 (MQ=255)
                                                              tCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCacac                                                                               >  1:1008947/1‑148 (MQ=255)
                                                              tCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCacac                                                                               >  2:628502/1‑148 (MQ=255)
                                                               cGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACAcc                                                                              <  2:3281196/148‑1 (MQ=255)
                                                                gCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCa                                                                             <  1:266765/148‑1 (MQ=255)
                                                                 ccGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCaa                                                                            >  2:353756/1‑148 (MQ=255)
                                                                 ccGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCaa                                                                            <  2:1064833/148‑1 (MQ=255)
                                                                       aaGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCg                                                                      <  1:2990005/148‑1 (MQ=255)
                                                                         gTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCGAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGa                                                                    <  1:1523608/148‑1 (MQ=255)
                                                                          ttGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAg                                                                   >  2:1842247/1‑148 (MQ=255)
                                                                          ttGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAg                                                                   <  2:2981544/148‑1 (MQ=255)
                                                                            gTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTc                                                                 <  1:1656589/148‑1 (MQ=255)
                                                                               gCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACg                                                              <  2:2737795/148‑1 (MQ=255)
                                                                               gCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACg                                                              >  2:2150667/1‑148 (MQ=255)
                                                                                 tGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGt                                                            <  2:2399756/148‑1 (MQ=255)
                                                                                 tGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGt                                                            <  2:2046038/148‑1 (MQ=255)
                                                                                 aggCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGt                                                            >  2:484265/2‑148 (MQ=255)
                                                                                      cTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATATCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGc                                                       <  1:2936421/148‑1 (MQ=255)
                                                                                      cTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGc                                                       >  1:2695788/1‑148 (MQ=255)
                                                                                      cTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCCCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTGTTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGc                                                       >  1:2257536/1‑148 (MQ=255)
                                                                                         aGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCAc                                                    >  2:2449051/1‑148 (MQ=255)
                                                                                          gCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCAcc                                                   >  1:2081040/1‑148 (MQ=255)
                                                                                           cTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCg                                                  >  1:61964/1‑148 (MQ=255)
                                                                                                ggTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGgt                                             <  2:2521045/148‑1 (MQ=255)
                                                                                                 gTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGgtg                                            >  1:2339081/1‑148 (MQ=255)
                                                                                                 gTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGgtg                                            >  2:2164323/1‑148 (MQ=255)
                                                                                                 gTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGgtg                                            >  1:1974908/1‑148 (MQ=255)
                                                                                                 gTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCGAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGACACGGTCAGGCCACCGTTGgtg                                            >  1:2691186/1‑148 (MQ=255)
                                                                                                   aaGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTa                                          <  2:3019833/148‑1 (MQ=255)
                                                                                                   aaGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTa                                          <  2:1424779/148‑1 (MQ=255)
                                                                                                   aaGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTa                                          <  2:3019827/148‑1 (MQ=255)
                                                                                                    aGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTaa                                         <  2:1626344/148‑1 (MQ=255)
                                                                                                    aGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTaa                                         <  1:1315343/148‑1 (MQ=255)
                                                                                                     gggCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAAc                                        >  1:2457688/1‑148 (MQ=255)
                                                                                                     gggCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAAc                                        >  1:798315/1‑148 (MQ=255)
                                                                                                     gggCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAAc                                        >  2:810451/1‑148 (MQ=255)
                                                                                                     gggCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCGGCTAGCGGACTCTTTGCGCGATACGCCAAAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACCGTCAGGCCACCGTTGGTGTAAc                                        >  1:908776/1‑148 (MQ=255)
                                                                                                      ggCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACg                                       <  2:2415531/148‑1 (MQ=255)
                                                                                                        cTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGAt                                     <  2:346670/148‑1 (MQ=255)
                                                                                                          ccAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATAc                                   >  2:1532232/1‑148 (MQ=255)
                                                                                                              cTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAAcc                               >  1:1867485/1‑148 (MQ=255)
                                                                                                              cTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAAcc                               >  1:196367/1‑148 (MQ=255)
                                                                                                                  cGTTCGGGTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAGGc                           <  2:941807/148‑1 (MQ=255)
                                                                                                                       gggTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAGGCCatat                      <  1:690264/148‑1 (MQ=255)
                                                                                                                        ggTTATAGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAGGCCATATc                     <  1:1212072/148‑1 (MQ=255)
                                                                                                                           tataGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAGGCCATATCCCa                  >  1:336073/1‑148 (MQ=255)
                                                                                                                             taGCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAGGCCATATCCCAGa                <  2:3075102/148‑1 (MQ=255)
                                                                                                                               gCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACGAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAGGCCATATCCCAGACg              <  1:3246735/148‑1 (MQ=255)
                                                                                                                               gCCACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAGGCCATATCCCAGACg              <  1:3246725/148‑1 (MQ=255)
                                                                                                                                ccACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAGGCCATATCCCAGACGa             <  2:680014/148‑1 (MQ=255)
                                                                                                                                ccACGCAGACCGCTGCGAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAGGCCATATCCCAGACGa             <  2:1269833/148‑1 (MQ=255)
                                                                                                                                 cACGCAGACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAGGCCATATCCCAGACGAt            >  1:823231/1‑148 (MQ=255)
                                                                                                                                       gACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAGGCCATATCCCAGACGATGCCGtt      <  2:2218903/148‑1 (MQ=255)
                                                                                                                                       gACCGCTGCTAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAGGCCATATCCCAGACGATGCCGtt      >  1:892408/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                              
TCGGGCTGTCAGTAACTTCATCAGACAGACGGGTGTAGCGTGCGGTACCGTAAACACTCCAGTCGCCGAGGAAGTTGTAGCTGGCGCTCAGCTCCAGGTAAGGGCTCCAGCTGTCGTTCGGGTTATAGCCACGCAGACCGCTGCGAGCGGACTCTTTGCGCGATACGCCATAATAGTATTCGTTCTGGTTTTCGCTGTTCCACTGCACACCAATACCCGGAGTCACGGTCAGGCCACCGTTGGTGTAACGATACAACCAGGCCATATCCCAGACGATGCCGTTGCTG  >  CP009273/1860039‑1860325

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: