Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,243,222 G→A 5.9% P433P (CCG→CCA)  PP_2839 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,243,2220GA 5.9% 208.4 / 2.5 68P433P (CCG→CCA) PP_2839hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (31/33);  new base A (2/2);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

AGCGGGGGCGCGGCATGTACAGCGCCAGCCTTGACCAGGGCCAGCTG‑CTGGCGTTGCGCCAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTGCAGCAGTGCGTCGGGGTGATTCAGGAACAGATGAACCAGGCACCCACGCACATCGACCCCCAGGTGTTGCAGTTCTACTATCAGGCCTTGCAGTTCTGCCG  >  NC_002947/3243080‑3243368
                                                                                                                                               |                                                                                                                                                  
aGCGGGGGCGCGGCATGTACAGCGCCAGCCTTGACCAGGGCCAGCTG‑CTGGCGTTGCGCCAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAg                                                                                                                                              >  8:268652/1‑149 (MQ=255)
aGCGGGGGCGCGGCATGTACAGCGCCAGCCTTGACCAGGGCCAGCTG‑CTGGCGTTGCGCCAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAg                                                                                                                                              <  2:33140/149‑1 (MQ=255)
  cGGGGGCGCGGCATGTACAGCGCCAGCCTTGACCAGGGCCAGCTG‑CTGGCGTTGCGCCAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGc                                                                                                                                            <  2:194163/149‑1 (MQ=255)
      ggCGCGGCATGTACAGCGCCAGCCTTGACCAGGGCCAGCTG‑CTGGCGTTGCGCCAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGc                                                                                                                                        <  7:10738/149‑1 (MQ=255)
      ggCGCGGCATGTACAGCGCCAGCCTTGACCAGGGCCAGCTG‑CTGGCGTTGCGCCAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGc                                                                                                                                        >  5:220849/1‑149 (MQ=255)
      ggCGCGGCATGTACAGCGCCAGCCTTGACCAGGGCCAGCTG‑CTGGCGTTGCGCCAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGc                                                                                                                                        >  8:229269/1‑149 (MQ=255)
       gcgcGGCATGTACAGCGCCAGCCTTGACCAGGGCCAGCTG‑CTGGCGTTGCGCCAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCg                                                                                                                                       <  2:348676/149‑1 (MQ=255)
              aTGTACAGCGCCAGCCTTGACCAGGGCCAGCTG‑CTGGCGTTGCGCCAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGTGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCt                                                                                                                                <  3:185584/149‑1 (MQ=255)
                 tACAGCGCCAGCCTTGACCAGGGCCAGCTG‑CTGGCGTTGCGCCAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGcc                                                                                                                             <  7:220685/149‑1 (MQ=255)
                 tACAGCGCCAGCCTTGACCAGGGCCAGCTG‑CTGGCGTTGCGCCAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGcc                                                                                                                             >  2:100659/1‑149 (MQ=255)
                  aCAGCGCCAGCCTTGACCAGGGCCAGCTGTCTGGCGTTGCGCCAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGcc                                                                                                                             >  1:30237/1‑149 (MQ=255)
                  aCAGCGCCAGCCTTGACCAGGGCCAGCTG‑CTGGCGTTGCGCCAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGccc                                                                                                                            >  3:123698/1‑149 (MQ=255)
                  aCAGCGCCAGCCTTGACCAGGGCCAGCTG‑CTGGCGTTGCGCCAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCATACCAGGCCAGCGAGTCGCTGccc                                                                                                                            >  5:419019/1‑149 (MQ=255)
                     gcgcCAGCCTTGACCAGGGCCAGCTG‑CTGGCGTTGCGCCAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGaa                                                                                                                         <  7:229269/149‑1 (MQ=255)
                     gcgcCAGCCTTGACCAGGGCCAGCTG‑CTGGCGTTGCGCCAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGaa                                                                                                                         <  3:54928/149‑1 (MQ=255)
                            ccTTGACCAGGGCCAGCTG‑CTGGCGTTGCGCCAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCctgc                                                                                                                  >  1:366177/1‑149 (MQ=255)
                             cTTGACCAGGGCCAGCTG‑CTGGCGTTGCGCCAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCctgct                                                                                                                 >  6:181114/1‑149 (MQ=255)
                              ttGACCAGGGCCAGCTG‑CTGGCGTTGCGCCAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCctgctg                                                                                                                >  3:187309/1‑149 (MQ=255)
                                  ccAGGGCCAGCTG‑CTGGCGTTGCGCCAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCa                                                                                                                                               <  4:233412/114‑1 (MQ=255)
                                  ccAGGGCCAGCTG‑CTGGCGTTGCGCCAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCa                                                                                                                                               >  3:233412/1‑114 (MQ=255)
                                     gggCCAGCTG‑CTGGCGTTGCGCCAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGccc                                                                                                         >  4:171075/1‑149 (MQ=255)
                                        ccAGCTG‑CTGGCGTTGCGCCAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTgc                                                                                                      <  2:30237/149‑1 (MQ=255)
                                        ccAGCTG‑CTGGCGTTGCGCCAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTgc                                                                                                      <  6:220849/149‑1 (MQ=255)
                                        ccAGCTG‑CTGGCGTTGCGCCAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTg                                                                                                       >  3:397963/1‑148 (MQ=255)
                                         tcaagg‑ctgGCGCTGCGCCAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTgca                                                                                                     <  2:108841/144‑1 (MQ=255)
                                         cAGCTG‑CTGGCGTTGCGCCAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGcc                                                                                                                             <  7:129519/125‑1 (MQ=255)
                                         cAGCTG‑CTGGCGTTGCGCCAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGcc                                                                                                                             >  8:129519/1‑125 (MQ=255)
                                                ctgGCGTTGCGCCAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCGGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTGCAGCAGt                                                                                                <  5:48443/148‑1 (MQ=255)
                                                ctgGCGTTGCGCCAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTGCAGCAGt                                                                                                <  3:171075/148‑1 (MQ=255)
                                                   gCGTTGCGCCAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGccc                                                                                                         <  8:53721/136‑1 (MQ=255)
                                                   gCGTTGCGCCAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGccc                                                                                                         >  7:53721/1‑136 (MQ=255)
                                                          gcCAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGc                                                                                                                                            >  8:224780/1‑94 (MQ=255)
                                                          gcCAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGc                                                                                                                                            <  7:224780/94‑1 (MQ=255)
                                                          gcCAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTgcagca                                                                                                  <  7:309235/136‑1 (MQ=255)
                                                          gcCAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTgcagca                                                                                                  >  8:309235/1‑136 (MQ=255)
                                                          gcCAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTGCAGCAGTGCGTCGGGGt                                                                                      >  1:134704/1‑148 (MQ=255)
                                                           ccAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTGCAGCAGTGCGTCGGGGTGa                                                                                    >  7:72628/1‑149 (MQ=255)
                                                            cAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTGCAGCAGTGCGTCGGGGTGAt                                                                                   >  8:143912/1‑149 (MQ=255)
                                                             agagCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTGCAGCAGTGCGTCGGGGTGAtt                                                                                  <  3:25372/149‑1 (MQ=255)
                                                                 cAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTGCAGCAGTGCGTCGGGGTGATTCAg                                                                               >  2:126457/1‑148 (MQ=255)
                                                                       gCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTGCAGCAGTGCGTCGGGGTGATTCAGGaa                                                                            >  6:75970/1‑145 (MQ=255)
                                                                       gCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTGCAGCAGTGCGTCGGGGTGATTCAGGAACAGa                                                                        <  1:64895/149‑1 (MQ=255)
                                                                        ccTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTGCAGCAGTGCGTCGGGGTGATTCAGGAACAGAt                                                                       <  7:268652/149‑1 (MQ=255)
                                                                           gggCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTGCAGCAGTGCGTCGGGGTGATTCAGGAACAGATGaa                                                                    >  6:309480/1‑149 (MQ=255)
                                                                              cTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTGCAGCAGTGCGTCGGGGTGATTCAGGAACAGATGAACCa                                                                 <  7:143912/149‑1 (MQ=255)
                                                                              cTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTGCAGCAGTGCGTCGGGGTGATTCAGGAACAGATGAACCa                                                                 <  3:143068/149‑1 (MQ=255)
                                                                               tGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTGCAGCAGTGCGTCGGGGTGATTCAGGAACAGATGAACCAg                                                                <  8:72628/149‑1 (MQ=255)
                                                                                 gTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTGCAGCAGTGCGTCGGGGTGATTCAGGAACAGATGAACCAgg                                                               <  5:75970/148‑1 (MQ=255)
                                                                                 gTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTGCAGCAGTGCGTCGGGGTGATTCAGGAACAGATGAACCAGGc                                                              <  5:309480/149‑1 (MQ=255)
                                                                                          cTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTGCAGCAGTGCGTCGGGGTGATTCAGGAACAGATGAACCAGGCACCCACGca                                                     >  3:224548/1‑149 (MQ=255)
                                                                                              aCCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTGCAGCAGTGCGTCGGGGTGATTCAGGAACAGATGAACCAGGCACCCACGCACATc                                                 >  6:354077/1‑149 (MQ=255)
                                                                                                 gCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCATACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTGCAGCAGTGCGTCGGGGTGATTCAGGAACAGATGAACCAGGCACCCACGCACATCGAc                                              >  5:113314/1‑149 (MQ=255)
                                                                                                        cAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTGCAGCAGTGCGTCGGGGTGATTCAGGAACAGATGAACCAGGCACCCACGCACATCGAcccccggt                                       >  7:203045/1‑146 (MQ=255)
                                                                                                          ggCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTGCAGCAGTGCGTCGGGGTGATTCAGGAACAGATGAACCAGGCACCCACGCACATCGAccccc                                          <  1:262509/144‑1 (MQ=255)
                                                                                                          ggCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTGCAGCAGTGCGTCGGGGTGATTCAGGAACAGATGAACCAGGCACCCACGCACATCGACCCCCAGgtg                                     >  6:108099/1‑149 (MQ=255)
                                                                                                          ggCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCAGCAGCAGTCCGTCGCGGAGAATCAGGAACAGATCAACCAGACACCAACGAACAACGACACCCaacat                                     >  2:262509/1‑145 (MQ=255)
                                                                                                             aGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTGCAGCAGTGCGTCGGGGTGATTCAGGAACAGATGAACCAGGCACCCACGCACATCGACCCCCAGGTGTTg                                  <  4:187309/149‑1 (MQ=255)
                                                                                                             aGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCATACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTGCAGCAGTGCGTCGGGGTGATTCAGGAACAGATGAACCAGGCACCCACGCACATCGACCCCCAGGTGTTg                                  <  6:419019/149‑1 (MQ=255)
                                                                                                              gTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTGCAGCAGTGCGTCGGGGTGATTCAGGAACAGATGAACCAGGCACCCACGCACATCGACCCCCAGGTGTTGc                                 <  4:397963/149‑1 (MQ=255)
                                                                                                                  aaCGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCACGCCCTGCAGCAGTGCGTCGGGGTGATTCAGGAACAGATGAACCAGGCACCCACGCACATCGACCCCCAGGTGTTGCAGtt                             <  4:254647/149‑1 (MQ=255)
                                                                                                                    cgcgCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTGCAGCAGTGCGTCGGGGTGATTCAGGAACAGATGAACCAGGCACCCACGCACATCGACCCCCAGGTGTTGCAGTTct                           <  4:224548/149‑1 (MQ=255)
                                                                                                                          gTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTGCAGCAGTGCGTCGGGGTGATTCAGGAACAGATGAACCAGGCACCCACGCACATCGACCCCCAGGTGTTGCAGTTCTACTAt                      >  2:149560/1‑148 (MQ=255)
                                                                                                                          gTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTGCAGCAGTGCGTCGGGGTGATTCAGGAACAGATGAACCAGGCACCCACGCACATCGACCCCCAGGTGTTGCAGTTCTACTAt                      >  1:297107/1‑148 (MQ=255)
                                                                                                                           tataAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTGCAGCAGTGCGTCGGGGTGATTCAGGAACAGATGAACCAGGCACCCACGCACATCGACCCCC‑GGTGTTGCAGTTCTACTATc                     <  8:203045/147‑1 (MQ=255)
                                                                                                                                   gcagcGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTGCAGCAGTGCGTCGGGGTGATTCAGGAACAGATGAACCAGGCCCCCACGCACATCGACCCCCAGGTGTTGCAGTTCTACTATCAGGCCTTg             >  1:408539/1‑148 (MQ=255)
                                                                                                                                   gcagcGTGCGCCATACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTGCAGCAGTGCGTCGGGGTGATTCAGGAACAGATGAACCAGGCACCCACGCACATCGACCCCCAGGTGTTGCAGTTCTACTATCAGGCCTTGc            <  6:113314/149‑1 (MQ=255)
                                                                                                                                           cgcCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTGCAGCAGTGCGTCGGGGTGATTCAGGAACAGATGAACCAGGCACCCACGCACATCGACCCCCAGGTGTTGCAGTTCTACTATCAGGCCTTGCAGTTCTGc    >  7:308536/1‑149 (MQ=255)
                                                                                                                                             ccGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTGCAGCAGTGCGTCGGGGTGATTCAGGAACAGATGAACCAGGCACCCACGCACATCGACCCCCAGGTGTTGCAGTTCTACTATCAGGCCTTGCAGTTCTGCCg  <  1:149560/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
AGCGGGGGCGCGGCATGTACAGCGCCAGCCTTGACCAGGGCCAGCTG‑CTGGCGTTGCGCCAGAGCAAACCGCCTGGGCTGGTCAGTGCGCTGGACCGCCTGAACAGGCAGTGGAACGCGCTGTATAAAGAGCAGCGTGCGCCGTACCAGGCCAGCGAGTCGCTGCCCGAAGGCCTGCTGCGCGCCCTGCAGCAGTGCGTCGGGGTGATTCAGGAACAGATGAACCAGGCACCCACGCACATCGACCCCCAGGTGTTGCAGTTCTACTATCAGGCCTTGCAGTTCTGCCG  >  NC_002947/3243080‑3243368

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: ‑