Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,776,469 A→T 6.4% L948Q (CTG→CAG)  pvdD ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,776,4690AT 6.4% 243.4 / 5.7 78L948Q (CTG→CAG) pvdDnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base A (40/33);  new base T (2/3);  total (42/36)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

GCGCCTTGCGGTCGAGCTTGCCGTTGGGGGTCAGCGGCATGCGTGGAATGAACACCAGGTGGGTCGGCACCATGTAGTCCGGCATCGACTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGCCCGTCGACGTCGATCACGGCCACTTCGCGCACCGCACCGAAGCCCTGCAGGCTGGCCTCGATCTCGCCCAACTCGATGCGGAAACCGCGAATCTTC  >  NC_002947/4776320‑4776618
                                                                                                                                                      |                                                                                                                                                     
gcgcCTTGCGGTCGAGCTTGCCGTTGGGGGTCCGCGGCATGCGTGGAATGAACACCAGGTGGGTCGGCACCATGTAGTCCGGCATCGACTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTgca                                                                                                                                                       >  2:1343393/1‑150 (MQ=255)
gcgcCTTGCGGTCGAGCTTGCCGTTGGGGGTCAGCGGCATGCGTGGCATGAACACCAGGTGGGTCGGCACCATGTAGTCCGGCATCGACTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTgca                                                                                                                                                       <  1:145074/150‑1 (MQ=255)
    cTTGCGGTCGAGCTTGCCGTTGGGGGTCAGCGGCATGCGTGGAATGAACACCAGGTGGGTCGGCACCATGTAGTCCGGCATCGACTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGcagcag                                                                                                                                                   >  2:309348/1‑150 (MQ=255)
     ttGCGGTCGAGCTTGCCGTTGGGGGTCAGCGGCATGCGTGGAATGAACACCAGGTGGGTCGGCACCATGTAGTCCGGCATCGACTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAgg                                                                                                                                                  <  1:1444812/150‑1 (MQ=255)
          gTCGAGCTTGCCGTTGGGGGTCAGCGGCATGCGTGGAATGAACACCAGGTGGGTCGGCACCATGTAGTCCGGCATCGACTTCAATGCGTGGGCCTTGAGAGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCTGCAGGCTgat                                                                                                                                             <  2:1908233/150‑3 (MQ=255)
           tCGAGCTTGCCGTTGGGGGTCAGCGGCATGCGTGGAATGAACACCAGGTGGGTCGGCACCATGTAGTCCGGCATCGACTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCa                                                                                                                                            >  2:1889002/1‑150 (MQ=255)
           tCGAGCTTGCCGTTGGGGGTCAGCGGCATGCGTGGAATGAACACCAGGTGGGTCGGCACCATGTAGTCCGGCATCGACTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCa                                                                                                                                            >  2:2165477/1‑150 (MQ=255)
             gAGCTTGCCGTTGGGGGTCAGCGGGATCGTCTTAATGAACACCAGGTGGGTCGGCACCATGTAGTCCGGCATCGACTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATc                                                                                                                                          <  1:1182230/150‑1 (MQ=255)
             gAGCTTGCCGTTGGGGGTCAGCGGCATGCGTGGAATGAACACCAGGTGGGTCGGCACCATGTAGTCCGGCATCGACTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATc                                                                                                                                          >  2:1783417/1‑150 (MQ=255)
             gAGCTTGCCGTTGGGGGTCAGCGGCATGCGTGGAATGAACACCAGGTGGGTCGGCACCATGTAGTCCGGCATCGACTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATc                                                                                                                                          >  1:1814498/1‑150 (MQ=255)
                cTTGCCGTTGGGGGTCAGCGGCATGCGTGGAATGAACACCAGGTGGGTCGGCACCATGTAGTCCGGCATCGACTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTa                                                                                                                                       <  2:2287777/150‑1 (MQ=255)
                   gCCGTTGGGGGTCAGCGGCATGCGTGGAATGAACACCAGGTGGGTCGGCACCATGTAGTCCGGCATCGACTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACgg                                                                                                                                    >  1:1725239/1‑150 (MQ=255)
                     cGTTGGGGGTCAGCGGCATGCGTGGAATGAACACCAGGTGGGTCGGCACCATGTAGTCCGGCATCGACTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGcc                                                                                                                                  >  1:1719908/1‑149 (MQ=255)
                       ttGGGGGTCAGCGGCATGCGTGGAATGAACACCAGGTGGGTCGGCACCATGTAGTCCGGCATCGACTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACg                                                                                                                                >  2:1025426/1‑150 (MQ=255)
                              tCAGCGGCATGCGTGGAATGAACACCAGGTGGGTCGGCACCATGTAGTCCGGCATCGACTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGc                                                                                                                         >  2:688149/1‑150 (MQ=255)
                                aGCGGCATGCGTGGAATGAACACCAGGTGGGTCGGCACCATGTAGTCCGGCATCGACTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGccc                                                                                                                       >  2:2329094/1‑149 (MQ=255)
                                    gCATGCGTGGAATGAACACCAGGTGGGTCGGCACCATGTAGTCCGGCATCGACTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCa                                                                                                                   <  1:201616/150‑1 (MQ=255)
                                       tGCGTGGAATGAACACCAGGTGGGTCGGCACCATGTAGTCCGGCATCGACTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGcagcag                                                                                                                                                   <  2:2044520/115‑1 (MQ=255)
                                       tGCGTGGAATGAACACCAGGTGGGTCGGCACCATGTAGTCCGGCATCGACTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGcagcag                                                                                                                                                   >  1:2044520/1‑115 (MQ=255)
                                              gAATGAACACCAGGTGGGTCGGCACCATGTAGTCCGGCATCGACTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGTCTGGCATCGTACGCCACGACGTAGCCGCCCAGCTGCTGGc                                                                                                          >  2:1321962/1‑150 (MQ=255)
                                                aTGAACACCAGGTGGGTCGGCACCATGTAGTCCGGCATCGACTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGccc                                                                                                        <  1:2165477/150‑1 (MQ=255)
                                                 tGAACACCAGGTGGGTCGGCACCATGTAGTCCGGCATCGACTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCa                                                                                                       <  2:1762032/150‑1 (MQ=255)
                                                    acacCAGGTGGGTCGGCACCATGTAGTCCGGCATCGACTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCg                                                                                                    <  1:309348/150‑1 (MQ=255)
                                                     cacCAGGTGGGTCGGCACCATGTAGTCCGGCATCGACTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCgg                                                                                                   <  2:1885950/150‑1 (MQ=255)
                                                      acCAGGTGGGTCGGCACCATGTAGTCCGGCATCGACTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGc                                                                                                  >  2:202924/1‑150 (MQ=255)
                                                        cAGGTGGGTCGGCACCATGTAGTCCGGCCTAGACGTCGATGTTTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGccc                                                                                                <  2:2017015/150‑1 (MQ=255)
                                                        cAGGTGGGTCGGCACCATGTAGTCCGGCATCGACTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGTCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGccc                                                                                                <  1:1321962/150‑1 (MQ=255)
                                                        cAGGTGGGTCGGCACCATGTAGTCCGGCATCGACTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGccc                                                                                                >  2:800399/1‑150 (MQ=255)
                                                           gTGGGTCGGCACCATGTAGTCCGGCATCGACTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGCCcgtc                                                                                             <  1:202924/150‑1 (MQ=255)
                                                                tCGGCACCATGTAGTCCGGCATCGACTTCAATGCGTGGGCCTTGAGTACCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGTCTGCGAGCTGCTTGCCCAGCGGCCcgtcgacgt                                                                                        >  2:1232882/1‑150 (MQ=255)
                                                                   gCACCATGTAGTCCGGCATCGACTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCTGCAGGCTgatg                                                                                                                                            <  2:1805274/95‑4 (MQ=255)
                                                                   gCACCATGTAGTCCGGCATCGACTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCTGCAGGCTgatg                                                                                                                                            >  1:1805274/1‑92 (MQ=255)
                                                                    cACCATGTAGTCCGGCATCGACTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATCGCCTGCCAGCTGCTTGCCCAGCGGCCCGTCGACGTCGAt                                                                                    <  1:1232882/150‑1 (MQ=255)
                                                                       cATGTAGTCCGGCATCGACTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGc                                                                                                  <  1:870872/133‑1 (MQ=255)
                                                                       cATGTAGTCCGGCATCGACTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGc                                                                                                  >  2:870872/1‑133 (MQ=255)
                                                                       cATGTAGTCCGGCATCGACTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGCCCGTCGACGTCGATCAc                                                                                 >  1:528860/1‑150 (MQ=255)
                                                                          gTAGTCCGGCATCGACTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGCCCGTCGACTTCGATCACGGc                                                                              <  1:1889002/150‑1 (MQ=255)
                                                                          gTAGTCCGGCATCGACTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGCCCGTCGACGTCGATCACGGc                                                                              <  1:67676/150‑1 (MQ=255)
                                                                                cGGCATCGACTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGGTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGCCCGTAGCCGTCGATCACGGcccattc                                                                        >  2:229670/1‑145 (MQ=255)
                                                                                    aTCGACTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGCCCGTCGACGTCGATCACGGCCACTTCgcgc                                                                    >  2:1315371/1‑150 (MQ=255)
                                                                                        aCTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGCCCGTCGACGTCGATCACg                                                                                >  2:2247658/1‑134 (MQ=255)
                                                                                        aCTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGCCCGTCGACGTCGATCACg                                                                                <  1:2247658/134‑1 (MQ=255)
                                                                                             aaTGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGCCCGTCGACGTCGATCACGGCCACTTCGCGcaccgcaccg                                                           >  2:59716/1‑150 (MQ=255)
                                                                                             aaTGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGCCCGTCGACGTCGATCACGGCCACTTCGCGcaccgcaccg                                                           >  2:1335762/1‑150 (MQ=255)
                                                                                                gCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGCCCGTCGACGTCGATCACGGCCACTTCGCGCACCGCACCGAAg                                                        <  2:715830/150‑1 (MQ=255)
                                                                                                   tGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGc                                                                                                     >  1:471756/1‑102 (MQ=255)
                                                                                                   tGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGc                                                                                                     <  2:471756/102‑1 (MQ=255)
                                                                                                    gggCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGCCCGTCGACGTCGATCACGGCCACTTCGCGCACCGAACCGAAGCCCt                                                    >  1:1236434/1‑150 (MQ=255)
                                                                                                      gCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGCCCGTCGACGTCGATCACGGCCACTTCGCGCACCGCACCGAAGCCCTGc                                                  <  1:1928243/150‑1 (MQ=255)
                                                                                                       ccTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGCCCGTCGACGTCGATCACGGCCACTTCGCGCACCGCACCGAAGCCCTGCa                                                 >  2:1572771/1‑150 (MQ=255)
                                                                                                           gAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGCCCGTCGACGTCGCTCCCGGCCACTTCGCGCGCCGCACCGAAGCCCTGCAGGCt                                             >  2:1331059/1‑150 (MQ=255)
                                                                                                            aGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGCCCGTCGACGTCGATCACGGCCACTTCGCGCACCGCACCGAAGCCCTGCAGGCTg                                            <  2:1821872/150‑1 (MQ=255)
                                                                                                              tGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGCCCGTCGACGTCGATCACGGCCACTTCGCGCACCGCACCGAAGCCCTGCAGGCTGGc                                          <  2:2328231/150‑1 (MQ=255)
                                                                                                                ccAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCa                                                                                                       >  2:2057577/1‑87 (MQ=255)
                                                                                                                ccAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCa                                                                                                       <  1:2057577/87‑1 (MQ=255)
                                                                                                                ccAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGGCCGTCGACGTCGATCACGGCCACTTCGCGCACCGAACCGAAGCCCTGCAGGCTGGCCt                                        <  2:1236434/150‑1 (MQ=255)
                                                                                                                 cAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGCCATCGTACGGCACGATATAGCCTGCCACGTGCTTGCCCAGCGGCCCGTCGACGTCGATCACGGCCACCTCGCGCACCGCACCGAAGCCCTGCAGGCTGGCCTc                                       <  1:1343393/150‑1 (MQ=255)
                                                                                                                  aGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGccc                                                                                                >  2:1143738/1‑92 (MQ=255)
                                                                                                                  aGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGccc                                                                                                <  1:1143738/92‑1 (MQ=255)
                                                                                                                  aGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGCCCGTCGACGTCGATCACGGCCACTTCGCGCACCGCACCGAAGCCCTGCAGGCTGGCCTCg                                      >  1:1976292/1‑150 (MQ=255)
                                                                                                                     cgcAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGCCCGTCGACGTCGATCACGGCCACTTCGCGCACCGCACCGAAGCCCTGCAGGCTGGCCTCGAtc                                   <  1:59716/150‑1 (MQ=255)
                                                                                                                      gcAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGCCCGTCGCCGTCGATCACGGCCACTTCGCGCACCGCACCGAAGCCCTGCAGGCTGGCCTCGAtct                                  >  1:648673/1‑150 (MQ=255)
                                                                                                                        aGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGCCCGTCGACGTCGATCACGGCCACTTCGCGCACCGAACCGTAGCCCTGCAGGCTGGCCTCGATCTCg                                <  1:1025426/150‑1 (MQ=255)
                                                                                                                           agagccTGGACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGCCCGTCGACGTCGATCACGGCCACTTCGcgcaccgcac                                                             <  1:2311567/114‑1 (MQ=255)
                                                                                                                           agagccTGGACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGCCCGTCGACGTCGATCACGGCCACTTCGcgcaccgcac                                                             >  2:2311567/5‑118 (MQ=255)
                                                                                                                            tttCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGCCCGTCGCCGTCGATCCCGGCCACTTCGCGCACCGAACCGAAGCCCCACAGGATGGCGTCGATCTCGcccc                            >  1:1158171/1‑149 (MQ=255)
                                                                                                                                cTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGCCCGTCGACGTCGATCACGGCCACTTCGCGCACAGCACCGCAGCCCTGCAGGCTGGCGTAGGTATCGCCCAAATc                        >  1:1532745/1‑150 (MQ=255)
                                                                                                                                cTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCACCTGCTTGCCCAGCGGCCCGTCGACGTCGACCACGGCCACTCCGCGCACCGCACCGAACCCCTGCAGGCTGGCCTCGATCTCGCCCCACTc                        >  1:1492886/1‑150 (MQ=255)
                                                                                                                                 tGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGCCCGTCGACGTCGATCACgg                                                                               <  1:654030/94‑1 (MQ=255)
                                                                                                                                 tGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGCCCGTCGACGTCGATCACgg                                                                               >  2:654030/1‑94 (MQ=255)
                                                                                                                                  gAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGCCCGTCGACGTCGATCACGGCCACTTCGCGCACCGCACCGAAGCCCTGCAGGCTGGCCTCGATCTCGCCCAACTCGa                      >  2:1439336/1‑150 (MQ=255)
                                                                                                                                  gAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGCCCGTCGACGTCGAGCACGGCCACTTCGCGCACCACACCGAACCCGGACGGGCTGGCGTCTGTCTCGCCCAACTCGa                      >  1:1978952/1‑150 (MQ=255)
                                                                                                                                   aaCCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGATGCTTGCCCAGCGGACCGTAGACGTCGATCACCGCAGCTCGGCGAACCGCACAACAGCGCGTCAGGCTGGCCTCGATCCAGCGCAACTCGag                     >  2:565080/1‑149 (MQ=255)
                                                                                                                                       tCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATCTAGCCTGCGAGCTGCTTGCACAGCGGCCCTCCGACGTCGATCACGGCCACTTGGCGCGCCGCACCGAAGCCCTGCAGGCTGGCCTCGATCCCGCCCAACTCGATGCgg                 >  1:1617460/1‑150 (MQ=255)
                                                                                                                                       tCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGCCCGTCGACGTCGATCACGGCCACTTCGCGCACCGCACCGAAGCCCTGCAGGCTGGCCTCGATCTCGCCCAACTCGATGccg                 >  1:2182564/1‑148 (MQ=255)
                                                                                                                                        cGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGCCCGTCGACGTCGATCACGGCCACTTCGCGCACCGCACCGAAGCCCTGCAGGCTGGCCTCGATCTCGCCCAACTCGATGCGGa                <  2:1725239/150‑1 (MQ=255)
                                                                                                                                        cGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCACCGGCCCGTCGACGTCGAt                                                                                    <  1:2437973/82‑1 (MQ=255)
                                                                                                                                        cGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCCGCTTGCCCACCGGCCCGTCGACGTCGAt                                                                                    >  2:2437973/1‑82 (MQ=255)
                                                                                                                                          atgccAGCCTGCTGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGCCCGTCGACGTCGATCACGGCCACTTCGCGCACCGCACCGAAGCCCTGCAGGCTGGCCTCGATCTCGc                               <  2:1277982/129‑1 (MQ=255)
                                                                                                                                          atgccAGCCTGCTGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGCCCGTCGACGTCGATCACGGCCACTTCGCGCACCGCACCGAAGCCCTGCAGGCTGGCCTCGATCTCGc                               >  1:1277982/5‑133 (MQ=255)
                                                                                                                                           cATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGCCCGTCGACGTCGATCACGGCCACTTCGCGCACCGCACCGAAGCCCTGCAGGCTGGCCTCGATCTCGCCCAACTCGATGCGGAAAc             <  1:1572771/150‑1 (MQ=255)
                                                                                                                                            aTCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGCCCGTCGACGTCGATCACGGCCACTTCGCGCATCGCACCGAAGCCCTGCAGGCTGGCCTCGATCTCGCCCAACTCGATGCGGAAAcc            >  1:1459426/1‑150 (MQ=255)
                                                                                                                                              cAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGCCCGTCGACGTCGATCACGGCCACTTCGCGCACCGCACCGAAGCCCTGCAGGCTGGCCTCGATCTCGCCCAACTCGATGCGGAAACcgc          <  2:648673/150‑1 (MQ=255)
                                                                                                                                                      agcagGCTGGCATCGTACGGCACGATATAGCCTGCCCGGTGGTTGCCCAGCGGCCCGTCGACGTCGATCACGGCCACTTCGCGCACCGCACCGAAGCCCTGCAGGCTGGCCTCGATCTCGCCCAACTCGATGCGGAAACCGCGAATCTTc  <  1:1204891/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GCGCCTTGCGGTCGAGCTTGCCGTTGGGGGTCAGCGGCATGCGTGGAATGAACACCAGGTGGGTCGGCACCATGTAGTCCGGCATCGACTTCAATGCGTGGGCCTTGAGTGCCAGGCGCAGGGTTTCCTGAACCTCGGCATCAGCCTGCAGCAGGCTGGCATCGTACGGCACGATATAGCCTGCCAGCTGCTTGCCCAGCGGCCCGTCGACGTCGATCACGGCCACTTCGCGCACCGCACCGAAGCCCTGCAGGCTGGCCTCGATCTCGCCCAACTCGATGCGGAAACCGCGAATCTTC  >  NC_002947/4776320‑4776618

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