Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,693,834 C→G 6.3% R87R (CGC→CGG)  PP_1489 → CheW domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,693,8340CG 6.3% 200.4 / 4.3 63R87R (CGC→CGG) PP_1489CheW domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (24/35);  new base G (2/2);  total (26/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

GTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGCCTTATCTGGGCCCGGTGCGCCAGGATG  >  NC_002947/1693688‑1693982
                                                                                                                                                   |                                                                                                                                                    
gtgtTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCctg                                                                                                                                                   <  1:497163/150‑1 (MQ=255)
 tgtTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCctgg                                                                                                                                                  <  1:179859/150‑1 (MQ=255)
     gCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCctggtgct                                                                                                                                              <  2:333884/150‑1 (MQ=255)
     gCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCctggt                                                                                                                                                 >  2:1645658/1‑147 (MQ=255)
     gCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCctggt                                                                                                                                                 <  1:1645658/147‑1 (MQ=255)
      ccGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCAGGcgggtgctg                                                                                                                                             <  2:1432896/150‑1 (MQ=255)
        gcgTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCctggtgctggt                                                                                                                                           <  1:1152880/150‑1 (MQ=255)
              gCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTa                                                                                                                                     >  2:1531638/1‑150 (MQ=255)
                     aaGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAggc                                                                                                                              >  1:385035/1‑150 (MQ=255)
                      aGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAggcg                                                                                                                             <  1:1030656/150‑1 (MQ=255)
                               cccAGGCAGGTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGc                                                                                                                    >  1:349579/1‑150 (MQ=255)
                                  aGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTg                                                                                                                 <  1:497361/150‑1 (MQ=255)
                                    gCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGc                                                                                                               <  1:1225622/150‑1 (MQ=255)
                                    gCACCCGGCTGGGGCCCGGGCTTCCTCGCACATCGCGGCCAGCTGGTTTCGGTCATCGACCCTGCC‑CAGTCGGCTAAGTCTCCGTCGGCGTTGCAACGCACCAGCACCCGCCTGGTGCTGGTACACTACCGAGGCGGCCTGCAGCTTGGc                                                                                                               <  2:470848/150‑1 (MQ=255)
                                      aCCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCt                                                                                                             >  2:573487/1‑150 (MQ=255)
                                       cctgcctgGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTg                                                                                                            <  2:255387/150‑1 (MQ=255)
                                        ctgcctgGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCctg                                                                                                                                                   <  2:63812/110‑1 (MQ=255)
                                        ctgcctgGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCctg                                                                                                                                                   >  1:63812/1‑110 (MQ=255)
                                         tgcctgGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGGTTGAt                                                                                                          <  2:813578/150‑1 (MQ=255)
                                            ctgGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCt                                                                                                       >  1:1281970/1‑150 (MQ=255)
                                             tgGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTg                                                                                                      <  1:1442587/150‑1 (MQ=255)
                                             tgGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTg                                                                                                      <  2:385035/150‑1 (MQ=255)
                                               ggTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGa                                                                                                    <  1:17587/150‑1 (MQ=255)
                                                 tCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGc                                                                                                  <  1:225787/150‑1 (MQ=255)
                                                   gcgGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGCTCCAGGATCAg                                                                                                <  2:1402591/150‑1 (MQ=255)
                                                   gcgGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCAGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAg                                                                                                <  2:499604/150‑1 (MQ=255)
                                                      ggCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCctggtgctgg                                                                                                                                            <  2:1167001/103‑1 (MQ=255)
                                                      ggCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCctggtgctgg                                                                                                                                            >  1:1167001/1‑103 (MQ=255)
                                                       gCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCa                                                                                            >  2:7365/1‑150 (MQ=255)
                                                         aTTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCAcc                                                                                          <  1:573487/150‑1 (MQ=255)
                                                           tctcGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCAc                                                                                                                                       <  2:391111/103‑1 (MQ=255)
                                                           tctcGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCAc                                                                                                                                       >  1:391111/1‑103 (MQ=255)
                                                           tctcGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGa                                                                                        >  1:1395069/1‑150 (MQ=255)
                                                             tcGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCct                                                                                                                                                    <  2:351611/88‑1 (MQ=255)
                                                             tcGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCct                                                                                                                                                    >  1:351611/1‑88 (MQ=255)
                                                                      gTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGtt                                                                             <  1:7365/150‑1 (MQ=255)
                                                                           cAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCaa                                                                        >  2:937124/1‑150 (MQ=255)
                                                                                       gTCATCGACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCCTGCTGCTGGTGCACTACCGAGGCCGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAACCCCTGCGTTGCCAACCGGACGAGTTc                                                            >  1:760390/1‑150 (MQ=255)
                                                                                         cagcGACATTGTCGCCTTGAGGTTTGGC‑CAGCCGGCGGTGCAACGTCCCAGCACCCGCCTGGTCCTAGAGCACTACCGAGGCGGCCTGCAGTTTGGCCTGCGCCTGGAGCAGGCCACCGCAACCCTGCGTTGCCAACCGGTCGAGTTCCa                                                          <  2:172251/147‑1 (MQ=255)
                                                                                         cATCGACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCa                                                          >  2:1568284/1‑150 (MQ=255)
                                                                                            agACCTTGCCACACGGAGCTCTCGA‑CGGCCGGCGGGTCATCGTACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGcc                                                       <  2:760390/149‑1 (MQ=255)
                                                                                             gACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCtt                                                                                                                  >  2:1171929/1‑90 (MQ=255)
                                                                                             gACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCtt                                                                                                                  <  1:1171929/90‑1 (MQ=255)
                                                                                                   gccgccTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATtgg                                                <  2:257754/150‑1 (MQ=255)
                                                                                                   gccgccTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATggg                                                <  2:1395069/150‑1 (MQ=255)
                                                                                                   gccgccTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATggg                                                <  1:1174925/150‑1 (MQ=255)
                                                                                                      gccTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGTCCTATggggag                                             <  2:1138087/150‑3 (MQ=255)
                                                                                                        cTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGa                                           >  2:856214/1‑150 (MQ=255)
                                                                                                         ttGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGac                                          >  1:564856/1‑150 (MQ=255)
                                                                                                           gAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGacac                                        >  1:869619/1‑150 (MQ=255)
                                                                                                               tttGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGc                                    <  2:1390703/150‑1 (MQ=255)
                                                                                                                              ggTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGCTCCAGCCCTATGGGGTGGACACAGGCGATGCGCCTTATCt                      <  1:1568284/150‑1 (MQ=255)
                                                                                                                                  cAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGCCTTATCTGGGc                  <  1:1531638/150‑1 (MQ=255)
                                                                                                                                        aCCAGCACCAGGCGGGTGCTGGTGCGCTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGa                                                                >  1:489258/1‑98 (MQ=38)
                                                                                                                                        aCCAGCACCAGGCGGGTGCTGGTGCGCTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGa                                                                <  2:489258/98‑1 (MQ=38)
                                                                                                                                        aCCAGCACCAGGCGGGTGCTGGTGCGCTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGCCTTATCTGGGCCCGGTg            >  1:1161007/1‑150 (MQ=255)
                                                                                                                                          cAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGCCTTATCTGGGCCCGGTgcg          >  2:269875/1‑150 (MQ=255)
                                                                                                                                                ccGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGTCCTATGGGGTGGACACAGGCGATGCGCCTTATCTGGGCCCGGTGCGCCAGGa    >  1:420643/1‑150 (MQ=255)
                                                                                                                                                ccGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGCCTTATCTGGGCCCGGTGCGCCAGGa    <  2:1281970/150‑1 (MQ=255)
                                                                                                                                                ccGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCAACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGCCTTATCTGGGCCCGGTGCGCCAGGa    >  1:570830/1‑150 (MQ=255)
                                                                                                                                                 cGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCt                                                                                  >  1:1337113/1‑71 (MQ=255)
                                                                                                                                                 cGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCt                                                                                  <  2:1337113/71‑1 (MQ=255)
                                                                                                                                                  gCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGTCTTATCTGGGCCCGGTGCGCCAGGATg  >  1:689523/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGC‑CAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTGGAGCAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGCCTTATCTGGGCCCGGTGCGCCAGGATG  >  NC_002947/1693688‑1693982

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