Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,831,323 T→C 6.8% K84R (AAA→AGA)  PP_2484 ← ArsR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,831,3230TC 6.8% 170.2 / 2.4 59K84R (AAA→AGA) PP_2484ArsR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (25/30);  new base C (2/2);  total (27/32)
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.50e-01

CGAAGCTTCGGATGATCTGATCGTTATTTGGTCAGATTGCTGTTTTTGACATCGGCAAAGACTGAAATCAGTCAGAAATACAAGGGTTGCAGCCCTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTGCGAGATAAAGCGATGTCGTTGACTGCGAGAGGCCCGTGCGCTCTTGGATGATACTGACGCACACACCTACCGTTTCTGGGTCGACCT  >  NC_002947/2831174‑2831470
                                                                                                                                                     |                                                                                                                                                   
cGAAGCTTCGGATGATCTGATCGTTATTTGGTCAGATTGCTGTTTTTGACATCGGCAAAGACTGAAATCAGTCAGAAATACAAGGGTTGCAGCCCTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGtt                                                                                                                                                     <  2:527768/150‑1 (MQ=255)
 gAAGCTTCGGATGATCTGATCGTTATTTGGTCAGATTGCTGTTTTTGACATCGGCAAAGACTGAAATCAGTCAGAAATACAAGGGTTGCAGCCCTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGttt                                                                                                                                                    >  1:955402/1‑150 (MQ=255)
     cTTCGGATGATCTGATCGTTATTTGGTCAGATTGCTGTTTTTGACATCGGCAAAGACTGAAATCAGTCAGAAATACAAGGGTTGCAGCCCTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTgtag                                                                                                                                                <  2:1103084/150‑1 (MQ=255)
     cTTCGGATGATCTGATCGTTATTTGGTCAGATTGCTGTTTTTGACATCGGCAAAGACTGAAATCAGTCAGAAATACAAGGGTTGCAGCCCTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTgtag                                                                                                                                                <  2:955402/150‑1 (MQ=255)
          gatgatCTGATCGTTATTTGGTCAGATTGCTGTTTTTGACATCGGCAAAGACTGAAATCAGTCAGAAATACAAGGGTTGCAGCCCTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGtagtag                                                                                                                                             <  2:2058/148‑1 (MQ=255)
          gatgatCTGATCGTTATTTGGTCAGATTGCTGTTTTTGACATCGGCAAAGACTGAAATCAGTCAGAAATACAAGGGTTGCAGCCCTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGtagtag                                                                                                                                             >  1:2058/1‑148 (MQ=255)
          gatgatCTGATCGTTATTTGGTCAGATTGCTGTTTTTGACATCGGCAAAGACTGAAATCAGTCAGAAATACAAGGGTTGCAGCCCTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGt                                                                                                                                           >  1:439451/1‑150 (MQ=255)
                gtggtcGTTGTTTGGTCAGATTGCTGTTTTTGACATCGGCAAAGACTGAAATCAGTCAGAAATACAAGGGTTGCAGCCCTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAggg                                                                                                                                     <  2:439451/146‑1 (MQ=255)
                     cGTTATTTGGTCAGATTGCTGTTTTTGACATCGGCAAAGACTGAAATCAGTCAGAAATACAAGGGTTGCAGCCCTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGa                                                                                                                                <  1:649932/150‑1 (MQ=255)
                                   aTTGCTGTTTTTGACATCGGCAAAGACTGAAATCAGTCAGAAATACAAGGGTTGCAGCCCTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCAAGGGGCCGATCCGCTGAAATGtt                                                                                                                  >  2:1370835/1‑149 (MQ=255)
                                        tGTTTTTGACATCGGCAAAGACTGAAATCAGTCAGAAATACAAGGGTTGCAGCCCTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTCTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCaa                                                                                                             >  2:966043/1‑150 (MQ=255)
                                        tGTTTTTGACATCGGCAAAGACTGAAATCAGTCAGAAATACAAGGGTTGCAGCCCTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTCTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCaa                                                                                                             >  2:965995/1‑150 (MQ=255)
                                                  aTCGGCAAAGACTGAAATCAGTCAGAAATACAAGGGTTGCAGCCCTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAAtt                                                                                                           >  2:604652/1‑142 (MQ=255)
                                                  aTCGGCAAAGACTGAAATCAGTCAGAAATACAAGGGTTGCAGCCCTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAAtt                                                                                                           <  1:604652/142‑1 (MQ=255)
                                                    cGGCAAAGACTGAAATCAGTCAGAAATACAAGGGTTGCAGCCCTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTg                                                                                                 <  2:141170/150‑1 (MQ=255)
                                                        aaaGACTGAAATCAGTCAGAAATACAAGGGTTGCAGCCCTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTa                                                                                               <  1:476258/148‑1 (MQ=255)
                                                        aaaGACTGAAATCAGTCAGAAATACAAGGGTTGCAGCCCTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTa                                                                                               >  2:476258/1‑148 (MQ=255)
                                                            aCTGAAATCAGTCAGAAATACAAGGGTTGCAGCCCTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTg                                                                                         >  2:619056/1‑150 (MQ=255)
                                                             cTGAAATCAGTCAGAAATACAAGGGTTGCAGCCCTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTgc                                                                                                         <  1:1104915/133‑1 (MQ=255)
                                                             cTGAAATCAGTCAGAAATACAAGGGTTGCAGCCCTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTgc                                                                                                         >  2:1104915/1‑133 (MQ=255)
                                                             cTGAAATCAGTCAGAAATACAAGGGTTGCAGCCCTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTGc                                                                                        <  1:1340172/150‑1 (MQ=255)
                                                                      gtcagAAATACAAGGGTTGCAGCCCTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTTTACTCGGCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTGCGAGATAAAg                                                                               <  1:885214/150‑1 (MQ=255)
                                                                          gAAATACAAGGGTTGCAGCCCTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTGCGAGATAAAGCGAt                                                                           <  1:1370835/150‑1 (MQ=255)
                                                                            aaTACAAGGGTTGCAGCCCTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTGCGAGATAAAGCGATGt                                                                         >  2:754936/1‑150 (MQ=255)
                                                                             aTACAAGGGTTGCAGCCGTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTGCGAGATAAAGCGATGTc                                                                        <  1:1237791/150‑1 (MQ=255)
                                                                             aTACAAGGGTTGCAGCCCTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTGCGAGATAAAGCGATGTc                                                                        <  1:341004/150‑1 (MQ=255)
                                                                                cAAGGGTTGCAGCCCTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTGCGAGATAAAGCGATGTCGtt                                                                     >  2:477419/1‑150 (MQ=255)
                                                                                cAAGGGTTGCAGCCCTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTCTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTGCGAGATAAAGCGATGTCGtt                                                                     <  1:965995/150‑1 (MQ=255)
                                                                                cAAGGGTTGCAGCCCTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTCTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTGCGAGATAAAGCGATGTCGtt                                                                     <  1:966043/150‑1 (MQ=255)
                                                                                 aaGGGTTGCAGCCCTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTGCGAGATAAAGCGATGTCGTTg                                                                    >  1:1411051/1‑150 (MQ=255)
                                                                                 aaGGGTTGCAGCCCTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTGCGAGATAAAGCGATGTCGTTg                                                                    >  2:159748/1‑150 (MQ=255)
                                                                                  aGGGTTGCAGCCCTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTGCGAGATAAAGCGATGTCGTTGa                                                                   <  1:754936/150‑1 (MQ=255)
                                                                                   gggTTGCAGCCCTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTgtagta                                                                                                                                              <  1:1164060/74‑1 (MQ=255)
                                                                                   gggTTGCAGCCCTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTgtagta                                                                                                                                              >  2:1164060/1‑74 (MQ=255)
                                                                                   gggTTGCAGCCCTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTGCGAGATAAAGCGATGTCGTTGAc                                                                  >  2:689053/1‑150 (MQ=255)
                                                                                           gCCCTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTGCGAGATAAAGCGATGTCGTTGACTGCGagag                                                          <  1:477419/150‑1 (MQ=255)
                                                                                            cccTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTGCga                                                                                      <  1:854173/121‑1 (MQ=255)
                                                                                            cccTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTGCga                                                                                      >  2:854325/1‑121 (MQ=255)
                                                                                            cccTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTGCga                                                                                      >  2:854173/1‑121 (MQ=255)
                                                                                            cccTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTGCga                                                                                      <  1:854325/121‑1 (MQ=255)
                                                                                                 gaagctttcttGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTCCGAGATAAAGCGATGTCGTTGACTGCGAGAGGCCCGt                                                    >  1:410322/5‑150 (MQ=255)
                                                                                                      tttcttGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTGCgaga                                                                                    <  2:188926/113‑1 (MQ=255)
                                                                                                      tttcttGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTGCgaga                                                                                    >  1:188926/1‑113 (MQ=255)
                                                                                                      tttcttGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTGCGAGATAAAGCGATGTCGTTGACTGCGAGAGGCCCGTGCGct                                               >  2:540494/1‑150 (MQ=255)
                                                                                                             cAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTGCGAGATAAAGCGATGTCGTTGACTGCGAGAGGCCCGTGCGCTCTTGgat                                        >  2:676971/1‑150 (MQ=255)
                                                                                                             cAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTGCGAGATAAAGCGATGTCGTTGACTGCGAGAGGCCCGTGCGCTCTTGgat                                        >  1:695156/1‑150 (MQ=255)
                                                                                                                 gcCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTGCGAGATAAAGCGATGTCGTTGACTGCGAGAGGCCCGTGCGCTCTTGGATGATa                                    <  1:619056/150‑1 (MQ=255)
                                                                                                                   cTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTGCGAGATAAAGCGATGTCGTTGACTGCGAGAGGCCCGTGCGCTCTTGGATGATa                                    <  1:1278372/148‑1 (MQ=255)
                                                                                                                   cTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTGCGAGATAAAGCGATGTCGTTGACTGCGAGAGGCCCGTGCGCTCTTGGATGATa                                    >  2:1278372/1‑148 (MQ=255)
                                                                                                                    ttGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTGCGAGATAAAGCGATGTCGTTGACTGCgaga                                                           >  2:805182/1‑124 (MQ=255)
                                                                                                                    ttGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTGCGAGATAAAGCGATGTCGTTGACTGCGAGAGGCCCGTGCGCTCTTGGATGATACTg                                 >  2:482125/1‑150 (MQ=255)
                                                                                                                    ttGAAGAAGGAATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTGCGAGATAAAGCGATGTCGTTGACTGCgaga                                                           <  1:805182/124‑1 (MQ=255)
                                                                                                                      gaagaaGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTGCGAGATAAAGCGATGTCGTTGACTGCGAGAGGCCCGTGCGCTCTTGGATGATACTGAc                               <  2:1411051/150‑1 (MQ=255)
                                                                                                                       aagaagGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTGCGAGATAAAGCGATGTCGTTGACTGCGAGAGGCCCGTGCGCTCTTGGATGATACTGACg                              <  2:410322/150‑1 (MQ=255)
                                                                                                                               aTCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTGCGAGATAAAGCGATGTCGTTGACTGCga                                                             >  1:1400535/1‑111 (MQ=255)
                                                                                                                               aTCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTGCGAGATAAAGCGATGTCGTTGACTGCga                                                             <  2:1400535/111‑1 (MQ=255)
                                                                                                                                             cGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTGCGAGATAAAGCGATGTCGTTGACTGCGAGAGGCCCGTGCGCTCTTGGATGATACTGACGCACACACCTACCGTTTCTGGGt        <  2:695156/150‑1 (MQ=255)
                                                                                                                                             cGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTGCGAGATAAAGCGATGTCGTTGACTGCGAGAGGCCCGTGCGCTCTTGGATGATACTGACGCACACACCTACCGTTTCTGGGt        <  1:569425/150‑1 (MQ=255)
                                                                                                                                             cGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTGCGAGATAAAGCGATGTCGTTGACTGCGAGAGGCCCGTGCGCTCTTGGATGATACTGACGCACACACCTACCGTTTCTGGGt        <  1:676971/150‑1 (MQ=255)
                                                                                                                                                 gCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTGCGAGATAAAGCGATGTCGTTGACTGCGAGAGGCCCGTGCGCTCTTGGATGATACTGACGCACACACCTACCGTTTCTGGGTCGAc    <  1:469666/150‑1 (MQ=255)
                                                                                                                                                   gTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAAGTGCGCGCGCTGTAGAGCTGCGAGATAAAGCGATGTCGTTGACTGCGAGAGGCCCGTACGCTCTTGGATGATACTGACGCACACACCTACCGTTTCTGGGTCGAcct  >  1:1248115/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CGAAGCTTCGGATGATCTGATCGTTATTTGGTCAGATTGCTGTTTTTGACATCGGCAAAGACTGAAATCAGTCAGAAATACAAGGGTTGCAGCCCTAGATTCTTTCTTGCAGCGCCTTGAAGAAGGCATCGATATTCCCATCGTGGCGTTTGTAGTAGGTCCAGGGGCCGATCCGCTGCGATGTCACCAATTGCGCGCGCTGTAGAGCTGCGAGATAAAGCGATGTCGTTGACTGCGAGAGGCCCGTGCGCTCTTGGATGATACTGACGCACACACCTACCGTTTCTGGGTCGACCT  >  NC_002947/2831174‑2831470

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: