Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,760,424 T→C 6.4% intergenic (+82/‑229) folEB → / → PP_3325 GTP cyclohydrolase I/ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,760,4240TC 6.4% 177.6 / 2.9 63intergenic (+82/‑229)folEB/PP_3325GTP cyclohydrolase I/ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base T (35/23);  new base C (2/2);  total (37/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.48e-01

CCCCGCAGGCCGCTGCGCCTGACCAACAAGAATCATTTGCAACCGGTAATGAACGCTAACCACTAGCGTTACGTTATTACATCTAATAACATGCGCATCCCAACACATGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGCGCGGCATCGATGTTGGCATCGAGTGCAGCCTGCACAGCCTACACCGTGGAGCCTGATGGTGTACCCGTTGTGAT  >  NC_002947/3760284‑3760571
                                                                                                                                            |                                                                                                                                                   
ccccGCAGGCCGCTGCGCCTGACCAACAAGAATCATTTGCAACCGGTAATGAACGCTAACCACTAGCGTTACGTTATTACATCTAATAACATGCGCATCCCAACACATGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGc                                                                                                                                              >  3:89455/1‑148 (MQ=255)
 cccGCAGGCCGCTGCGCCTGACCAACAAGAATCATTTGCAACCGGTAATGAACGCTAACCACTAGCGTTACGTTATTACATCTAATAACATGCGCATCCCAACACATGCACATTTCCTCAGCAGGATTCATTTGTCAGGTCTGCCGc                                                                                                                                              >  6:157477/1‑147 (MQ=255)
  ccGCAGGCCGCTGCGCCTGACCAACAAGAATCATTTGCAACCGGTAATGAACGCTAACCACTAGCGTTACGTTATTACATCTAATAACATGCGCATCCCAACACATGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAg                                                                                                                                           >  7:82145/1‑149 (MQ=255)
            cTGCGCCTGACCAACAAGAATCATTTGCAACCGGTAATGAACGCTAACCACTAGCGTTACGTTATTACATCTAATAACATGCGCATCCCAACACATGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGct                                                                                                                                 >  4:284727/1‑149 (MQ=255)
               cgcCTGACCAACAAGAATCATTTGCAACCGGTAATGAACGCTAACCACTAGCGTTACGTTATTACATCTAATAACATGCGCATCCCAACACATGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGctct                                                                                                                               <  3:275334/148‑1 (MQ=255)
                gcCTGACCAACAAGAATCATTTGCAACCGGTAATGAACGCTAACCACTAGCGTTACGTTATTACATCTAATAACATGCGCATCCCAACACATGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTAc                                                                                                                             <  7:141434/149‑1 (MQ=255)
                 ccTGACCAACAAGAATCATTTGCAACCGGTAATGAACGCTAACCACTAGCGTTACGTTATTACATCTAATAACATGCGCATCCCAACACATGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACt                                                                                                                            <  2:222735/149‑1 (MQ=255)
                    gACCAACAAGAATCATTTGCAACCGGTAATGAACGCTAACCACTAGCGTTACGTTATTACATCTAATAACATGCGCATCCCAACACATGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCg                                                                                                                          >  7:111638/1‑148 (MQ=255)
                             gAATCATTTGCAACCGGTAATGAACGCTAACCACTAGCGTTACGTTATTACATCTAATAACATGCGCATCCCAACACATGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTc                                                                                                                >  3:254982/1‑149 (MQ=255)
                                           cGGTAATGAACGCTAACCACTAGCGTTACGTTATTACATCTAATAACATGCGCATCCCAACACATGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGAt                                                                                                  >  2:215442/1‑149 (MQ=255)
                                                aTGAACGCTAACCACTAGCGTTACGTTATTACATCTAATAACATGCGCATCCCAACACATGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATgcagc                                                                                             >  3:17203/1‑149 (MQ=255)
                                                aTGAACGCTAACCACTAGCGTTACGTTATTACATCTAATAACATGCGCATCCCAACACATGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATgcag                                                                                              <  4:89455/148‑1 (MQ=255)
                                                        tAACCACTAGCGTTACGTTATTACATCTAATAACATGCGCATCCCAACACATGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGa                                                                                                   >  3:49189/1‑135 (MQ=255)
                                                        tAACCACTAGCGTTACGTTATTACATCTAATAACATGCGCATCCCAACACATGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGa                                                                                                   <  4:49189/135‑1 (MQ=255)
                                                        tAACCACTAGCGTTACGTTATTACATCTAATAACATGCGCATCCCAACACATGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGccc                                                                                      >  3:151784/1‑148 (MQ=255)
                                                            cACTAGCGTTACGTTATTACATCTAATAACATGCGCATCCCAACACATGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAg                                                                                  >  4:50033/1‑148 (MQ=255)
                                                             aCTAGCGTTACGTTATTACATCTAATAACATGCGCATCCCAACACATGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATg                                                                                                 <  3:180773/132‑1 (MQ=255)
                                                             aCTAGCGTTACGTTATTACATCTAATAACATGCGCATCCCAACACATGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATg                                                                                                 >  4:180773/1‑132 (MQ=255)
                                                             aCTAGCGTTACGTTATTACATCTAATAACATGCGCATCCCAACACATGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCa                                                                                <  4:254982/149‑1 (MQ=255)
                                                             aCTAGCGTTACGTTATTACATCTAATAACATGCGCATCCCAACACATGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCa                                                                                <  5:305983/149‑1 (MQ=255)
                                                                 gCGTTACGTTATTACATCTAATAACATGCGCATCCCAACACATGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGgc                                                                            <  1:215442/149‑1 (MQ=255)
                                                                    ttacgttaTTACATCTAATAACATGCGCATCCCAACACATGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGcgcg                                                                         <  8:82145/149‑1 (MQ=255)
                                                                      acgttaTTACATCTAATAACATGCGCATCCCAACACATGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGCGCgg                                                                        >  7:49430/1‑148 (MQ=255)
                                                                        gttaTTACATCTAATAACATGCGCATCCCAACACATGAACATTTCGTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGCGCGGCAt                                                                     >  1:222429/1‑149 (MQ=255)
                                                                         ttattaCATCTAATAACATGCGCATCCCAACACATGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGCGCGGCATc                                                                    >  5:99880/1‑149 (MQ=255)
                                                                              aCATCTAATAACATGCGCATCCCAACACATGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGCGCGGCATCGATg                                                                <  3:284727/148‑1 (MQ=255)
                                                                                 tCTAATAACATGCGCATCCCAACACATGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATTTGTGATGCAGCACGGCCCGAAGCAGGGCGCGGCATCGATGTTgg                                                            >  2:221805/1‑149 (MQ=255)
                                                                                   taataaCATGCGCATCCCAACACATGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGCGCGGCATCGATGTTggaa                                                          >  7:70110/1‑147 (MQ=255)
                                                                                    aataaCATGCGCATCCCAACACATGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGCGCGGCATCGATGTTGGCAt                                                         >  3:299174/1‑149 (MQ=255)
                                                                                     ataaCATGCGCATCCCAACACATGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGCGCGGCATCGATGTTGGCATc                                                        <  4:17203/149‑1 (MQ=255)
                                                                                       aaCATGCGCATCCCAACACATGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGCGCGGCATCGATGTTGGCATc                                                        >  4:266345/1‑147 (MQ=255)
                                                                                         cATGCGCATCCCAACACATGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGCGCGGCATCGATGTTGGCATCGAGt                                                    <  4:299174/149‑1 (MQ=255)
                                                                                            gcgcATCCCAACACATGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGCGCGGCATCGATGTTGGCATCGAGTg                                                   >  6:212239/1‑147 (MQ=255)
                                                                                             cgcATCCCAACACATGAACATTTCCTCATCAGGATTCATTTGTCAGGGCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGCGCGGCATCGATGTTGGCATCGAGTGCa                                                 <  5:157477/148‑1 (MQ=255)
                                                                                                aTCCCAACACATGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTAATCTGTGATGCAGCACGGCACGAAGCAGGGCGCGGCATCGATGTTGGCATCGAGTGCAGCCt                                             >  2:308157/1‑149 (MQ=255)
                                                                                                   ccAACACATGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGCGCGGCATCGATGTTGGCATCGAGTGCAGCCTGca                                          <  8:6880/149‑1 (MQ=255)
                                                                                                    cAACACATGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGCGCGGCATCGATGTTg                                                             >  8:67280/1‑129 (MQ=255)
                                                                                                       cacaTGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGCGCGGCATCGATGTTGGCATCGAGTGCAGCCTGCACAGc                                      <  6:99880/149‑1 (MQ=255)
                                                                                                        acaTGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGCGCGGCATCGATGTTGGCATCGAGTGCAGCCTGCACAGc                                      >  8:271552/1‑148 (MQ=255)
                                                                                                          aTGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGCGCGGCATCGATGTTGGc                                                           >  7:125514/1‑125 (MQ=255)
                                                                                                          aTGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGCGCGGCATCGATGTTGGc                                                           <  8:125514/125‑1 (MQ=255)
                                                                                                                aTTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGCGCGGCATCGATGTTGGCATCGAGTGCAGCCTGCACAGCCTACACCGt                             >  8:95445/1‑149 (MQ=255)
                                                                                                                 tttCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGCGCGGCATCGATGTTGGCATCGAGTGCAGCCTGCACAGCCTACACCGTg                            <  8:226895/149‑1 (MQ=255)
                                                                                                                     cTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGCGCGGCATCGATGTTGGCATCGAGTGCAGCCTGCACAGCCTACACCGTGGAGc                        >  6:111961/1‑149 (MQ=255)
                                                                                                                     cTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGCGCGGCATCGATGTTGGCATCGAGTGCAGCCTGCACAGCCTACACCGTGGAGc                        <  1:308157/149‑1 (MQ=255)
                                                                                                                       catcaGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGCGCGGCATCGATGTTGGCATCGAGTGCAGCCTGCACAGCCTACACCGTGGAGCCt                      >  2:163721/1‑149 (MQ=255)
                                                                                                                        atcaGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGCGCGGCATCGATGTTGGCATCGAGTGCAGCCTGcaca                                        <  5:75660/130‑1 (MQ=255)
                                                                                                                        atcaGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGCGCGGCATCGATGTTGGCATCGAGTGCAGCCTGcaca                                        >  6:75660/1‑130 (MQ=255)
                                                                                                                          caGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGCGCGGCATCGATGTTGGCATCGAGt                                                    <  1:7337/116‑1 (MQ=255)
                                                                                                                          caGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGCGCGGCATCGATGTTGGCATCGAGt                                                    >  2:7337/1‑116 (MQ=255)
                                                                                                                          caGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGCGCGGCATCGATGTTGGCATCGAGTGCAGCCTGCACAGCCTACACCGTGGAGCCTGAt                   >  8:100704/1‑149 (MQ=255)
                                                                                                                           aGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGCGCGGCATCGATGTTGGCATCGAGTGCAGCCTGCACAGCCTACACCGTGGAGCCTGATg                  <  3:50033/149‑1 (MQ=255)
                                                                                                                            ggATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGCGCGGCATCGATGTTGGCATCGAGTGCAGCCTGCACAGCCTACACCGTGGAGCCTGATgg                 >  1:125091/1‑149 (MQ=255)
                                                                                                                              aTTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGCGCGGCATCGATGTTGGCATCGAGTGCAGCCTGCACAGCCTACACCGTGGAGCCTGATGgtg               >  6:313479/1‑149 (MQ=255)
                                                                                                                                 cttgcgTCAGACCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGCGCGGCATCGATGTTGGCATCGAGTGCAGCCTGCACAGCCTacac                                <  2:93061/124‑1 (MQ=255)
                                                                                                                                 cttgcgTCAGACCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGCGCGGCATCGATGTTGGCATCGAGTGCAGCCTGCACAGCCTacac                                <  2:96835/124‑1 (MQ=255)
                                                                                                                                 cttgcgTCAGACCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGCGCGGCATCGATGTTGGCATCGAGTGCAGCCTGCACAGCCTacac                                >  1:96835/6‑129 (MQ=255)
                                                                                                                                 cttgcgTCAGACCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGCGCGGCATCGATGTTGGCATCGAGTGCAGCCTGCACAGCCTacac                                >  1:93061/6‑129 (MQ=255)
                                                                                                                                   tttGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGCGCGGCATCGATGTTGGCATCGAGTGCAGCCTGCACAGCCTACACCGTGGAGCCCGATGGTGTAccc          >  2:293792/1‑149 (MQ=255)
                                                                                                                                       tCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGa                                                                                                   <  6:46831/56‑1 (MQ=255)
                                                                                                                                       tCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGa                                                                                                   >  5:46831/1‑56 (MQ=255)
                                                                                                                                       tCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGCGCGGCATCGATGTTGGCATCGAGTGCAGCCTGCACAGCCTACACCGTGGAGCCTGATGGTGTACCCGtt       <  3:300582/148‑1 (MQ=255)
                                                                                                                                           gTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGCGCGGCATCGATGTTGGCCTCGAGTGCAGCCTGCACAGCCTACACCGTGGAGCCTGATGGTGTACCCGTTGTGAt  >  4:276150/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
CCCCGCAGGCCGCTGCGCCTGACCAACAAGAATCATTTGCAACCGGTAATGAACGCTAACCACTAGCGTTACGTTATTACATCTAATAACATGCGCATCCCAACACATGAACATTTCCTCATCAGGATTCATTTGTCAGGTCTGACGCAAGGCGGGCCGCTCTACTCGCAAGGGTTTCTTGTCATCTGTGATGCAGCACGGCCCGAAGCAGGGCGCGGCATCGATGTTGGCATCGAGTGCAGCCTGCACAGCCTACACCGTGGAGCCTGATGGTGTACCCGTTGTGAT  >  NC_002947/3760284‑3760571

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: