Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,639,173 C→T 6.5% A94T (GCG→ACG)  czcR‑III ← response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,639,1730CT 6.5% 177.7 / 2.9 61A94T (GCG→ACG) czcR‑IIIresponse regulator
Reads supporting (aligned to +/- strand):  ref base C (31/26);  new base T (2/2);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.18e-01

GTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGCATCATGATGCGTGTGGCCGATTCGGCACGTAGCCGCTGCAGCACGGTCCAGCCATCGATTTCGGGCAGGTTGACGTCGAGGATGACCAGGTCA  >  NC_002947/1639036‑1639318
                                                                                                                                         |                                                                                                                                                 
gTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGAGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAg                                                                                                                                       <  1:1182271/150‑1 (MQ=255)
           ggTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTAccc                                                                                                                            <  2:122331/150‑1 (MQ=255)
            gTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGACTTTTACACg                                                                                                                           >  2:161602/1‑150 (MQ=255)
             tCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCgg                                                                                                                          >  1:1198494/1‑150 (MQ=255)
               cAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGTTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTc                                                                                                                        <  2:170181/150‑1 (MQ=255)
               cAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAGCGTCTTGAGCGGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTc                                                                                                                        <  1:161602/150‑1 (MQ=255)
                           aGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCTAAg                                                                                                                 >  1:260668/1‑145 (MQ=255)
                           aGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCTAAg                                                                                                                 <  2:260668/145‑1 (MQ=255)
                             cGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGt                                                                                                          <  1:501070/150‑1 (MQ=255)
                                 gcgcAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGc                                                                                                      >  1:1207141/1‑150 (MQ=255)
                                   gcAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCgg                                                                                                    >  2:1160056/1‑150 (MQ=255)
                                       ggTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCACCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAg                                                                                                >  2:492515/1‑150 (MQ=255)
                                           cTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGcatc                                                                                            <  2:517767/150‑1 (MQ=255)
                                            tGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGGTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGcatca                                                                                           >  2:642116/1‑150 (MQ=255)
                                             gggCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCgcgc                                                                                                                                                 >  1:35524/1‑95 (MQ=255)
                                             gggCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCgcgc                                                                                                                                                 <  2:35524/95‑1 (MQ=255)
                                             gggCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGcatcat                                                                                          <  2:1207141/150‑1 (MQ=255)
                                              ggCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCTGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGCATCatg                                                                                         <  2:1142191/150‑1 (MQ=255)
                                                     gTTGCTGGTCGTTGCGGCGTAGCAGTCTTCGGATGCTCGCCAGTAGCTCATGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGCATCatgatgagca                                                                                  <  1:685720/150‑5 (MQ=255)
                                                      ttGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGCATCATGATGCgtgt                                                                                 >  2:104859/1‑150 (MQ=255)
                                                      ttGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGCATCATGATGCgtgt                                                                                 <  1:828343/150‑1 (MQ=255)
                                                         cTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGCATCATGATGCGTGTGGc                                                                              <  2:1198494/150‑1 (MQ=255)
                                                            gTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGCATCATGATGCGTGTGGCCgg                                                                           >  2:85699/1‑149 (MQ=255)
                                                            gTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGCATCATGATGCGTGTGGCCGa                                                                           >  1:107946/1‑150 (MQ=255)
                                                                     gcgTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGCATCATGATGCGTGTGGCCGATTCGGCACg                                                                  >  2:886791/1‑150 (MQ=255)
                                                                        tAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGTGACCGTgg                                                                                                        <  2:476320/109‑1 (MQ=255)
                                                                        tAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGTGACCGTgg                                                                                                        >  1:476320/1‑109 (MQ=255)
                                                                         aggaGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGCATCATGATGCGTGTGGCCGATTCGGCACGTAGc                                                              <  2:603714/147‑1 (MQ=255)
                                                                         aGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGt                                                                                                   >  1:773049/1‑113 (MQ=255)
                                                                         aGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGt                                                                                                   <  2:773049/113‑1 (MQ=255)
                                                                          gCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGCATCATGATGCGTGTGGCCGATTCGGCACGTAGcc                                                             >  2:698977/1‑150 (MQ=255)
                                                                          gCAGGCTGCGGATGCGCACCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGCATCATGATGCGTGTGGCCGATTCGGCACGTAGcc                                                             <  1:492515/150‑1 (MQ=255)
                                                                                   ggATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGCATCATGATGCGTGGGGCCGATTAGGCACGTAGCCGCTgcagca                                                    >  1:556071/1‑150 (MQ=255)
                                                                                      tGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGCATCATGATGCGTGTGGCCGATTCGGCACGTAGCCGCTgcagca                                                    >  1:733012/1‑147 (MQ=255)
                                                                                      tGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGCATCATGATGCGTGTGGCCGATTCGGCACGTAGCCGCTgcagca                                                    <  2:733012/147‑1 (MQ=255)
                                                                                      tGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGCATCATGATGCGTGTGGCCGATTCGGCACGTAGCCGCTGCAGCACgg                                                 >  1:950719/1‑150 (MQ=255)
                                                                                                       gggAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGCATCATGATGCGTGTGGCCGATTCGGCACGTAGCCGCTGCAGCACGGTCCCGCCATCGATTTCg                                <  1:642116/150‑1 (MQ=255)
                                                                                                        ggAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGCATCATGATGCGTGTGGCCGATTCGGCACGTAGCCGCTGCAGCACGGTCCAGCCATCGATTTCgg                               >  2:855549/1‑150 (MQ=255)
                                                                                                         gAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGCATCATGATGCGTGTGGCCGATTCGGCACGTAGCCGCTGCAGCACGGTCCAGCCATCGATTTCggg                              <  1:886791/150‑1 (MQ=255)
                                                                                                             tCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGCATCATGATGCGTGTGGCCGATTCGGCACGTAGCCGCTGCAGCACGGTCCAGCCATCGATTTCGGGCAgg                          >  2:345199/1‑150 (MQ=255)
                                                                                                             tCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGCATCATGATGCGTGTGGCCGATTCGGCACGTAGCCGCTGCAGCACGGTCCAGCCATCGATTTCGGGCAgg                          >  1:71649/1‑150 (MQ=255)
                                                                                                             tCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGCATCATGATGCGTGTGGCCGATTCGGCACGTAGCCGCTGCAGCACGGTCCAGCCATCGATTTCGGGCAgg                          >  1:633350/1‑150 (MQ=255)
                                                                                                               aaaGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGCATCATGATGCGTGTGGCCGATTCGGCACGTAGCCGCTGCAGCACGGTCCAGCCATCGATTTCGGGCAGGtt                        >  1:1256689/1‑150 (MQ=255)
                                                                                                                  gggCTTGAGCAGGTAGTCATCCGTGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGCATCATGATGCGTGTGGCCGATTCGGCACGTAGcc                                                             <  2:310987/110‑1 (MQ=255)
                                                                                                                  gggCTTGAGCAGGTAGTCATCCGTGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGCATCATGATGCGTGTGGCCGATTCGGCACGTAGcc                                                             >  1:310822/1‑110 (MQ=255)
                                                                                                                  gggCTTGAGCAGGTAGTCATCCGTGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGCATCATGATGCGTGTGGCCGATTCGGCACGTAGcc                                                             >  1:310987/1‑110 (MQ=255)
                                                                                                                  gggCTTGAGCAGGTAGTCATCCGTGCCCAGGTCCAGGCCTTTTACCCGGGCGGCCAAGCGACCGTGGGCGGTCAGCATCATGATGCGTGTGGCCGATTCGGCACGTAGcc                                                             <  2:310822/110‑1 (MQ=255)
                                                                                                                  gggCGTGGGCATTTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGCATCATGATGCGTGTGGCCGATTCGGCACGTAGCCGCTGCAGCACGGTCCAGCCATCGATTTCGGGCAGGTTGAc                     <  2:71649/150‑1 (MQ=255)
                                                                                                                    gCTTGAGCAGGTAGTCATCCGCACCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGc                                                                                                                <  1:269088/57‑1 (MQ=255)
                                                                                                                    gCTTGAGCAGGTAGTCATCCGCACCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGc                                                                                                                >  2:269088/1‑57 (MQ=255)
                                                                                                                           cAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGCATCATGATGCGTGTGGCCGAtt                                                                         >  1:522899/1‑89 (MQ=255)
                                                                                                                           cAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGCATCATGATGCGTGTGGCCGAtt                                                                         <  2:522899/89‑1 (MQ=255)
                                                                                                                             ggTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGCATCATGATGCGTGTGGCCGATTCGGCACGTAGCCGCTGCAGCACGGTCCAGCCATCGATTTCGGGCAGGTTGACGTCGAGGATGa          <  2:1256689/150‑1 (MQ=255)
                                                                                                                              gTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGAcc                                                                                                            >  2:757408/1‑51 (MQ=255)
                                                                                                                              gTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGAcc                                                                                                            <  1:757408/51‑1 (MQ=255)
                                                                                                                              gTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGCATCATGATGCGTGTGGCCGATTCGGCACGTAGCCGCTGCAGCACGGTCCAGCCATCGATTTc                                 >  1:911629/1‑126 (MQ=255)
                                                                                                                              gTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGCATCATGATGCGTGTGGCCGATTCGGCACGTAGCCGCTGCAGCACGGTCCAGCCATCGATTTc                                 <  2:911629/126‑1 (MQ=255)
                                                                                                                              gTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGCATCATGATGCGTGTGGCCGATTCGGCACGTAGCCGCTGCAGCACGGTCCAGCCATCGATTTCGGGCAGGTTGACGTCGAGGATGAc         >  2:611661/1‑150 (MQ=255)
                                                                                                                                  tCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGCATCATGATGCGTGTGGCCGATTCGGCACGTAGCCGCTGCAGCACGGTCCAGCCATCGATTTCGGGCAGGTTGACGTCGAGGATGACCAgg     >  2:884752/1‑150 (MQ=255)
                                                                                                                                    aTCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGCATCATGATGCGTGTGGCCGATTCGGCACGTAGCCGCTGCAGCACGGTCCAGCCATCGATTTCGGGCAGGTTGACGTCGAGGATGACCAGGTc   >  2:607236/1‑150 (MQ=255)
                                                                                                                                    aTCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGCATCATGATGCGTGTGGCCGATTCGGCACGTAGCCGCTGCAGCACGGTCCAGCCATCGATTTCGGGCAGGTTGACGTCGAGGATGACCAGGTc   >  1:106537/1‑150 (MQ=255)
                                                                                                                                    aTCCGCGCCCAGGGCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGCATCATGATGCGTGTGGCCGATTCGGCACGTAGCCGCTGCAGCACGGTCCAGCCATCGATTTCGGGCAGGTTGACGTCGAGGATGACCAGGTc   >  2:261956/1‑150 (MQ=255)
                                                                                                                                     tCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGCATCATGATGCGTGTGGCCGATTCGGCACGTAGCCGCTGCAGCACGGTCCAGCCATCGATTTCGGGCAGGTTGACGTCGAGGATGACCAGGTCa  <  1:1160056/150‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                 
GTGGCGGCCGGGGTCCAGTTCCAGGTCAGCGACGCGCAGGGTGCTGGGCTGCAGTTGCTGGTCGTTGCGGCGTAGCAGGCTGCGGATGCGCGCCAGTAGCTCAGGGAACTCAAAGGGCTTGAGCAGGTAGTCATCCGCGCCCAGGTCCAGGCCTTTTACCCGGTCGGCCAAGCGACCGTGGGCGGTCAGCATCATGATGCGTGTGGCCGATTCGGCACGTAGCCGCTGCAGCACGGTCCAGCCATCGATTTCGGGCAGGTTGACGTCGAGGATGACCAGGTCA  >  NC_002947/1639036‑1639318

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