Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,722,806 T→A 6.9% Y21* (TAT→TAA) ‡ PP_3288 → universal stress protein family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,722,8060TA 6.9% 212.9 / 5.3 72Y21* (TAT→TAA) ‡PP_3288universal stress protein family protein
Reads supporting (aligned to +/- strand):  ref base T (35/32);  new base A (2/3);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.67e-01

GTGCCTTGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTG  >  NC_002947/3722661‑3722951
                                                                                                                                                 |                                                                                                                                                 
gTGCCTTGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGAGCCGAACCCCGGTCTATgcc                                                                                                                                                >  7:308697/1‑147 (MQ=255)
gTGCCTTGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCATCCCCCCTGATGCGCCGAACCCCGGTCTATGAc                                                                                                                                                >  3:28072/1‑149 (MQ=255)
    cTTGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCggg                                                                                                                                            <  7:54087/149‑1 (MQ=255)
        aGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGgctgc                                                                                                                                        <  8:11961/149‑1 (MQ=255)
         gTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCATAGACCggggttc                                                                                                                                        <  2:191047/148‑5 (MQ=255)
          ttGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCGAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGctgctg                                                                                                                                      >  8:311681/1‑149 (MQ=255)
            gATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTgcg                                                                                                                                    >  6:179381/1‑149 (MQ=255)
             aTAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTgcg                                                                                                                                    <  5:179381/148‑1 (MQ=255)
              tAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCt                                                                                                                                  <  2:304424/149‑1 (MQ=255)
               aaaCATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCATAGACCggggttcg                                                                                                                                       >  4:231046/1‑138 (MQ=255)
               aaaCATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCATAGACCggggttcg                                                                                                                                       <  3:231046/143‑6 (MQ=255)
                aaCATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGCTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTCATGCGACGAACCCCGGTCTATGACCGCGCTGCTGCGCTgg                                                                                                                                >  1:81974/1‑149 (MQ=255)
                 aCATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTggc                                                                                                                               >  5:87192/1‑149 (MQ=255)
                       aGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTggccaaggc                                                                                                                         >  6:245709/1‑149 (MQ=255)
                        gTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTggccaaggcc                                                                                                                        <  6:87192/149‑1 (MQ=255)
                             cAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGccaaggccaagg                                                                                                                    >  1:338148/1‑148 (MQ=255)
                                cGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCGAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAAt                                                                                                                <  7:311681/149‑1 (MQ=255)
                                   gtAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGgc                                                                                                             >  8:202265/1‑149 (MQ=255)
                                      gCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGgcgc                                                                                                           <  5:143502/148‑1 (MQ=255)
                                       cTGGCAGATACTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCTTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTg                                                                                                         <  8:266202/149‑1 (MQ=255)
                                        tGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTg                                                                                                         >  6:294607/1‑148 (MQ=255)
                                                gCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATg                                                                                                               <  3:4110/134‑1 (MQ=255)
                                                gCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATg                                                                                                               >  4:4110/1‑134 (MQ=255)
                                                              tCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCAGCAGCCCGGTCATAGACCGGGGTTCGgcgcatcag                                                                                                                              >  3:38247/1‑100 (MQ=21)
                                                              tCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCAGCAGCCCGGTCATAGACCGGGGTTCGgcgcatcag                                                                                                                              <  4:38247/105‑6 (MQ=21)
                                                                gCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTg                                                                                 >  5:198535/1‑148 (MQ=255)
                                                                 cAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTgg                                                                                >  6:15109/1‑148 (MQ=255)
                                                                  aaCCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTgg                                                                                >  7:145478/1‑147 (MQ=255)
                                                                      tCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGc                                                                          >  3:82243/1‑149 (MQ=255)
                                                                       cAGAGGTGTGCCATGGAAACTCCGCAAACACTGCTCCTGATCGCCTCCCACCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAACGCCAAGGGAATGGCGCTGCACATCGACGATTTCGACTACGTGGAAGGGCt                                                                         >  6:22026/1‑149 (MQ=255)
                                                                        agagGTGTCCCAGGGAAACGCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCAGAACCCCGGTCTATGACCGGGCTGCTGCGATGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCtt                                                                        <  8:308697/149‑1 (MQ=255)
                                                                           ggTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGcc                                                                     <  5:15109/149‑1 (MQ=255)
                                                                                 ccATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCAcggcg                                                               <  4:214171/149‑1 (MQ=255)
                                                                                  cATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCAcggcgg                                                              <  7:202265/149‑1 (MQ=255)
                                                                                   aTGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCggg                                                             >  4:120605/1‑149 (MQ=255)
                                                                                       aaaCTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTg                                                         >  3:34009/1‑149 (MQ=255)
                                                                                       aaaCTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTg                                                         >  2:100230/1‑149 (MQ=255)
                                                                                         aCTCCGCAAAGACTGCTGCTGATCTCCTCGCCCGTGATGCGCAGAACGCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTGGACTACGTGGAAGGGCTTGCCACGGCTGTCCTGGt                                                       <  5:22026/149‑1 (MQ=255)
                                                                                           tCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTg                                                      >  7:56859/1‑148 (MQ=255)
                                                                                           tCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGa                                                     <  8:135316/149‑1 (MQ=255)
                                                                                              gCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCggg                                                             >  1:23254/1‑138 (MQ=255)
                                                                                              gCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCggg                                                             <  2:23254/138‑1 (MQ=255)
                                                                                                aaaGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGAGTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGAc                                                >  2:94719/1‑149 (MQ=255)
                                                                                                    aCTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAgg                                            <  7:135682/149‑1 (MQ=255)
                                                                                                    aCTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAgg                                            <  7:140264/149‑1 (MQ=255)
                                                                                                      tgctgctgATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTg                                                         <  5:318217/134‑1 (MQ=255)
                                                                                                      tgctgctgATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTg                                                         >  6:318217/1‑134 (MQ=255)
                                                                                                      tgctgctgATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGgcg                                          >  8:93639/1‑149 (MQ=255)
                                                                                                      tgctgctgATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGgcg                                          >  4:186000/1‑149 (MQ=255)
                                                                                                           ctgATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTgg                                      >  6:310216/1‑148 (MQ=255)
                                                                                                              aTCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCa                                    <  3:120605/147‑1 (MQ=255)
                                                                                                                 gCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCAt                               <  8:145478/149‑1 (MQ=255)
                                                                                                                 gCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCAt                               <  7:93639/149‑1 (MQ=255)
                                                                                                                 gCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCAt                               <  3:186000/149‑1 (MQ=255)
                                                                                                                   cTCCCCCCTGATGCGCCGAACCCCGGTCTATGACGGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATgc                             <  2:338148/149‑1 (MQ=255)
                                                                                                                       ccccTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGa                         <  4:34009/149‑1 (MQ=255)
                                                                                                                         ccTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAgg                       <  6:198535/149‑1 (MQ=255)
                                                                                                                         ccTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCAGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGTCTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAgg                       <  5:294607/149‑1 (MQ=255)
                                                                                                                           tGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGt                     >  5:169944/1‑149 (MQ=255)
                                                                                                                            gATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGt                     >  3:32263/1‑148 (MQ=255)
                                                                                                                                 gcCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTg               <  1:100230/149‑1 (MQ=255)
                                                                                                                                 gcCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGCGCTTGCCACGGCGGGCCTGGTGAATGAAAAGGCGCTGGCAGTAATGCGCGAGGCTTACGtt               >  3:308283/1‑148 (MQ=255)
                                                                                                                                  ccGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCAcgg                                                                 >  2:199575/1‑98 (MQ=255)
                                                                                                                                  ccGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCAcgg                                                                 <  1:199575/98‑1 (MQ=255)
                                                                                                                                  ccGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTg               >  8:278084/1‑148 (MQ=255)
                                                                                                                                    gAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCGAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTg               >  2:282044/1‑146 (MQ=255)
                                                                                                                                    gAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTgcag            >  1:187921/1‑149 (MQ=255)
                                                                                                                                     aaCCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTgcag            >  1:319782/1‑148 (MQ=255)
                                                                                                                                      aCCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTgcagca          <  4:28072/149‑1 (MQ=255)
                                                                                                                                         caggTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTgcagcagca       <  4:82243/147‑1 (MQ=255)
                                                                                                                                          cGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGa                                                     >  8:231327/1‑102 (MQ=255)
                                                                                                                                             tCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGt   <  8:56859/149‑1 (MQ=255)
                                                                                                                                                 tGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGt   >  8:270030/1‑145 (MQ=255)
                                                                                                                                                 tGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGt   <  7:270030/145‑1 (MQ=255)
                                                                                                                                                 tGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTg  <  5:245709/146‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GTGCCTTGAGTTGATAAACATCAAGTTCGCAACGTGTAGCTGGCAGATGCTCTACCTACGTATCGCAACCTCAGAGGTGTGCCATGGAAACTCCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTG  >  NC_002947/3722661‑3722951

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