Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,615,345 G→C 7.6% V572V (GTG→GTC) ‡ PP_4085 → rhs family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,615,3450GC 7.6% 241.6 / 9.1 79V572V (GTG→GTC) ‡PP_4085rhs family protein
Reads supporting (aligned to +/- strand):  ref base G (34/39);  new base C (4/2);  total (38/41)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

ATCAACTGGCTAGCGGTGCTGTCCGCGAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAGGTCTATAAGCACATCGACCCTCAGG  >  NC_002947/4615200‑4615493
                                                                                                                                                  |                                                                                                                                                     
aTCAACTGGCTAGCGGTGCTGTCCGCGAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGAtt                                                                                                                                                    >  3:327272/1‑149 (MQ=255)
aTCAACTGGCTAGCGGTGCTGTCCGCGAATACGAATGGGCGCGCTTCGCCATTACTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGAtt                                                                                                                                                    >  8:288066/1‑149 (MQ=255)
 tCAACTGGCTAGCGGTGCTGTCCGCGAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTc                                                                                                                                                   <  6:315074/149‑1 (MQ=255)
 tCAACTGGCTAGCGGTGCTGTCCGCGAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTc                                                                                                                                                   <  3:50652/149‑1 (MQ=255)
 tCAACTGGCTAGCGGTGCTGTCCGCGAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTc                                                                                                                                                   <  7:90464/149‑1 (MQ=255)
         cTAGCGGTGCTGTCCGCGAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAATCACTCGCCattcg                                                                                                                                           <  7:304046/149‑5 (MQ=255)
          tAGCGGTGCTGTCCGCGAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTggg                                                                                                                                          >  8:281296/1‑149 (MQ=255)
              ggTGCTGTCCGCGAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAg                                                                                                                                      <  7:288066/149‑1 (MQ=255)
               gTGCTGTCCGCGAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGc                                                                                                                                     >  2:129369/1‑149 (MQ=255)
                  cTGTCCGCGAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GAt                                                                                                                                 >  1:210214/1‑149 (MQ=255)
                     tCCGCGAATACGAATGGGCTCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGTGAGTGATTCGCCACTGGGGCAGC‑GATg                                                                                                                                <  7:181612/147‑1 (MQ=255)
                     tCCGCGAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCATGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATg                                                                                                                                <  3:351269/147‑1 (MQ=255)
                     tCCGCGAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGt                                                                                                                              <  5:258181/149‑1 (MQ=255)
                      ccGCGAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCCCTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGt                                                                                                                              >  5:82997/1‑148 (MQ=255)
                        gcgAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGaa                                                                                                                           <  6:138660/149‑1 (MQ=255)
                         cgAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAATCACTCGCCat                                                                                                                                              <  8:49211/130‑2 (MQ=255)
                         cgAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAATCACTCGCCat                                                                                                                                              >  7:49211/1‑129 (MQ=255)
                          gAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAaga                                                                                                                         <  7:61217/149‑1 (MQ=255)
                           aaTACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAagag                                                                                                                        <  8:238921/149‑1 (MQ=255)
                               cGAATGGGCGGGCTTCGCAATGCCTGAAACGCGCCCTACTCCCAAGCGCCGGGACGCAAGGCCTTACCGTTTGCTTCAAGTGCTGTTGACCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTAc                                                                                                                     <  1:309782/149‑1 (MQ=255)
                                gAATGGGCGCGCTTCGCCATTCCTGGAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCg                                                                                                                   >  1:77157/1‑149 (MQ=255)
                                gAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCg                                                                                                                   >  8:351185/1‑149 (MQ=255)
                                gAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCg                                                                                                                   >  5:388168/1‑149 (MQ=255)
                                gAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCg                                                                                                                   >  2:73278/1‑149 (MQ=255)
                                      gcgcgcTTCGCCATTCCTGGAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGa                                                                                                             <  2:77157/149‑1 (MQ=255)
                                      gcgcgcTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGa                                                                                                             <  2:210214/149‑1 (MQ=255)
                                      gcgcgcTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGa                                                                                                             <  6:58559/149‑1 (MQ=255)
                                      gcgcgcTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGa                                                                                                             <  6:312860/149‑1 (MQ=255)
                                      gcgcgcTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGa                                                                                                             <  4:327272/149‑1 (MQ=255)
                                       cgcgcTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAg                                                                                                            <  2:209224/149‑1 (MQ=255)
                                       cgcgcTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAg                                                                                                            <  6:82997/149‑1 (MQ=255)
                                       cgcgcTTAGCCATTACTGAAACGCGACCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGCTGAAGAGTACCgctaagag                                                                                                            >  4:122437/1‑142 (MQ=255)
                                          gcTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCATGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGCTTCGAGTa                                                                                                          >  5:360915/1‑148 (MQ=255)
                                            ttCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACaa                                                                                                       >  3:85740/1‑149 (MQ=255)
                                                   ttCCTGAAACGCGCCCTACTCCCAAGCGCCTTGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGTGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGaa                                                                                                <  1:56805/149‑1 (MQ=255)
                                                   ttCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGt                                                                                                                       <  8:46407/126‑1 (MQ=255)
                                                   ttCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGt                                                                                                                       >  7:46407/1‑126 (MQ=255)
                                                       tGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCg                                                                                            >  6:203605/1‑149 (MQ=255)
                                                         aaaCGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCgg                                                                                           <  6:388168/148‑1 (MQ=255)
                                                         aaaCGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGt                                                                                          <  8:258807/149‑1 (MQ=255)
                                                                     cTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTc                                                                              <  7:281296/149‑1 (MQ=255)
                                                                        ccAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCAc                                                                            <  7:324114/148‑1 (MQ=255)
                                                                        ccAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACg                                                                           >  7:243269/1‑149 (MQ=255)
                                                                              gcCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAAc                                                                     >  7:354904/1‑149 (MQ=255)
                                                                               ccTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAAcc                                                                    <  3:320539/149‑1 (MQ=255)
                                                                                    aCGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGgc                                                               >  4:115975/1‑149 (MQ=255)
                                                                                      ggCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGgccg                                                             >  7:214889/1‑149 (MQ=255)
                                                                                        cACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCg                                                                 >  4:375525/1‑143 (MQ=255)
                                                                                        cACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCg                                                                 <  3:375525/143‑1 (MQ=255)
                                                                                        cACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGgccgcc                                                           <  2:311247/149‑1 (MQ=255)
                                                                                             ccTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGAt                                                      >  7:46446/1‑149 (MQ=255)
                                                                                              cTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATc                                                     <  1:129369/149‑1 (MQ=255)
                                                                                                  cGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCacta                                                 <  4:85740/149‑1 (MQ=255)
                                                                                                  cGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCacta                                                 >  2:27893/1‑149 (MQ=255)
                                                                                                  cGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCacta                                                 >  7:38784/1‑149 (MQ=255)
                                                                                                         cTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCac                                                   >  4:5755/1‑140 (MQ=255)
                                                                                                         cTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCac                                                   <  3:5755/140‑1 (MQ=255)
                                                                                                          tCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGt                                                                                          >  2:185997/1‑100 (MQ=255)
                                                                                                          tCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGAt                                                      >  2:126365/1‑136 (MQ=255)
                                                                                                          tCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑CATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGAt                                                      <  1:126365/136‑1 (MQ=255)
                                                                                                            cACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTAGGCCGCCGTGATCACTACTACTggggg                                       >  8:232574/1‑149 (MQ=255)
                                                                                                                   ggAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGCGGATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTg                                 >  1:294005/1‑149 (MQ=255)
                                                                                                                          ccAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAgg                          >  2:304863/1‑149 (MQ=255)
                                                                                                                            aaCCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAGGt                         >  6:362721/1‑148 (MQ=255)
                                                                                                                            aaCCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAGGt                         <  5:362721/148‑1 (MQ=255)
                                                                                                                              cccGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAGGTCta                      <  7:351185/149‑1 (MQ=255)
                                                                                                                                 gATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAGGTCtat                     >  2:293337/1‑147 (MQ=255)
                                                                                                                                    ccccagtGGCGAATCACTCGCCATTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAGGTCTATAAGca                >  2:304031/7‑149 (MQ=255)
                                                                                                                                     aCGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAGGTCTATAAGca                <  1:73278/148‑1 (MQ=255)
                                                                                                                                      ccagtGGCGAATCACTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAGGTCTATAAGcaca              >  3:382560/5‑149 (MQ=255)
                                                                                                                                      ccagtGGCGAATCACTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAGGTCTATAAGcaca              >  4:130837/5‑149 (MQ=255)
                                                                                                                                      cGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTAc                              >  2:375305/1‑133 (MQ=255)
                                                                                                                                      cGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTAc                              <  1:375305/133‑1 (MQ=255)
                                                                                                                                         aTGGCGAGTGATTCGCCACTGGGGCAGC‑GATGTTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTc                                                                                 >  2:220045/1‑79 (MQ=255)
                                                                                                                                         aTGGCGAGTGATTCGCCACTGGGGCAGC‑GATGTTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTc                                                                                 <  1:220045/79‑1 (MQ=255)
                                                                                                                                               aGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAGGTCTATAAGCACATCGACCCt      >  7:188065/1‑148 (MQ=255)
                                                                                                                                                gTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAGGTCTATAAGCACATCGACCCTCa    <  1:27893/149‑1 (MQ=255)
                                                                                                                                                 tGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAGGTCTATAAGCACATCGACCCTCAg   <  3:122437/149‑1 (MQ=255)
                                                                                                                                                 tGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTAAAACCTGGAAGCCGCTGAAACCCTCGTCACGTAGTTCCTCGGGCGCCGTGATCACTACTACTGGGGGCCACTGTACTAGGTCTATAAGCACATCGACCCTCAg   <  6:360915/149‑1 (MQ=255)
                                                                                                                                                  gATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAGGTCTATAAGCACATCGACCCTCAgg  >  7:266364/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
ATCAACTGGCTAGCGGTGCTGTCCGCGAATACGAATGGGCGCGCTTCGCCATTCCTGAAACGCGCCCTACTCCCAAGCGCCTGGACGGCACGCCCTACCGCTTGCCTCCACTGCTGGAG‑CCCCAACCCGATCACGAATGGCGAGTGATTCGCCACTGGGGCAGC‑GATGGTGAAGAGTACCGGTTCGAGTACAACCTGGAAGCCGGTGAAACCCTCGTCACGGACAACCTCGGCCGCCGTGATCACTACTACTGGGGGCCACTGTACGAGGTCTATAAGCACATCGACCCTCAGG  >  NC_002947/4615200‑4615493

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