Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,361,428 T→C 6.6% H164H (CAT→CAC)  phoP → two‑component system DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,361,4280TC 6.6% 174.8 / 2.7 60H164H (CAT→CAC) phoPtwo‑component system DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (28/28);  new base C (2/2);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.69e-01

TGTTGCGCCGCTCCAGTGGTTTCACCCAGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGCGAGCGCGGCTTCAAACCCATCGACACGGT  >  NC_002947/1361280‑1361574
                                                                                                                                                    |                                                                                                                                                  
tGTTGCGCCGCTCCAGTGGTTTCACCCAGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACATCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGTCCCTGACCTCCTACGAATTCCGCATCCTCTAATCACTAATGCGACAt                                                                                                                                                    >  8:348766/1‑149 (MQ=255)
     cgccgcTCCAGTGGTTTCACCCAGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCAcc                                                                                                                                                <  8:4585/148‑1 (MQ=255)
       ccgcTCCAGTGGTTTCACCCAGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACcag                                                                                                                                              >  8:92452/1‑148 (MQ=255)
        cgcTCCAGTGGTTTCACCCAGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACcagc                                                                                                                                             <  4:109579/148‑1 (MQ=255)
        cgcTCCAGTGGTTTCACCCAGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACcagc                                                                                                                                             >  3:109579/1‑148 (MQ=255)
            ccAGTGGCTTCACCCAGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAggtg                                                                                                                                        <  4:68977/149‑1 (MQ=255)
                  gTTTCACCCAGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGc                                                                                                                                   >  7:196533/1‑148 (MQ=255)
                   tttCACCCAGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCACCACCTGCTGGTGATGGcgc                                                                                                                                 <  1:196478/149‑3 (MQ=255)
                     tCACCCAGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCATCCGCAAGCAGGCCACCGTGGACGAAGAAGCCTTGGCCCTGTCCGCCTAAGAATACCGCATCCTCTATTTCCTCATGCGCCATCACCTGCAGGTGGTGGCCaa                                                                                                                                <  6:225949/148‑1 (MQ=255)
                       aCCCAGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGa                                                                                                                             >  6:202275/1‑149 (MQ=255)
                       aCCCAGTCGACCATCCCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCACCACCTGCTGGTGATGGcgcatga                                                                                                                             <  8:335607/149‑7 (MQ=255)
                         ccAGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATc                                                                                                                           <  7:92452/149‑1 (MQ=255)
                         ccAGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATc                                                                                                                           <  5:252182/149‑1 (MQ=255)
                           aGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGc                                                                                                                         >  2:189671/1‑149 (MQ=255)
                           aGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGc                                                                                                                         >  2:110887/1‑149 (MQ=255)
                            gTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACTAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATTTTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGcc                                                                                                                        <  1:87488/149‑1 (MQ=255)
                               gACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGa                                                                                                                     <  6:101513/149‑1 (MQ=255)
                                 ccATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATg                                                                                                                   <  1:189671/149‑1 (MQ=255)
                                 ccATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATg                                                                                                                   >  3:257397/1‑149 (MQ=255)
                                   aTCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCTAGCAGGCCACCCGCGACGAACAAGCCTTGGCCCTGCCCGCCTACGAATACCGCAGCCTCGCATACCTCCTGCGCCCTCACCCGCAGGTGGTGGCCAAGGATCGCCGGATg                                                                                                                   >  2:125589/1‑147 (MQ=255)
                                   aTCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGa                                                                                                                 <  1:110887/149‑1 (MQ=255)
                                   aTCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGa                                                                                                                 <  6:186335/149‑1 (MQ=255)
                                      gccgcTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCa                                                                                                              <  7:343842/149‑1 (MQ=255)
                                      gccgcTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCa                                                                                                              <  7:348766/149‑1 (MQ=255)
                                        cgcTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGt                                                                                                            >  8:321223/1‑149 (MQ=255)
                                         gcTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGtt                                                                                                           <  6:188133/149‑1 (MQ=255)
                                         gcTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGtt                                                                                                           <  8:303599/149‑1 (MQ=255)
                                         gcTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGtt                                                                                                           <  5:202275/149‑1 (MQ=255)
                                             gCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTa                                                                                                        <  1:12223/148‑1 (MQ=255)
                                               cccTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTAcc                                                                                                      >  4:253944/1‑148 (MQ=255)
                                                ccTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCa                                                                                                    <  2:24158/149‑1 (MQ=255)
                                                           gACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGaa                                                                                           >  6:178101/1‑147 (MQ=255)
                                                               tCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGCCCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGCAGCACTTGAACCCAGGAGACGAAGAgcgc                                                                                     >  3:231957/1‑149 (MQ=255)
                                                                 aaCCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCATGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAaga                                                                                         >  8:338945/1‑143 (MQ=255)
                                                                 aaCCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAaga                                                                                         <  7:338945/143‑1 (MQ=255)
                                                                     gCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGAccc                                                                                >  2:251624/1‑148 (MQ=255)
                                                                              ccACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATc                                                                      <  2:186459/149‑1 (MQ=255)
                                                                                 cccTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGaa                                                                   >  8:355260/1‑149 (MQ=255)
                                                                                 cccTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGaa                                                                   <  8:196533/149‑1 (MQ=255)
                                                                                      gacgaACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCt                                                              >  3:171285/1‑149 (MQ=255)
                                                                                               gCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGgccg                                                     >  1:257902/1‑149 (MQ=255)
                                                                                                    ggCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCt                                                  <  3:165640/147‑1 (MQ=255)
                                                                                                       ccTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGcgcc                                             >  7:193786/1‑149 (MQ=255)
                                                                                                       ccTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGcgcc                                             >  8:192463/1‑149 (MQ=255)
                                                                                                         tGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGcgccg                                            >  3:112722/1‑148 (MQ=255)
                                                                                                           aCCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGcgccgc                                           <  1:67838/147‑1 (MQ=255)
                                                                                                           aCCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCaa                                         <  8:302368/149‑1 (MQ=255)
                                                                                                               ccTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCtt                                     >  2:219441/1‑149 (MQ=255)
                                                                                                                 tACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTg                                    >  3:276635/1‑148 (MQ=255)
                                                                                                                    gAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTAGAAgg                                >  5:317326/1‑149 (MQ=255)
                                                                                                                     aaTACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGc                               <  8:193786/149‑1 (MQ=255)
                                                                                                                       tACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGc                               >  7:243505/1‑147 (MQ=255)
                                                                                                                         ccGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGCGAgc                           >  8:266687/1‑149 (MQ=255)
                                                                                                                                         cTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGc                               <  8:315618/129‑1 (MQ=255)
                                                                                                                                         cTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGc                               >  7:315618/1‑129 (MQ=255)
                                                                                                                                           ccttggcCACCACCTGCTGGTGATGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGCGAGCGCGGCTTCAAACCCATc          >  3:332365/6‑148 (MQ=255)
                                                                                                                                             ttggcCACCACCTGCTGGTGATGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGCGAGCGCGGCTTCAAACCCATCGa        >  3:249929/4‑148 (MQ=255)
                                                                                                                                               cgcCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCt                                                              <  6:321558/92‑1 (MQ=255)
                                                                                                                                               cgcCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCt                                                              >  5:321558/1‑92 (MQ=255)
                                                                                                                                                 ccATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGCGAGCGCGGCTTCAAACCCATCGACACgg   <  4:112722/149‑1 (MQ=255)
                                                                                                                                                   aTCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGCGAGCGCGGCTTCAAACCCATCGACACGGt  >  5:109048/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TGTTGCGCCGCTCCAGTGGTTTCACCCAGTCGACCATCGCCGCTGGCCCCTTGGTGCTTGACCTCAACCGCAAGCAGGCCACCCTCGACGAACAAGCCTTGGCCCTGACCGCCTACGAATACCGCATCCTCGAATACCTCATGCGCCATCACCAGCAGGTGGTGGCCAAGGATCGCCTGATGGAGCAGTTGTACCCAGGAGACGAAGAGCGCGACCCCAATGTCATCGAAGTGCTGGTTGGCCGCCTGCGCCGCAAGCTTGAAGGCGAGCGCGGCTTCAAACCCATCGACACGGT  >  NC_002947/1361280‑1361574

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