Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 544,652 G→T 6.4% K964N (AAG→AAT) ‡ rpoC → DNA‑directed RNA polymerase subunit beta'

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947544,6520GT 6.4% 248.4 / 4.7 78K964N (AAG→AAT) ‡rpoCDNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base G (37/36);  new base T (2/3);  total (39/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

TCGGCGTTATCGCTGCCCAGTCGATCGGTGAGCCGGGTACCCAGCTGACCATGCGTACGTTCCACATCGGTGGTGCTGCAAGCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGTTCGGTCGTGAGCGTGAGCGCTACAAGCTGCCTTACGGTGCGGTGATTTCGGTCAAGGAAGGTGAAAAGGTCGAAGCTG  >  NC_002947/544504‑544800
                                                                                                                                                     |                                                                                                                                                    
tCGGCGTTATCGCTGCCCAGTCGATCTGTGAGCCGGGTACCCAGCTGACCATGCGTACGTTCCACATCGGTGTTGCTGCAAGCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCTGTATGGTGCGTCTGCACAACCTGAAg                                                                                                                                                      <  8:66300/149‑1 (MQ=255)
  ggCGTTATCGCTGCCCAGTCGATCGGTGAGCCGGGTACCCAGCTGACCATGCGTACGTTCCACATCGGTGGTGCTGCAAGCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCa                                                                                                                                                    <  3:105934/149‑1 (MQ=255)
      ttATCGCTGCACAGTCGATCGGTGAGCCGGGTACCCAGCTGACCATGCGTACGTTCCACATCGGTGGTGCTGCAAGCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGtt                                                                                                                                                <  1:354241/149‑1 (MQ=255)
         tCGCTGCCCAGTCGATCGGTGAGCCGGGTACCCAGCTGACCATGCGTACGTTCCACATCGGTGGTGCTGCAAGCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAg                                                                                                                                             <  2:342335/149‑1 (MQ=255)
           gCTGCCCAGTCGATCGGTGAGCCGGGTACCCAGCTGACCATGCATACGTTCCACATCGGTGGTGCTGCAAGCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCg                                                                                                                                           >  1:72109/1‑149 (MQ=255)
           gCTGCCCAGTCGATCGCTGAGCCGTGTACCCAGCTGACCATGCGTACGTTACACATCGGTGGTGCTTCAAGCAGTACCTCGGCTGCCGAGAGCGTCCAGGTTAAGAATTGCGGTATGGTGCGTCTGCACAACCTGCTTCAGGTTGTGca                                                                                                                                           <  6:377637/149‑2 (MQ=255)
           gCTGCCCAGCCGATCGGTGAGCCGGGTACCCAGCTGACCATGCGTACGTTCCACATCGGTGGTGCTGCAAGCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCTCAACCTGCTTCAGGTTGTGc                                                                                                                                            <  3:264725/148‑1 (MQ=255)
               cccAGTCGATCGGTGAGCCGGGTACCCAGCTGACCATGCGTACGTTCCACATCGGTGGTGCTGCAAGCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGcc                                                                                                                                       >  1:272425/1‑149 (MQ=255)
                ccAGTCGATCGGTGAGCCGGGTACCCAGCTGACCATGCGTACGTTCCACATCGGTGGTGCTGCAAGCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGcc                                                                                                                                       >  4:37002/1‑148 (MQ=255)
                        tCGGTGAGCCGGGTACCCAGCTGACCATGCGTACGTTCCACATCGGTGGTGCTGCAAGCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGCTTCAGGTTgtg                                                                                                                                             <  2:342883/134‑3 (MQ=255)
                        tCGGTGAGCCGGGTACCCAGCTGACCATGCGTACGTTCCACATCGGTGGTGCTGCAAGCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGCTTCAGGTTgtg                                                                                                                                             >  1:342883/1‑132 (MQ=255)
                          ggTGAGCCGGGTACCCAGCTGACCATGCGTACGTTCCACATCGGTGGTGCTGCAAGCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAAcc                                                                                                                             <  1:241537/148‑1 (MQ=255)
                           gTGAGCCGGGTACCCAGCTGACAATGCGTACGTTCAACATCGGTGGTGCTGCAATCCGTACCTCGGCTGCCGACAGCGACAAGGTTAAGAATGGCGGTATGGTGAGTCTGCACAACATGAAGCAGGTTGAGAGTGCCGACGGTACCCTg                                                                                                                           >  8:170887/1‑149 (MQ=255)
                             gAGCCGGGTACCCAGCTGACCATGCGTACGTTCCACATCGGTGGTGCTGCAAGCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGt                                                                                                                         <  3:223161/149‑1 (MQ=255)
                              aGCCGGGTACCCAGCTGACCATGCGTACGTTCCACATCGGTGGTGCTGCAAGCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGtt                                                                                                                        <  5:219449/149‑1 (MQ=255)
                                      aCCCAGCTGACCATGCGTACGTTCCACATCGGTGGTGCTGCAAGCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATc                                                                                                                >  8:39183/1‑149 (MQ=255)
                                      aCCCAGCTGACCATGCGTACGTTCCACATCGGTGGTGCTGCAAGCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATc                                                                                                                >  1:276966/1‑149 (MQ=255)
                                        ccAGCTGACCATGCGTACGTTCCACATCGGTGGTGCTGCAAGCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATcgc                                                                                                              >  4:87313/1‑149 (MQ=255)
                                            cTGACCATGCGTACGTTCCACATCGGTGGTGCTGCAAGCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTc                                                                                                          >  1:272033/1‑149 (MQ=255)
                                             tGACCATGGGTACGTTCCACATCGGTTGTGCTGCAAGCCGTTACTCGGCTTCCGTAAGGGTCCAGGTTAAGAATGGCGGTATGGTGCGTCCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTc                                                                                                          <  1:108605/149‑1 (MQ=255)
                                               aCCATGCGTACGTTCCACATCGGTGGTGCTGCAAGCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCgg                                                                                                        <  4:93731/148‑1 (MQ=255)
                                               aCCATGCGTACGTTCCACATCGGTGGTGCTGCAAGCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCgg                                                                                                        <  3:37002/148‑1 (MQ=255)
                                                 cATGCGTACGTTCCACATCGGTGGTGCTGCAAGCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACg                                                                                                                                   <  6:80193/119‑1 (MQ=255)
                                                 cATGCGTACGTTCCACATCGGTGGTGCTGCAAGCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACg                                                                                                                                   >  5:80193/1‑119 (MQ=255)
                                                      gTACGTTCCACATCGGTGGTGCTGCAAGCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTg                                                                                                >  5:80219/1‑149 (MQ=255)
                                                       tACGTTCCACATCGGTGGTGCTGCAAGCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTgg                                                                                               <  6:322131/149‑1 (MQ=255)
                                                       tACGTTCCACATCGGTGGTGCTGCAAGCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTgg                                                                                               >  5:39106/1‑149 (MQ=255)
                                                          gTTCCACATCGGTGGTGCTGCAAGCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGa                                                                                            >  5:368228/1‑149 (MQ=255)
                                                             ccACATCGCTGGTGCTGAAAGCCGAACCTCGGCTGACGACAGCGTCAAGGTTAATAATCGCGTTATTTTCCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCATATCGCGTTAGGGTGAGTTGTCTAtt                                                                                          >  8:196120/1‑148 (MQ=255)
                                                              cacaTCGGTGGTGCTGCAAGCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTg                                                                                         <  7:39183/148‑1 (MQ=255)
                                                              cacaTCGGTGGTGCTGCAAGCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTg                                                                                         <  3:109187/148‑1 (MQ=255)
                                                               acaTCGGTGGTGCTGCAAGCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGAAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCt                                                                                       >  1:349042/1‑149 (MQ=255)
                                                                caTCGGTGGTGCTGCAAGCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGc                                                                                                                                            <  1:346929/95‑1 (MQ=255)
                                                                caTCGGTGGTGCTGCAAGCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGc                                                                                                                                            >  2:346929/1‑95 (MQ=255)
                                                                 aTCGGTGGTGCTGCAAGCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTATCCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGa                                                                                     >  8:324589/1‑149 (MQ=255)
                                                                  tCGGTGGTGCTGCAAGCCGTACCTCGGCTGCCGACTGCGGCCAGGTTTAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTTAGCGTGACGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGAc                                                                                    <  7:170887/149‑1 (MQ=255)
                                                                      tggtgCTGCAAGCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGt                                                                                >  4:77314/1‑149 (MQ=255)
                                                                       ggtgCTGCAAGCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGtt                                                                               >  6:199384/1‑149 (MQ=255)
                                                                       ggtgCTGCAAGCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGtt                                                                               >  8:218280/1‑149 (MQ=255)
                                                                         tgctgcAAGCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGTTCg                                                                             >  1:196645/1‑149 (MQ=255)
                                                                              cAAGCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGTTCGGTc                                                                          <  4:149644/147‑1 (MQ=255)
                                                                                aGCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGTTCGGTcgt                                                                        >  6:289018/1‑147 (MQ=255)
                                                                                 gCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGTTCGGTcgtgag                                                                     >  7:102004/1‑149 (MQ=255)
                                                                                  ccGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACTTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGTTCGGTcgtgagc                                                                    <  8:102004/149‑1 (MQ=255)
                                                                                  ccGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGAc                                                                                    <  5:207003/133‑1 (MQ=255)
                                                                                  ccGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGAc                                                                                    >  6:207003/1‑133 (MQ=255)
                                                                                    gTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGATGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGt                                                                                                      >  7:190007/1‑113 (MQ=255)
                                                                                    gTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGATGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGt                                                                                                      <  8:190007/113‑1 (MQ=255)
                                                                                      aCCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGGGTTCGGTcgtgagcgt                                                                  >  1:345003/1‑147 (MQ=255)
                                                                                       ccTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGTTCGGTcgtgagcgtgag                                                               <  6:39106/149‑1 (MQ=255)
                                                                                       ccTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGTTCGGTcgtgagcgtga                                                                <  2:276966/148‑1 (MQ=255)
                                                                                        cTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGTTCGGTCgtgagcgtgagc                                                              >  2:313216/1‑149 (MQ=255)
                                                                                         tCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGTTCGGTCGtgagcgtgagcg                                                             >  7:143663/1‑149 (MQ=255)
                                                                                             cTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGTTCGGTCGTGAGCGTGAGCGCTa                                                          >  6:171105/1‑148 (MQ=255)
                                                                                                    cAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGTTCGGTCGTGAGCGTGAGCGCTACAAGCTGc                                                  <  3:77314/149‑1 (MQ=255)
                                                                                                        gTCCAGGTTAAGACTGGCGGTATGGTGCGTCTGCACAACCTGCAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGTTCGGTCGTGAGCGTGAGCGCTACAAGCTGCCTTa                                              <  5:287509/149‑1 (MQ=255)
                                                                                                           cAGGTTCAGATTGGCGGTATGGTGCGTCTGGACAACCTGAAGCAGGGTGAGCGTGCCGACGGTAACGTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGTTCGGTCGTGAGCGTGAGCGCTACAAGCTGCCTTAcg                                            <  5:13111/148‑1 (MQ=255)
                                                                                                           cAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGTTCGGTCGTGAGCGTGAGCGCTACAAGCTGCCTTAcgg                                           >  5:176735/1‑149 (MQ=255)
                                                                                                            aGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGTTCGGTCGTGAGCGTGAGCGCTACAAGCTGCCTTAcggt                                          <  2:272425/149‑1 (MQ=255)
                                                                                                              gTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGtt                                                                               >  8:124191/1‑110 (MQ=255)
                                                                                                                tAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGTTCGGTCGTGAGCGTGAGCGCTACAAGCTGCCTTAcggtgcgg                                      >  7:52700/1‑149 (MQ=255)
                                                                                                                      tGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGTTCGGTCGTGAGCGTGAGCGCTACAAGCTGCCTTACGGTGCGGTGAttt                                >  6:114224/1‑149 (MQ=255)
                                                                                                                        gCGGTATGGTGCGTCTGCACAACCTGAAGCAGGGTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGTTCGGTCGTGAGCGTGAGCGCTACAAGCTGCCTTACGGTGCGGTGATTTCg                              <  6:367745/149‑1 (MQ=255)
                                                                                                                            tATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGTTCGGTCGTGAGCGTGAGCGCTACAAGCTGCCTTACGGTGCGGTGATTTCGGTc                           >  3:319579/1‑148 (MQ=255)
                                                                                                                               ggTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGTTCGGTCGTGAGCGTGAGCGCTACAAGCTGCCTTACGGTGCGGTGATTTCGGTCaag                        <  6:80219/148‑1 (MQ=255)
                                                                                                                                  gCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGTTCGGTCGTGAGCGTGAGCGCTACAAGCTGCCTTACGGTGCGGTGATTTCGGTCaaggaag                    <  4:319579/149‑1 (MQ=255)
                                                                                                                                   cGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGTTCGGTCGTGAGCGTGAGCGCTACAAGCTGCCTTACGGTGCGGTGATTTCGGTCaaggaag                    <  2:72109/148‑1 (MQ=255)
                                                                                                                                   cGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGTTCGGTCGTGAGCGTGAGCGCTACAAGCTGCCTTACGGTGCGGTGATTTCGGTCaaggaag                    <  7:196120/148‑1 (MQ=255)
                                                                                                                                    gTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGTTCGGTCGTGAGCGTGAGCGCTACAAGCTGCCTTACGGTGCGGTGATTTCGGTCAAGGAAGGt                  >  1:54801/1‑149 (MQ=255)
                                                                                                                                         gctcaACCTGCTTCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGTTCGGTCGTGAGCGTGAGCGCTACAAGCTGCCTTACGGTGCGGTGATTTCGGTCAAGGAAGGTGaa               >  1:223042/4‑148 (MQ=255)
                                                                                                                                         gCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGTTCGGTCGTGAGCGTGAGCGCTACAAGCTGCCTTACGGTGCGGTGATTTCGGTCAAGGAAGGTGaa               >  4:361281/1‑148 (MQ=255)
                                                                                                                                          cacaACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGTTCGGTCGTGAGCGTGAGCGCTACAAGCTGCCTTACGGTGCGGTGAttt                                <  3:283869/130‑1 (MQ=255)
                                                                                                                                          cacaACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGTTCGGTCGTGAGCGTGAGCGCTACAAGCTGCCTTACGGTGCGGTGAttt                                >  4:283869/1‑130 (MQ=255)
                                                                                                                                          cacaACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGTTCGGTCGTGAGCGTGAGCGCTACAAGCTGCCTTACGGTGCGGTGATTTCGGTCAAGGAAGGTGaaaa             <  5:199384/149‑1 (MQ=255)
                                                                                                                                                cTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGTTCGGTCGTGAGCGTGAGCGCTACAAGCTGCCTTACGGTGCGGTGATTTCGGTCaa                         >  6:320357/1‑131 (MQ=255)
                                                                                                                                                cTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGTTCGGTCGTGAGCGTGAGCGCTACAAGCTGCCTTACGGTGCGGTGATTTCGGTCaa                         <  5:320357/131‑1 (MQ=255)
                                                                                                                                                   aaGCTGGTTGAGGGTGCCGACGGTACCCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGTTCGGTCGTGAGCGTGAGCGCTACAAGCTGACTTACGGTGCGGTGATTTCGGTCAAGGAAGGTGAAAAGGTCGaa      <  2:349042/147‑1 (MQ=255)
                                                                                                                                                   aaGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGTTCGGTCGTGAGCGTGAGCGCTACAAGCTGCCTTACGGTGCGGTGATTTCGGTCAAGGAAGGTGAAAAGGTCGaa      <  7:218280/147‑1 (MQ=255)
                                                                                                                                                   aaGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGTTCGGTCGTGAGCGTGAGCGCTACAAGCTGCCTTACGGTGCGGTGATTTCGGTCAAGGAAGGTGAAAAGGTCGaa      <  1:313216/147‑1 (MQ=255)
                                                                                                                                                     gCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGTTCGGTCGTGAGCGTGAGCGCTACAAGCTGCCTTACGGTGCGGTGAttt                                >  6:83604/1‑119 (MQ=255)
                                                                                                                                                     gCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGTTCGGTCGTGAGCGTGAGCGCTACAAGCTGCCTTACGGTGCGGTGAttt                                <  5:83604/119‑1 (MQ=255)
                                                                                                                                                     gCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGTTCGGTCGTGAGCGTGAGCGCTACAAGCTGCCTTACGGTGCGGTGATTTCGGTCAAGGAAGGTGAAAAGGTCGAAGCTg  >  7:274835/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCGGCGTTATCGCTGCCCAGTCGATCGGTGAGCCGGGTACCCAGCTGACCATGCGTACGTTCCACATCGGTGGTGCTGCAAGCCGTACCTCGGCTGCCGACAGCGTCCAGGTTAAGAATGGCGGTATGGTGCGTCTGCACAACCTGAAGCAGGTTGAGCGTGCCGACGGTAACCTGGTTGCCGTATCGCGTTCGGGTGAGTTGGCGATTGCTGACGAGTTCGGTCGTGAGCGTGAGCGCTACAAGCTGCCTTACGGTGCGGTGATTTCGGTCAAGGAAGGTGAAAAGGTCGAAGCTG  >  NC_002947/544504‑544800

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: