Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,845,651:1 +A 11.7% intergenic (+20/‑82) coaD → / → yfhL phosphopantetheine adenylyltransferase/4Fe‑4S cluster‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,845,6511.A11.7% 198.5 / 24.1 60intergenic (+20/‑82)coaD/yfhLphosphopantetheine adenylyltransferase/4Fe‑4S cluster‑containing protein
Reads supporting (aligned to +/- strand):  ref base . (24/29);  new base A (3/4);  total (27/33)
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.96e-01

GTCGGAGCGTTATTCGTTCATTTCCTCGACCCTGGTCCGGGAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGATCATCACCGACGATTGCATCAATTGCGACGTCTGCGAACCCGAGTGCCCGA  >  NC_002947/5845502‑5845793
                                                                                                                                                      |                                                                                                                                              
gTCGGAGCGTTATTCGTTCATTTCCTCGACCCTGGTCCGGGAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCcgc                                                                                                                                                 >  2:936299/1‑150 (MQ=255)
gTCGGAGCGTTATTCGTTCATTTCCTCGACCCTGGTCCGGGAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCcgc                                                                                                                                                 >  2:1024010/1‑150 (MQ=255)
  cGGAGCGTTATTCGTTCATTTCCTCGACCCTGGTCCGGGAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGACGGTATCCTGACGCCTGATTGTCGCGCCCGCGt                                                                                                                                              <  1:874425/150‑1 (MQ=255)
            ttCGTTCATTTCCTCGACCCTGGTCCGGGAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCCCCCAGTGTTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCggg                                                                                                                                    <  1:477598/150‑1 (MQ=255)
              cGTTCATTTCCTCGACCCTGGTCCGGGAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGgcg                                                                                                                                  >  2:811614/1‑150 (MQ=255)
               gTTCATTTCGTCGAGCCTGGTCCGGGAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGgcgc                                                                                                                                 >  2:983064/1‑150 (MQ=255)
               gTTCATTTCCTCGACCCTGGTCCGGGAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGgcgc                                                                                                                                 >  1:530394/1‑150 (MQ=255)
                ttCATTTCCTCGACCCTGGTCCGGGAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGGcgcg                                                                                                                                >  2:194293/1‑150 (MQ=255)
                ttCATTTCCTCGACCCTGGTCCGGGAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGGcgcg                                                                                                                                >  2:318427/1‑150 (MQ=255)
                 tCATTTCCTCGACCCTGGTCCGGGAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGTCTGACGCCTGATTGTCGCGCCCGCAGTGCACGCGGGCGCGa                                                                                                                               <  1:517536/150‑1 (MQ=255)
                    tttCCTCGACCCTGGTCCGGGAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCNAGCCTGNCGCCTGATTGTCNCGCCCGCGTGCACTGCGGGCGCGAATg                                                                                                                            <  1:983064/150‑1 (MQ=255)
                    tttCCTCGACCCTGGTCCGGGAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATg                                                                                                                            >  2:833934/1‑150 (MQ=255)
                      tCCTCGACCCTGGTCCGGGAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCg                                                                                                                          <  1:1024010/150‑1 (MQ=255)
                      tCCTCGACCCTGGTCCGGGAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCg                                                                                                                          <  1:464134/150‑1 (MQ=255)
                       ccTCGACCCTGGTCCGGGAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCgg                                                                                                                                     <  2:584484/138‑1 (MQ=255)
                       ccTCGACCCTGGTCCGGGAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCgg                                                                                                                                     >  1:584484/1‑138 (MQ=255)
                       ccTCGACCCTGGTCCGGGAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCgg                                                                                                                         <  2:785619/150‑1 (MQ=255)
                       ccTCGACCCTGGTCCGGGAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCgg                                                                                                                         >  1:86921/1‑150 (MQ=255)
                         tCGACCCTGGTCCGGGAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCGGca                                                                                                                       >  1:307524/1‑150 (MQ=255)
                                tGGTCCGGGAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCGGCACAATtgt                                                                                                                >  1:883322/1‑150 (MQ=255)
                                    ccGGGAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCGGCACAATTGTGccc                                                                                                            <  2:883322/150‑1 (MQ=255)
                                      gggAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCAt                                                                                                          >  1:367489/1‑150 (MQ=255)
                                       ggAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCAtt                                                                                                         <  2:668861/150‑1 (MQ=255)
                                         agaTTGCAGCGCTGGGCGGTGATATAAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGc                                                                                                       <  1:318427/150‑1 (MQ=255)
                                           aTTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGt                                                                                                     >  2:414019/1‑150 (MQ=255)
                                           aTTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGt                                                                                                     >  2:1119214/1‑150 (MQ=255)
                                            ttGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGtt                                                                                                    <  1:936299/150‑1 (MQ=255)
                                               cAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGaa                                                                                                 <  2:846472/150‑1 (MQ=255)
                                                aGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCAGTGCACGCGGGCGCGa                                                                                                                               <  2:130512/119‑1 (MQ=255)
                                                aGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCAGTGCACGCGGGCGCGa                                                                                                                               >  1:130512/1‑119 (MQ=255)
                                                 gcgcTGGGCGGTTATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGATGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACa                                                                                               <  1:311626/150‑1 (MQ=255)
                                                    cTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCAGTGCACGCGGgcgc                                                                                                                                 >  2:825238/1‑113 (MQ=255)
                                                    cTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCAGTGCACGCGGgcgc                                                                                                                                 <  1:825238/113‑1 (MQ=255)
                                                             gATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCtgt                                                                                   <  2:650165/150‑1 (MQ=255)
                                                                   aGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCGGCACAATtgt                                                                                                                <  1:807347/115‑1 (MQ=255)
                                                                   aGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCGGCACAATtgt                                                                                                                >  2:807347/1‑115 (MQ=255)
                                                                   aGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGccc                                                                             >  2:1130543/1‑150 (MQ=255)
                                                                    gCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGcccc                                                                            <  1:650215/150‑1 (MQ=255)
                                                                            gTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGATGCCTGATTGTCGCGCCCGCGTGCACTGCGGGAGCGAATGCGGAACAATTGTGCCCATTGCGTTGCACATGCCCGGGCTGTGAGCCccggcgcgcg                                                                    >  2:401050/1‑145 (MQ=255)
                                                                               cACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTcc                                                                 <  1:194293/150‑1 (MQ=255)
                                                                                       ggtggCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTcc                                                         <  2:1130452/150‑1 (MQ=255)
                                                                                             cGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGatc                                                   <  1:62373/150‑1 (MQ=255)
                                                                                             cGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGatc                                                   <  1:1130543/150‑1 (MQ=255)
                                                                                              gACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGatca                                                  <  1:297060/150‑1 (MQ=255)
                                                                                                     tgggggaACGCTGCGATCAGCACGCCAGACGCATGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGATCATCACcga                                           <  2:530394/145‑1 (MQ=255)
                                                                                                       actgaACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGATCATCACcga                                           <  1:505377/148‑1 (MQ=255)
                                                                                                       actgaACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGATCATCACcga                                           >  2:505377/1‑148 (MQ=255)
                                                                                                          gaACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGATCATCACCGACGAtt                                      <  2:367489/150‑1 (MQ=255)
                                                                                                          gaACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGATCATCACCGACGAtc                                      <  1:745461/150‑2 (MQ=255)
                                                                                                                     gaagCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGt                                                                   >  1:787418/1‑110 (MQ=255)
                                                                                                                     gaagCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGt                                                                   <  2:787418/110‑1 (MQ=255)
                                                                                                                     gaagCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAACCCCGCACGGAGGACCCATGGCCCTTAGCATGACCGACGATTGCATCAATTGc                           >  1:408627/1‑150 (MQ=255)
                                                                                                                                 aCGCCTGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGATCATCACCGACGATTGCATc                                 <  2:1129847/132‑1 (MQ=255)
                                                                                                                                 aCGCCTGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGATCATCACCGACGATTGCATc                                 >  1:1129847/1‑132 (MQ=255)
                                                                                                                                      tGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCtgt                                                                                   >  2:963688/1‑77 (MQ=255)
                                                                                                                                      tGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCtgt                                                                                   <  1:963688/77‑1 (MQ=255)
                                                                                                                                        aTTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGATCATCACCGACGATTGCATCAATTGCGACGTCTGCGAACCCGAGt        <  2:307524/150‑1 (MQ=255)
                                                                                                                                         ttGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGATCATCACCGACGATTGCATCAATTGCGACGTCTGCGAACCCGAGTg       >  2:962880/1‑150 (MQ=255)
                                                                                                                                          attCGCGCCCGCAGTGCACGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGATCATCACCGACGATTGCATCAATTGCGAc                        <  2:692559/130‑1 (MQ=255)
                                                                                                                                          attCGCGCCCGCAGTGCACGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGATCATCACCGACGATTGCATCAATTGCGAc                        >  1:692559/3‑132 (MQ=255)
                                                                                                                                              gcgcCCGCGTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGATCATCACCGACGATTGCATCAATTGCGACGTCTGCGAACCCGAGTGCCCGa  >  1:842282/1‑150 (MQ=255)
                                                                                                                                              gcgcCCGCGTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGATCATCACCGACGATTGCATCAATTGCGACGTCTGCGAACCCGAGTGCCCGa  >  2:1099365/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                              
GTCGGAGCGTTATTCGTTCATTTCCTCGACCCTGGTCCGGGAGATTGCAGCGCTGGGCGGTGATATCAGCAAGTTCGTCCACCCGGTGGTGGCCGACGCACTGACTGAACGCTTCAAGAAGCAAGCCTGACGCCTGATTGTCGCGCCCGCGTGCACTGCGGGCGCGAATGCGGCACAATTGTGCCCATTGCGTTGAACATGCCCGGGCTGTGAGCCCCGGCCGGAGTCCCCATGTCCCTGATCATCACCGACGATTGCATCAATTGCGACGTCTGCGAACCCGAGTGCCCGA  >  NC_002947/5845502‑5845793

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