Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,205,133 G→C 6.5% R31G (CGG→GGG)  PP_1056 ← iron‑chelator utilization protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,205,1330GC 6.5% 192.2 / 2.8 61R31G (CGG→GGG) PP_1056iron‑chelator utilization protein
Reads supporting (aligned to +/- strand):  ref base G (29/28);  new base C (2/2);  total (31/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 = 5.55e-01

ATGGTGGGCCGCGCGCCACCGTCCTTGCCCAGATTACGGGCCTCGATGGCCTGTTGTTGCTCCGGCGTTTCGGCGAACAGCAGCTTGATATGGTCGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCGCAGCACCTGCAGACGGCGCTGGCGGATCTCGTGGTTTACGCGGTGAATGGTGTCACTCATGGGGCGTCTCCCTGAAAGGTTCGGTGGGCCCGGAGGCGATGGCCTTGGCGGTAT  >  NC_002947/1204985‑1205277
                                                                                                                                                     |                                                                                                                                                
aTGGTGGGCCGCGCGCCACCGTCCTTGCCCAGATTACGGGCCTCGATGGCCTGTTGTTGCTCCGGCGTTTCGGCGAACAGCAGCTTGATATGGTCGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCg                                                                                                                                                  >  2:122582/1‑149 (MQ=255)
 tggtggGCCGCGCGCCACCGTCCTTGCCCAGATTACGGGCCTCGATGGCCTGTTGTTGCTCCGGGCGTTTCGGCGAACAGCAGCTTGATATGGTCGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCg                                                                                                                                                  >  1:286263/1‑149 (MQ=255)
   gtggGCCGCGCGCCACCGTCCTTGCCCAGATTACGGGCCTCGATGGCCTGTTGTTGCTCCGGCGTTTCGGCGAACAGCAGCTTGATATGGTCGTCGCTGCCCTTGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGgc                                                                                                                                                <  1:253023/148‑1 (MQ=255)
      ggCCGCGCGCCACCGTCCTTGCCCAGATTACGGACCTCGATGGCCTGTTGTTGCTCCGGCTTTTCGGCGAACAGCAGCTTGATATGGTCGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCAt                                                                                                                                            <  2:238618/149‑1 (MQ=255)
            gcgcCACCGTCCTTGCCCAGATTACGGGCCTCGATGGCCTGTTGTTGCTCCGGCGTTTCGGCGAACAGCAGCTTGATATGGTCGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCgg                                                                                                                                      <  3:7398/149‑1 (MQ=255)
                cACCGTCCTTGCCCAGATTACGGGCCTCGATGGCCTGTTGTTGCTCCGGCGTTTCGGCGAACAGCAGCTTGATATGGTCGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCAGGGt                                                                                                                                   >  5:313176/1‑148 (MQ=255)
                cACCGTCCTTGCCCAGATTACGGGCCTCGATGGCCTGTTGTTGCTCCGGCGTTTCGGCGAACAGCAGCTTGATATGGACGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGt                                                                                                                                   >  7:138959/1‑148 (MQ=255)
                 aCCGTCCTTGCCCAGATTACGGGCCTCGATGGCCTGTTGTTGCTCCGGCGTTTCGGCGAACAGCAGCTTGATATGGTCGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCggg                                                                                                                                     <  7:94884/145‑1 (MQ=255)
                 aCCGTCCTTGCCCAGATTACGGGCCTCGATGGCCTGTTGTTGCTCCGGCGTTTCGGCGAACAGCAGCTTGATATGGTCGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCggg                                                                                                                                     >  8:94884/1‑145 (MQ=255)
                          gCCCAGATTACGGGCCTCGATGGCCTGTTGTTGCTCCGGCGTTTCGGAGAACAGCAGCTTGATATGGTCGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGCTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGa                                                                                                                        >  3:30246/1‑149 (MQ=255)
                            ccAGATTACGGGCCTCGATGGCCTGTTGTTGCTCCGGCGTTTCGGCGAACAGCAGCTTGATATGGTCGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAg                                                                                                                                >  4:222818/1‑139 (MQ=255)
                            ccAGATTACGGGCCTCGATGGCCTGTTGTTGCTCCGGCGTTTCGGCGAACAGCAGCTTGATATGGTCGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAg                                                                                                                                <  3:222818/139‑1 (MQ=255)
                            ccAGATTACGGGCCTCGATGGCCTGTTGTTGCTCCGGCGTTTCGGCGAACAGCAGCTTGATATGGTCGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGAcc                                                                                                                      <  8:45577/149‑1 (MQ=255)
                                aTTACGGGCCTCGATGGCCTGTTGTTGCTCCGGCGTTTCGGCGAACAGCAGCTTGATATGGTCGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCgca                                                                                                                  <  3:280056/149‑1 (MQ=255)
                                       gCCTCGATGGCCTGTTGTTGCTCCGGCGTTTCGGCGAACAGCAGCTTGATATGGTCGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCACGGTGATCCGGCGCATGCGCGGCGTCAGCTCGGTGACCCGCCGCACCTg                                                                                                           >  4:223362/1‑149 (MQ=255)
                                             aTGGCCTGTTGTTGCTCCGGCGTTTCGGCGAACAGCAGCTTGATATGGTCGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCGCAGCACCTGCAGACg                                                                                                     >  6:86372/1‑149 (MQ=255)
                                                      tgttgCTCCGGCGTTTCGGCGAACAGCAGCTTGATATGGTCGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCGCAGCACCTGCAGACGGCGCTGGCg                                                                                            >  2:34747/1‑149 (MQ=255)
                                                       gttgCTCCGGCGTTTCGGCGAACAGCAGCTTGATATGGTCGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCGCAGCACCTGCAGACGGCGCTGGCg                                                                                            >  7:187664/1‑148 (MQ=255)
                                                       gttgCTCCGGCGTTTCGGCGAACAGCAGCTTGATATGGTCGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCGCAGCACCTGCAGACGGCGCTGGCg                                                                                            >  6:143856/1‑148 (MQ=255)
                                                            tCCGGCGTTTCGGCGAACAGCAGCTTGATATGGTCGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCGCAGCACCTGCAGACGGCGCTGGCGGAtctc                                                                                      >  4:344286/1‑149 (MQ=255)
                                                              cGGCGTTTCGGCGAACAGCAGCTTGATATGGTCGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCGCAGCACCTGCAGACGGCGCTGGCg                                                                                            >  4:59939/1‑141 (MQ=255)
                                                              cGGCGTTTCGGCGAACAGCAGCTTGATATGGTCGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCGCAGCACCTGCAGACGGCGCTGGCg                                                                                            <  3:59939/141‑1 (MQ=255)
                                                              cGGCGTTTCGGCGAACAGCAGCTTGATATGGTCGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCGCAGCACCTGCAGACGGCGCTGGCGGATCTCg                                                                                     <  5:290673/148‑1 (MQ=255)
                                                                 gCGTTTCGGCGAACAGCAGCTTGATATGGTCGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCGCAGCACCTGCAGACGGCGCTGGCGGATCTCGTgg                                                                                  >  1:358551/1‑149 (MQ=255)
                                                                  cGTTTCGGCGAACAGCAGCTTGATATGGTCGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCGCAGCACCTGCAGACGGCGCTGGCGGATCTCGTGGt                                                                                 <  8:187664/149‑1 (MQ=255)
                                                                       cGGCGAACAGCAGCTTGATATGGTCGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGATCCGCAGCACCTGCAGACGGCGCTGGCGGATCTCGTGGTTTAcg                                                                            >  6:210672/1‑149 (MQ=255)
                                                                          cGAACAGCAGCTTGATATGGTCGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCGCAGCACCTGCAGACGGCGCTGGCGGATCTCGTGGTTTACGCgg                                                                         <  2:358551/149‑1 (MQ=255)
                                                                           gAACAGCAGCTTGATATGGTCGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCGCAGCACCTGCAGACGGCGCTGGCGGATCTCGTGGTTTACGCGGt                                                                        <  5:143856/149‑1 (MQ=255)
                                                                            aaCAGCAGCTTGATATGGTCGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCGCAGCACCTGCAGACGGCGCTGGCGGATCTCGTGGTTTACGCGGt                                                                        >  6:355202/1‑148 (MQ=255)
                                                                            aaCAGCAGCTTGATATGGTCGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCGCAGCACCTGCAGACGGCGCTGGCGGATCTCGTGGTTTACGCGGTg                                                                       <  1:122582/149‑1 (MQ=255)
                                                                               agcagcTTGATATGGTCGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCGCAGCACCTGCAGACGGCGCTGGCGGATCTCGTGGTTTACGCGGTGAAt                                                                    <  1:34747/149‑1 (MQ=255)
                                                                                  agcTTGATATGGTCGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGATCCGCAGCACCTGCAGACGGCGCTGGCGGATCTCGTGGTTTACGCGGTGAATGgt                                                                 <  5:210672/149‑1 (MQ=255)
                                                                                       gATATGGTCGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCGCAGCACCTGCAGACGGCGCTGGCGGATCTCGTGGTTTACGCGGTGAATGGTGTCa                                                             <  5:86372/148‑1 (MQ=255)
                                                                                       gATATGGTCGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGACCCCGCGCATGCGCCGGATCac                                                                                                                                <  1:14750/81‑2 (MQ=38)
                                                                                       gATATGGTCGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGACCCCGCGCATGCGCCGGATCac                                                                                                                                <  1:20360/81‑2 (MQ=38)
                                                                                       gATATGGTCGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGACCCCGCGCATGCGCCGGATCac                                                                                                                                >  2:14750/1‑80 (MQ=38)
                                                                                       gATATGGTCGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGACCCCGCGCATGCGCCGGATCac                                                                                                                                >  2:20360/1‑80 (MQ=38)
                                                                                          atGGTCGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGAc                                                                                                                       >  1:353635/1‑87 (MQ=255)
                                                                                          atGGTCGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGAc                                                                                                                       <  2:353635/87‑1 (MQ=255)
                                                                                             gtcgtcGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCGCAGCACCTGCAGACGGCGCTGGCGGATCTCGTGGTTTACGCGGTGAATGGTGTCACTCATgg                                                      <  6:119765/149‑1 (MQ=255)
                                                                                                gtcgCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCGCAGCACCTGCAGACGGCGCTGGCGGATCTCGTGGTTTACGCGGTGAATGgtg                                                                <  3:279651/136‑1 (MQ=255)
                                                                                                gtcgCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCGCAGCACCTGCAGACGGCGCTGGCGGATCTCGTGGTTTACGCGGTGAATGgtg                                                                >  4:279651/1‑136 (MQ=255)
                                                                                                 tcgCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCGCAGCACCTGCAGACGGCGCTGGCGGATCTCGTGGTTTACGCGGTGAATGGTGTCACTCATGGGGCg                                                  >  3:64598/1‑149 (MQ=255)
                                                                                                          aTGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCGCAGCACCTGCAGACGGCGCTGGCGGATCTCGTGGTTTACGCGGTGAATGGGGTCACTCATGGGGCGTCTCCCTGa                                         <  6:293623/149‑1 (MQ=255)
                                                                                                             cTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCGCAGCACCTGCAGACGGCGCTGGCGGATCTCGTGGTTTACGCGGTGAATGGTGTCACTCATGGGGCGTCTCCCTGaaa                                       <  5:355202/148‑1 (MQ=255)
                                                                                                                gTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCGCAGCACCTGCAGACGGCGCTGGCgg                                                                                           >  6:230471/1‑93 (MQ=255)
                                                                                                                gTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCGCAGCACCTGCAGACGGCGCTGGCgg                                                                                           <  5:230471/93‑1 (MQ=255)
                                                                                                                 tAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCGCAGCACCTGCAGACGGCGCTGGCGGATCTCGTGGTTTACGCGGTGAATGGTGTCACTCATGGGGCGTCTCCCTGAAAGGTTc                                  >  1:41495/1‑149 (MQ=255)
                                                                                                                  aaaGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAg                                                                                                                                >  8:320933/1‑54 (MQ=255)
                                                                                                                  aaaGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAg                                                                                                                                <  7:320933/54‑1 (MQ=255)
                                                                                                                   aaGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCGCAGCACCTGCAGACGGCGCTGGCGGATCTCGTGGTTTACGCGGTGAATGGTGTCACTCATGGGGCGTCTCCCTGAAAGGTTCg                                 <  2:286263/148‑1 (MQ=255)
                                                                                                                      ccTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCGCAGCACCTGCAGACGGCGCTGGCGGATCTCGTGGTTTACGCGGTGAATGGTGTCACTCATGGGGCGTCTCCCTGAAAAGTTCGGTgg                             >  7:46051/1‑149 (MQ=255)
                                                                                                                          gCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCGCAGCACCTGCAGACGGCGCTGGCGGATCTCGTGGTTTACGCGGTGAATGGTGTCACTCATGGGGCGTCTCCCTGAAAGGTTCGGTGGGccc                         >  4:212612/1‑149 (MQ=255)
                                                                                                                          gCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCGCAGCACCTGCAGACGGCGCTGGCGGATCTCGTGGTTTACGCGGTGAATGGTGTCACTCATGGGGCGTCTCCCTGAAAGGTTCGGTGGGcc                          <  8:138959/148‑1 (MQ=255)
                                                                                                                             gCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCGCAGCACCTGCAGACGGCGCTGGCGGATCTCGTGGTTTACGCGGTGAATGGTGTCACTCATGGGGCGTCTCCCTGAAAGGTTCGGTGGGCCCg                        >  5:29644/1‑147 (MQ=255)
                                                                                                                             gCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGCTCAGCTCGGTGACCCCCAGCACCTGAACACGGCGCTGGCGGATCTCGTCCTATACGCGCTGAATGGTGTCACTCATGGGGCGTCTCCCTGAAAGGAAAGGTGGGCCCg                        >  1:300601/1‑147 (MQ=255)
                                                                                                                              cTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCGCAGCACCTGCAGACGGCGCTGGCGGATCTCGTGGTTTACGCGGTGAATGGTGTCACTCATGGGGCGTCTCCCTGAAAGGTTCGGTGGGCCCGGAg                     >  5:198232/1‑149 (MQ=255)
                                                                                                                                 ggCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCGCAGCACCTGCAGACGGCGCTGGCGGATCTCGGGGTTTACGCGGTGAATGGTGTCACTCATGGGGCGTCTCCCTGAAAGGTTCGGTGGGCCCGGAgg                    <  3:223362/147‑1 (MQ=255)
                                                                                                                                    aCCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCGCAGCACCTGCAGACGGCGCTGGCGGATCTCGTGGTTTACGCGGTGAATGGTGTCACTCATGGGGCGTCTCCCTGAAAGGTTCGGTGGGCCCGGAGGCGATg               <  2:41495/149‑1 (MQ=255)
                                                                                                                                       ccccAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCGCAGCACCTGCAGACGGCGCTGGCGGATCTCGTGGTTTACGCGGTGAATGGTGTCACTCATGGGGCGTCTCCCTGAAAGGTTCGGTGGGCCCGGAGGCGATGGcc            >  6:18991/1‑149 (MQ=255)
                                                                                                                                         ccAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCGCAGCACCTGCAGACGGCGCTGGCGGATCTCGTGGTTTACGCGGTGAATGGTGTCACTCATGGGGCGTCTCCCTGAAAGGtt                                   <  5:295641/124‑1 (MQ=255)
                                                                                                                                         ccAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCGCAGCACCTGCAGACGGCGCTGGCGGATCTCGTGGTTTACGCGGTGAATGGTGTCACTCATGGGGCGTCTCCCTGAAAGGtt                                   >  6:295641/1‑124 (MQ=255)
                                                                                                                                                 aTCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCGCAGCACCTGCAGACGGCGCTGGCGGATCTCGTGGTTTACGCGGTGAATGGTGTCACTCATGGGGCGTCTCCCTGAAAGGTTCGGTGGGCCCGGAGGCGATGGCCTTGGCGGTAt  >  3:78050/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
ATGGTGGGCCGCGCGCCACCGTCCTTGCCCAGATTACGGGCCTCGATGGCCTGTTGTTGCTCCGGCGTTTCGGCGAACAGCAGCTTGATATGGTCGTCGCTGCCCATGCTGGTAAAGCCTTGCAGCTCGGCACCCCCCAGGGTGATCCGGCGCATGCGCGGGGTCAGCTCGGTGACCCGCAGCACCTGCAGACGGCGCTGGCGGATCTCGTGGTTTACGCGGTGAATGGTGTCACTCATGGGGCGTCTCCCTGAAAGGTTCGGTGGGCCCGGAGGCGATGGCCTTGGCGGTAT  >  NC_002947/1204985‑1205277

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: