Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,079,174 T→A 14.5% intergenic (‑371/+26) mgtE ← / ← PP_t56 magnesium transporter/tRNA‑Arg

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,079,1740TA14.5% 164.5 / 21.4 62intergenic (‑371/+26)mgtE/PP_t56magnesium transporter/tRNA‑Arg
Reads supporting (aligned to +/- strand):  ref base T (28/25);  new base A (5/4);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

GCTATCAGTGCGACGATTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGAAAT  >  NC_002947/5079027‑5079319
                                                                                                                                                   |                                                                                                                                                 
gCTATCAGTGCGACGATTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTAt                                                                                                                                                 >  1:1584534/1‑150 (MQ=255)
   aTCAGTGCGACGATTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGa                                                                                                                                              >  1:990643/1‑150 (MQ=255)
     cAGTGCGACGATTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGc                                                                                                                                            <  2:990643/150‑1 (MQ=255)
      aGTGCGACGATTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCg                                                                                                                                           <  1:1001411/150‑1 (MQ=255)
         gCGACGATTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGc                                                                                                                                        <  1:1126981/150‑1 (MQ=255)
          cgacgaTTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCt                                                                                                                                       >  2:37815/1‑150 (MQ=255)
            acgaTTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCttt                                                                                                                                     >  2:797870/1‑150 (MQ=255)
                ttgttgTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTgt                                                                                                                                 <  1:797870/150‑2 (MQ=255)
                      tgGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCAATAGAGCGGGCttttt                                                                                                                                   <  2:827258/142‑1 (MQ=255)
                      tgGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCAATAGAGCGGGCttttt                                                                                                                                   >  1:827258/1‑142 (MQ=255)
                             gAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCACACGAATTCATTGTGGCACAGACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCgg                                                                                                                    >  1:367635/1‑150 (MQ=255)
                              aaaGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGa                                                                                                                   >  1:286450/1‑150 (MQ=255)
                              aaaGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGa                                                                                                                   <  2:418718/150‑1 (MQ=255)
                                               tGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAAc                                                                                                  <  1:1632730/150‑1 (MQ=255)
                                                   gCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTcag                                                                                                              <  1:238903/134‑1 (MQ=255)
                                                   gCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTcag                                                                                                              >  2:238903/1‑134 (MQ=255)
                                                   gCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGc                                                                                              <  1:1364499/150‑1 (MQ=255)
                                                   gCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGc                                                                                              <  1:37815/150‑1 (MQ=255)
                                                     cGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCgg                                                                                                                                          >  1:1490230/1‑104 (MQ=255)
                                                     cGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCgg                                                                                                                                          <  2:1490230/104‑1 (MQ=255)
                                                        tCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCg                                                                                         >  1:1009653/1‑150 (MQ=255)
                                                         cAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCAATAGAGCGGGCttttt                                                                                                                                   >  2:750285/1‑107 (MQ=255)
                                                         cAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCAATAGAGCGGGCttttt                                                                                                                                   <  1:750285/107‑1 (MQ=255)
                                                          aGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCgg                                                                                                                                          >  1:1490148/1‑99 (MQ=255)
                                                          aGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCgg                                                                                                                                          <  2:1490148/99‑1 (MQ=255)
                                                           gCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTc                                                                                      >  2:34670/1‑150 (MQ=255)
                                                                  cAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGt                                                                               >  1:583035/1‑150 (MQ=255)
                                                                  cAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGt                                                                               >  2:1313936/1‑150 (MQ=255)
                                                                    aaCTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCAATAGAGCGGGCTTTTTgttt                                                                                                                               >  2:805666/1‑97 (MQ=38)
                                                                    aaCTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCAATAGAGCGGGCTTTTTgttt                                                                                                                               <  1:805666/100‑4 (MQ=38)
                                                                        acacCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGa                                                                         <  1:1313936/150‑1 (MQ=255)
                                                                           ccGAAATCATTACAGAAGAGCGTAGATCGACAGGCTAAAGAATTCATTGGCGAACAAACAAAAAGCGCGCTCTATGGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGCACCTGCGACCGCTCGGTTCGTAGCCGAGTa                                                                      <  1:300138/150‑1 (MQ=38)
                                                                              cAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTGTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTActc                                                                   >  2:132501/1‑150 (MQ=255)
                                                                                  cATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGAc                                                                                           <  2:1341923/122‑1 (MQ=255)
                                                                                  cATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGAc                                                                                           >  1:1341923/1‑122 (MQ=255)
                                                                                   aTTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATgg                                                                                                                       >  1:57522/1‑93 (MQ=255)
                                                                                   aTTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATgg                                                                                                                       <  2:57522/93‑1 (MQ=255)
                                                                                       cAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCt                                                          <  2:286450/150‑1 (MQ=255)
                                                                                            gagCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCt                                                     >  2:447105/1‑150 (MQ=255)
                                                                                                       gACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTg                                          >  2:1195776/1‑150 (MQ=255)
                                                                                                        aCAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGa                                                        <  1:426980/135‑1 (MQ=255)
                                                                                                        aCAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGa                                                        >  2:426980/1‑135 (MQ=255)
                                                                                                        aCAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGa                                         <  1:447105/150‑1 (MQ=255)
                                                                                                         cAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGAt                                        >  2:448465/1‑150 (MQ=255)
                                                                                                             cTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTa                                    >  1:1082444/1‑150 (MQ=255)
                                                                                                                aaGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATggg                                                                                                                             <  2:955105/58‑1 (MQ=255)
                                                                                                                aaGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATggg                                                                                                                             >  1:955105/1‑58 (MQ=255)
                                                                                                                aaGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGtt                                 >  2:555788/1‑150 (MQ=255)
                                                                                                                     ttCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGa                            <  2:583035/150‑1 (MQ=255)
                                                                                                                     ttCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGa                            <  2:594283/150‑1 (MQ=255)
                                                                                                                         ttGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGa                                                 <  2:416090/125‑1 (MQ=255)
                                                                                                                         ttGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGa                                                 >  1:416090/1‑125 (MQ=255)
                                                                                                                                   aacaaaAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTg              <  1:448465/150‑1 (MQ=255)
                                                                                                                                   aacaaaAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCACGTTACCACTg              >  2:5457/1‑150 (MQ=255)
                                                                                                                                      aaaaaGCCCGCTCAATAGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAg                                                <  1:1036262/113‑1 (MQ=255)
                                                                                                                                      aaaaaGCCCGCTCAATAGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAg                                                >  2:1036262/1‑113 (MQ=255)
                                                                                                                                           gCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTa         <  1:519916/147‑1 (MQ=255)
                                                                                                                                           gCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTa         >  2:519916/1‑147 (MQ=255)
                                                                                                                                            cccGCTCAATAGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGa     >  2:962793/1‑150 (MQ=255)
                                                                                                                                              cGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGaaa   >  2:474856/1‑150 (MQ=255)
                                                                                                                                               gCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGAAAt  <  1:555788/150‑1 (MQ=255)
                                                                                                                                               gCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTAGCCAACTAAGCTATGCGTCCGCGATGGGAATTTTATCACAGGGTAACAGTGTCTTAATGAAAt  >  2:645983/1‑150 (MQ=18)
                                                                                                                                                ctctATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTAt                                                   <  2:348566/100‑1 (MQ=255)
                                                                                                                                                ctctATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTAt                                                   >  1:348566/1‑100 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GCTATCAGTGCGACGATTGTTGTGGATGGGAAAGCACATTTGATTTTTGAGGCCGCTCAGCGGCTCCAAACTACACCGCAATCATTACAGAAGAGCGTAGATCGACAGGCTCAAGAATTCATTGGCGCACAAACAAAAAGCCCGCTCTATTGAGCGGGCTTTTTGATGGGGTATGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGATGGTAGTTTTAGACCAGGTTACCACTGGCTTACTGAAAT  >  NC_002947/5079027‑5079319

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