Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 519,677 T→C 5.7% intergenic (+31/‑69) argC → / → erpA N‑acetyl‑gamma‑glutamyl‑phosphate reductase/iron‑sulfur cluster insertion protein ErpA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947519,6770TC 5.7% 258.5 / 4.2 87intergenic (+31/‑69)argC/erpAN‑acetyl‑gamma‑glutamyl‑phosphate reductase/iron‑sulfur cluster insertion protein ErpA
Reads supporting (aligned to +/- strand):  ref base T (39/43);  new base C (3/2);  total (42/45)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

GGTGATCGACAACCTGGTCAAGGGTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCTACCGCGCCGGGAGAATCAACATGAGCGTCGAAACCTTCACCCCTACCGGT  >  NC_002947/519535‑519823
                                                                                                                                              |                                                                                                                                                   
ggTGATCGACAACCTGGTCAAGGGTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGc                                                                                                                                              >  2:241363/1‑149 (MQ=255)
    aTCGACAACCTGGTCAAGGGTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCttt                                                                                                                                          >  7:151004/1‑149 (MQ=255)
      atccgatcTGTTCAAGGGTGCCTCCGTCCATGCGGTGCAGAAACTCAACATTCTGTTCGGTGTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGTCCGCCAGGATTGGAAAAGGCCCTCGTCGTGCGGGCCtt                                                                                                                                           <  5:213648/139‑1 (MQ=255)
       gACAACCTGGTCAAGGGTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTg                                                                                                                                       <  1:241363/149‑1 (MQ=255)
       gACAACCTGGTCAAGGGTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTATGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGACAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTg                                                                                                                                       >  8:167842/1‑149 (MQ=255)
         cAACCTGGTCAAGGGTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCACGACGCGGGCCttttccaa                                                                                                                                     <  3:222509/149‑5 (MQ=255)
          aaCCTGGTCAAGGGTGCCTCCGGCCAGTCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGCTGAGCGCATGGGCCTGTCCCATGGCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGTCCCGCGTCGTGCGGGCCTTTTTGCAg                                                                                                                                    <  7:145855/149‑1 (MQ=255)
          aaCCTGGTCAAGGGTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAg                                                                                                                                    >  7:34230/1‑149 (MQ=255)
            ccTGGTCAAGGGTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTGGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCACGACGCGGGCCttttc                                                                                                                                        <  2:155654/143‑2 (MQ=255)
            ccTGGTCAAGGGTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAgg                                                                                                                                   >  2:221011/1‑148 (MQ=255)
            ccTGGTCAAGGGTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGt                                                                                                                                  <  2:249167/149‑1 (MQ=255)
            ccTGGTCAAGGGTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCACGACGCGGGCCttttc                                                                                                                                        >  1:155654/1‑142 (MQ=255)
             cTGGTCAAGGGTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTc                                                                                                                                 >  2:292803/1‑149 (MQ=255)
              tGGTCAAGGGTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTcg                                                                                                                                >  1:93410/1‑149 (MQ=255)
                gTCAAGGGTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCttt                                                                                                                                          >  1:10078/1‑137 (MQ=255)
                gTCAAGGGTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCttt                                                                                                                                          <  2:10078/137‑1 (MQ=255)
                gTCAAGGGTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTcgcg                                                                                                                              <  8:199651/149‑1 (MQ=255)
                gTCAAGGGTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTcgc                                                                                                                               <  1:162922/148‑1 (MQ=255)
                gTCAAGGGTGCCTCCGGCCAGGAGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGACGCCAGAATTGGAAAAGTCCCGCGTCGTGCGGGCCtt                                                                                                                                           >  6:213648/1‑136 (MQ=255)
                      ggTGCCTCCGGCCAGGCGGTGCAGAACCTCAAGATTCTTTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTaaa                                                                                                                        <  4:171851/149‑1 (MQ=255)
                       gTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAg                                                                                                                       >  2:336422/1‑149 (MQ=255)
                       gTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAg                                                                                                                       <  8:34230/149‑1 (MQ=255)
                           cTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGa                                                                                                                             >  4:5960/1‑139 (MQ=255)
                           cTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGa                                                                                                                             <  3:5960/139‑1 (MQ=255)
                            tCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGca                                                                                                                  >  7:56190/1‑149 (MQ=255)
                            tCCGGCCAGGCGGTGCAGAACCTCAACATTCTCTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGca                                                                                                                  >  7:63808/1‑149 (MQ=255)
                                   aGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTc                                                                                                            <  1:221011/149‑1 (MQ=255)
                                     gCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCtt                                                                                                          >  2:244095/1‑148 (MQ=255)
                                           cAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGAt                                                                                                   >  6:90397/1‑149 (MQ=255)
                                            aGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAgg                                                                                                                                   <  7:118229/116‑1 (MQ=255)
                                            aGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAgg                                                                                                                                   >  8:118229/1‑116 (MQ=255)
                                             gAACCTCAACATCCTGTTCGGGCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCGGCCGTAAGGCCGCCAGCATTCGAAACGGCCCGCGTCGTGCGGGCCTTTTTGCAGGGCGCGACTAAAGCCGCACAATTCTTGACCGAttt                                                                                                 <  5:340762/149‑1 (MQ=255)
                                                 cTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTtctc                                                                                             <  8:27784/149‑1 (MQ=255)
                                                 cTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTtctc                                                                                             >  7:69324/1‑149 (MQ=255)
                                                    taCATTCTGTTCGGTCCGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCg                                                                                            <  4:111994/146‑1 (MQ=255)
                                                    aaCATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGTCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTc                                                                                                            >  3:288695/1‑131 (MQ=255)
                                                    aaCATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGTCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTc                                                                                                            <  4:288695/131‑1 (MQ=255)
                                                    aaCATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCgg                                                                                           <  8:151004/148‑1 (MQ=255)
                                                    aaCATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCg                                                                                            <  1:292803/147‑1 (MQ=255)
                                                       aTTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCgg                                                                                           <  6:218215/145‑1 (MQ=255)
                                                       aTTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCgg                                                                                           >  5:218215/1‑145 (MQ=255)
                                                          cTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGc                                                                                     <  5:90397/148‑1 (MQ=255)
                                                          cTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGc                                                                                     <  8:56190/148‑1 (MQ=255)
                                                          cTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGc                                                                                     <  1:108599/148‑1 (MQ=255)
                                                             ttCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGa                                                                                  >  3:46392/1‑148 (MQ=255)
                                                                 gTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCg                                                                                            >  3:14160/1‑134 (MQ=255)
                                                                 gTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCg                                                                                            <  4:14160/134‑1 (MQ=255)
                                                                 gTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATg                                                                             <  4:46392/149‑1 (MQ=255)
                                                                 gTCTGGATGAGCGCATGGGCCTGTCCCATGCCCGCTTGCTGCCATAAGGCCGCCAGGATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCTACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATg                                                                             >  6:166511/1‑149 (MQ=255)
                                                                  tCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATgc                                                                            <  6:58392/149‑1 (MQ=255)
                                                                      gATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGc                                                                                     >  6:290491/1‑136 (MQ=255)
                                                                      gATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGc                                                                                     <  5:290491/136‑1 (MQ=255)
                                                                            cgcATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGa                                                                  <  4:186966/149‑1 (MQ=255)
                                                                             gcATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATgcgc                                                                          >  1:96126/1‑140 (MQ=255)
                                                                             gcATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATgcgc                                                                          <  2:96126/140‑1 (MQ=255)
                                                                              cATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGCGt                                                                <  4:157513/149‑1 (MQ=255)
                                                                               aTGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGtt                                                               <  2:93410/149‑1 (MQ=255)
                                                                                tGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCATGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGttt                                                              <  1:244095/149‑1 (MQ=255)
                                                                                    gcTGTCCCATGCCGGCTTGCTGCCGTAAGGCATCCAGCATTGGAAAAGGCCGGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCTCAATTCTTGAACGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTAt                                                           <  5:313965/147‑1 (MQ=255)
                                                                                      tGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATggc                                                        <  8:63808/149‑1 (MQ=255)
                                                                                          ccATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATggcggc                                                     >  3:255361/1‑148 (MQ=255)
                                                                                           cATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCtt                                                                                                          >  7:53608/1‑94 (MQ=255)
                                                                                           cATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCtt                                                                                                          <  8:53608/94‑1 (MQ=255)
                                                                                           cATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTAt                                                           <  2:185380/141‑1 (MQ=255)
                                                                                           cATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTAt                                                           >  1:185380/1‑141 (MQ=255)
                                                                                             tGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCTAcc                                                 >  3:170666/1‑149 (MQ=255)
                                                                                               ccGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCTACcgc                                               >  1:83157/1‑149 (MQ=255)
                                                                                                cGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTAt                                                           >  6:247179/1‑136 (MQ=255)
                                                                                                cGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTAt                                                           <  5:247179/136‑1 (MQ=255)
                                                                                                           cGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCTACCGCGCCGGGAGAAt                                    <  7:167842/148‑1 (MQ=255)
                                                                                                             tAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCTACCGCGCCGGGAGAATCaa                                 >  6:127938/1‑149 (MQ=255)
                                                                                                             tAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCTACCGCGCCGGGAGAATCaa                                 <  2:83157/149‑1 (MQ=255)
                                                                                                                 gccgccAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCTACCGCGCCGGGAGAATCAACAt                              >  8:262001/1‑148 (MQ=255)
                                                                                                                          aTTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCTACCGCGCCGGGAGAATCAACATGAGCGTCGaa                    >  5:218282/1‑149 (MQ=255)
                                                                                                                          aTTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCTACCGCGCCGGGAGAATCAACATGAGCGTCGa                     >  5:213392/1‑148 (MQ=255)
                                                                                                                           tgcaaaaaGGCCCGCACGACGCGGGCCTTTTCGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAGTGCGCGTCATCGAGTTTTATGGCGGCTCCCGCGCCGGGAGAATCCGCATGAGCGTCgcgc                   >  6:167741/5‑146 (MQ=255)
                                                                                                                            tGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCTACCGCGCCGGGAGAATCAACATGAGCGTCGAAAc                  >  1:325168/1‑149 (MQ=255)
                                                                                                                             ggAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCTACCGCGCCGGGAGAATCAACATGAGCGTCGAAAcc                 >  5:195904/1‑149 (MQ=255)
                                                                                                                              gAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCt                                                    >  3:282827/1‑113 (MQ=255)
                                                                                                                              gAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCt                                                    <  4:282827/113‑1 (MQ=255)
                                                                                                                              gAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCTACCGCGCCGGGAGAATCAACATGAGCGTCGAAAcc                 >  4:298844/1‑148 (MQ=255)
                                                                                                                              gAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCTACCGCGCCGGGAGAATCAACATGAGCGTCGAAACCt                <  7:78324/149‑1 (MQ=255)
                                                                                                                              aaaaaGGCCCGCACGACGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCTACCGCGCCGGGAGAATCAACATGAGCGTCGAAACCt                >  7:94558/2‑149 (MQ=255)
                                                                                                                               aaaaGGGCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCTACCGCGCCGGGAGAATCAACATGAGCGTCGAAACCtt               <  3:298844/149‑1 (MQ=255)
                                                                                                                                aaaGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCTACCGCGCCGGGAGAATCAACATGAGCGTCGAAACCTTc              <  6:195904/149‑1 (MQ=255)
                                                                                                                                    gCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCTACCGCGCCGGGAGAATCAACATGAGCGTCGAAACCTTCAccc          <  1:336422/149‑1 (MQ=255)
                                                                                                                                      ccGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCTACCGCGCCGGGAGAATCAACATGAGCGTCGAAACCTTCAcccc         >  5:198406/1‑148 (MQ=255)
                                                                                                                                             gtGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCTACCGCGCCGGGAGAATCAACATGAGCGTCGAAACCTTCACCCCTACCGGt  <  6:198406/148‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
GGTGATCGACAACCTGGTCAAGGGTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCTACCGCGCCGGGAGAATCAACATGAGCGTCGAAACCTTCACCCCTACCGGT  >  NC_002947/519535‑519823

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: