Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,851,650 C→T 94.1% R232W (CGG→TGG)  relA → ATP:GTP 3'‑pyrophosphotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,851,6500CT94.1% 199.7 / 2.1 68R232W (CGG→TGG) relAATP:GTP 3'‑pyrophosphotransferase
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (32/32);  total (34/34)
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.51e-01

TGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGA  >  NC_002947/1851508‑1851795
                                                                                                                                              |                                                                                                                                                 
tgtCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGTGGCGGCt                                                                                                                                            <  1:824715/150‑1 (MQ=255)
     cgcgTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTTCAAGCAGATCGCCAAGCTGTTGCATGAGTGGCGGCTGGAcc                                                                                                                                       >  2:1321802/1‑150 (MQ=255)
            aGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGTggcggc                                                                                                                                             >  2:874595/1‑137 (MQ=255)
            aGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGTggcggc                                                                                                                                             <  1:874595/137‑1 (MQ=255)
                ttCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTtc                                                                                                                            <  1:926738/150‑1 (MQ=255)
                ttCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCAAGAGTCGCGGCTAGACGGCGaacccctc                                                                                                                            >  2:1079518/1‑143 (MQ=255)
                             cgccgcTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGtga                                                                                                               <  1:412954/150‑1 (MQ=255)
                                cgcTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGtgatga                                                                                                            >  2:1282887/1‑150 (MQ=255)
                                cgcTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGtgatga                                                                                                            <  1:1871278/150‑1 (MQ=255)
                                  cTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAAc                                                                                                          <  1:568758/150‑1 (MQ=255)
                                   tCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAAcc                                                                                                         >  2:525225/1‑150 (MQ=255)
                                     gcgcACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAg                                                                                                       >  2:276034/1‑150 (MQ=255)
                                     gcgcACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAg                                                                                                       <  1:293382/150‑1 (MQ=255)
                                            gCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGc                                                                                                      <  2:411173/144‑1 (MQ=255)
                                            gCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGc                                                                                                      >  1:411173/1‑144 (MQ=255)
                                                 ggTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAg                                                                                           <  1:525225/150‑1 (MQ=255)
                                                  gTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGt                                                                                          <  2:195521/150‑1 (MQ=255)
                                                   tATCGATCACATCAAGTGGATGCTGGAAGGCCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGtt                                                                                         <  1:592357/150‑1 (MQ=255)
                                                       gATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGc                                                                                       <  2:588124/148‑1 (MQ=255)
                                                       gATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGc                                                                                       >  1:588124/1‑148 (MQ=255)
                                                       gATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTc                                                                                     >  1:372933/1‑150 (MQ=255)
                                                       gATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTc                                                                                     >  2:437239/1‑150 (MQ=255)
                                                           acaTCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCa                                                                                 <  1:1079518/150‑1 (MQ=255)
                                                                  gTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGt                                                                          >  1:70030/1‑150 (MQ=255)
                                                                         cTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGcc                                                                   >  1:1428583/1‑150 (MQ=255)
                                                                          tGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCtt                                                                                                                             <  2:1620081/91‑1 (MQ=255)
                                                                          tGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCtt                                                                                                                             >  1:1620081/1‑91 (MQ=255)
                                                                                  cTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGGGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGc                                                                    <  2:1486112/140‑1 (MQ=255)
                                                                                  cTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGGGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGc                                                                    >  1:1486112/1‑140 (MQ=255)
                                                                                    gTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCa                                                                                 >  2:1360239/1‑125 (MQ=255)
                                                                                    gTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCa                                                                                 <  1:1360239/125‑1 (MQ=255)
                                                                                       cTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTtca                                                                                                                           >  1:871265/1‑80 (MQ=255)
                                                                                       cTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTtca                                                                                                                           <  2:871265/80‑1 (MQ=255)
                                                                                       cTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCg                                                     <  1:1918258/150‑1 (MQ=255)
                                                                                        ttCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGa                                                                                                                                   <  1:1264896/71‑1 (MQ=255)
                                                                                        ttCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGa                                                                                                                                   >  2:1264896/1‑71 (MQ=255)
                                                                                        ttCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCgg                                                    >  2:541438/1‑150 (MQ=255)
                                                                                        ttCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCgg                                                    >  2:1153793/1‑150 (MQ=255)
                                                                                                   ggAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGGGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACatat                                         >  1:132919/1‑150 (MQ=255)
                                                                                                     aGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCt                                       <  2:372933/150‑1 (MQ=255)
                                                                                                       cccGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTAt                                     >  2:1770772/1‑150 (MQ=255)
                                                                                                         cGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTTGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTc                                   <  2:1486892/150‑1 (MQ=255)
                                                                                                         cGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTc                                   >  2:480460/1‑150 (MQ=255)
                                                                                                         cGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTc                                   >  1:1063029/1‑150 (MQ=255)
                                                                                                           aCCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTCGa                                 <  2:713910/150‑1 (MQ=255)
                                                                                                           aCCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACACAAGCGGCCGGGCGAAACATATCTATTCGa                                 <  1:1321802/150‑1 (MQ=255)
                                                                                                                   aaGCAGATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTACAGAACGAGTTGCTCGCCACGGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTCGATCTGgcgc                         >  2:462172/1‑150 (MQ=255)
                                                                                                                        gATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTCGATc                               >  2:1840788/1‑139 (MQ=255)
                                                                                                                        gATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTCGATc                               <  1:1840788/139‑1 (MQ=255)
                                                                                                                        gATCGCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGAt                    >  1:245235/1‑150 (MQ=255)
                                                                                                                            gCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGaa                                                                       >  2:10259/1‑95 (MQ=255)
                                                                                                                            gCCAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGaa                                                                       <  1:10259/95‑1 (MQ=255)
                                                                                                                             ccAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACatat                                         <  2:764326/124‑1 (MQ=255)
                                                                                                                             ccAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACatat                                         >  1:764326/1‑124 (MQ=255)
                                                                                                                             ccAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTAt                                     >  1:1662026/1‑128 (MQ=255)
                                                                                                                             ccAAGCTGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTAt                                     <  2:1662026/128‑1 (MQ=255)
                                                                                                                                   tGTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGATGCAGCGCAAAg         >  1:56765/1‑150 (MQ=255)
                                                                                                                                    gTTGCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACAGATCTATTCGATCTGGCGCAAGATGCAGCGCAAAgg        >  2:1039248/1‑150 (MQ=255)
                                                                                                                                      tGCATGGGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGcc      <  2:132919/150‑1 (MQ=255)
                                                                                                                                       gCATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCt     <  1:276034/150‑1 (MQ=255)
                                                                                                                                        cATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGAc                                                                <  1:171000/90‑1 (MQ=255)
                                                                                                                                        cATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGAc                                                                >  2:171000/1‑90 (MQ=255)
                                                                                                                                        cATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACatat                                         >  1:1249607/1‑113 (MQ=255)
                                                                                                                                        cATGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACatat                                         <  2:1249607/113‑1 (MQ=255)
                                                                                                                                        cATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAg                                                                      <  1:639955/84‑1 (MQ=255)
                                                                                                                                        cATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAg                                                                      >  2:639955/1‑84 (MQ=255)
                                                                                                                                          tGAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGa  <  1:480460/150‑1 (MQ=255)
                                                                                                                                          ggAGTGGCGGCTGGACCGCGAGCGCTTCATCAGCGAGGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGa  <  2:255643/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
TGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGA  >  NC_002947/1851508‑1851795

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: