Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,778,570 A→G 6.0% N542D (AAC→GAC)  PP_3340 → TonB‑dependent receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,778,5700AG 6.0% 192.5 / 2.4 66N542D (AAC→GAC) PP_3340TonB‑dependent receptor
Reads supporting (aligned to +/- strand):  ref base A (25/37);  new base G (2/2);  total (27/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

AGTTCACCCCGGCCAGCTGGGTCGATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGAAATCGATCACGTGCCACACCATCTGTGGAACACCGG  >  NC_002947/3778424‑3778709
                                                                                                                                                  |                                                                                                                                           
aGTTCACCCCGGCCAGCTGGGTCGATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAAc                                                                                                                                           <  4:223331/149‑1 (MQ=255)
 gTTCACCCCGGCCAGCTGGGTCGATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAAcc                                                                                                                                          <  6:178394/149‑1 (MQ=255)
  ttCACCCCGGCCAGCTGGGTCGATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCt                                                                                                                                         >  3:258329/1‑149 (MQ=255)
  ttCACCCCGGCCAGCTGGGTCGATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCt                                                                                                                                         <  3:234791/149‑1 (MQ=255)
    cACCCCGGCCAGCTGGGTCGATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTg                                                                                                                                        <  6:309226/148‑1 (MQ=255)
     aCCCCGGCCAGCTGGGTCGATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGc                                                                                                                                       <  8:71042/148‑1 (MQ=255)
      ccccGGCCAGCTGGGTCGATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTTGACCTGCag                                                                                                                                     <  2:233348/149‑2 (MQ=255)
       cccGGCCAGCTGGGTCGATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTTGACCTGCagg                                                                                                                                    >  4:3837/1‑147 (MQ=255)
       cccGGCCAGCTGGGTCGATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCAcc                                                                                                                                    >  4:122010/1‑149 (MQ=255)
       cccGGCCAGCTGGGTCGATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCAcc                                                                                                                                    >  5:12952/1‑149 (MQ=255)
        ccGGCCAGCTGGGTCGATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTTGACCTGCaggt                                                                                                                                   <  3:3837/149‑4 (MQ=255)
        ccGGCCAGCTGGGTCGATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCAccc                                                                                                                                   >  2:102332/1‑149 (MQ=255)
                cTGGGTCGATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCGGGTCAACCTGCACCCGGacgaac                                                                                                                           <  6:275006/149‑1 (MQ=255)
                cTGGGTCGATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGacgaac                                                                                                                           <  6:338365/149‑1 (MQ=255)
                    gTCGATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACgaacgaac                                                                                                                       >  8:190881/1‑149 (MQ=255)
                       gATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCGAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGa                                                                                                                    <  5:364106/149‑1 (MQ=255)
                            ccGCGAGACCTATTGCCGTCAGCACTCCTCTCCCGAAGTCAGCCGTAGCCAGAACGACACCACAGGCGACTCCGACAACTTCGCCAAAACACCCCGCTGGTGCTATCACCTGCAGGTCAACCTGCACACCGACGCACGAACCCAGgtgt                                                                                                               >  4:195978/1‑149 (MQ=255)
                                gTGGCCTATTGGCGTCATGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGAt                                                                                                           <  3:116094/149‑1 (MQ=255)
                                     cTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTcc                                                                                                       >  1:79908/1‑148 (MQ=255)
                                       aTTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTa                                                                                                     >  4:40346/1‑148 (MQ=255)
                                        ttGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCTGCGAGTCGGACAACTTCGGCGAAACACGCCGATGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACTAACCGAGGTGTGGATGTCCTa                                                                                                     <  4:353771/147‑1 (MQ=255)
                                          ggCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTc                                                                                                  <  4:258329/148‑1 (MQ=255)
                                           gCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCt                                                                                                >  2:1311/1‑149 (MQ=255)
                                           gCGTCAGGACGCCTCTGGCGAAGTGAGCCGGCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTcc                                                                                                 <  3:293795/148‑1 (MQ=255)
                                             gTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTgg                                                                                              <  2:79908/149‑1 (MQ=255)
                                             gTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCACTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTgg                                                                                              <  8:173982/149‑1 (MQ=255)
                                                aGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAg                                                                                           <  1:102332/149‑1 (MQ=255)
                                                aGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAg                                                                                           <  1:137004/149‑1 (MQ=255)
                                                  gACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGacgaacg                                                                                                                          <  4:352042/116‑1 (MQ=255)
                                                     gCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTc                                                                                      >  5:118769/1‑149 (MQ=255)
                                                            gCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGAt                                                                                <  4:196772/148‑1 (MQ=255)
                                                            gCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATc                                                                               >  8:53779/1‑149 (MQ=255)
                                                               aaGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCtt                                                                            >  2:218844/1‑149 (MQ=255)
                                                                 gTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGa                                                                          <  4:259353/149‑1 (MQ=255)
                                                                    aGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGc                                                                        <  3:351176/148‑1 (MQ=255)
                                                                     gCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGcc                                                                       <  3:40346/148‑1 (MQ=255)
                                                                         tCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCcagca                                                                  <  3:196018/149‑1 (MQ=255)
                                                                         tCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCcagca                                                                  <  5:40574/149‑1 (MQ=255)
                                                                            ccTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGa                                                                          >  4:336875/1‑138 (MQ=255)
                                                                            ccTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGa                                                                          <  3:336875/138‑1 (MQ=255)
                                                                             cTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAgcgc                                                              <  1:218844/149‑1 (MQ=255)
                                                                              tGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGcgcg                                                             >  5:158112/1‑149 (MQ=255)
                                                                                  cGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGc                                                         <  6:12952/149‑1 (MQ=255)
                                                                                  cGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGc                                                         <  1:1311/149‑1 (MQ=255)
                                                                                   gACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTcc                                                                                     <  8:247382/120‑1 (MQ=255)
                                                                                   gACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTcc                                                                                     >  7:247382/1‑120 (MQ=255)
                                                                                      cccAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGc                                                                        >  1:288520/1‑130 (MQ=255)
                                                                                      cccAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGc                                                                        <  2:288520/130‑1 (MQ=255)
                                                                                      cccAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAgg                                                     >  1:337899/1‑149 (MQ=255)
                                                                                      cccAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAgg                                                     <  3:122010/149‑1 (MQ=255)
                                                                                       ccAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAgg                                                     >  5:381531/1‑148 (MQ=255)
                                                                                              cgAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAg                                             <  7:190881/149‑1 (MQ=255)
                                                                                              cgAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAg                                             <  6:158112/149‑1 (MQ=255)
                                                                                                     gACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAgg                                       >  5:129735/1‑148 (MQ=255)
                                                                                                     gACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAgg                                       >  1:238299/1‑148 (MQ=255)
                                                                                                          cGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGAAATCg                                 <  6:381531/149‑1 (MQ=255)
                                                                                                              ggCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGAAATCGATc                              <  8:25236/148‑1 (MQ=255)
                                                                                                              ggCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGAAATCGATCa                             <  2:337899/149‑1 (MQ=255)
                                                                                                                        gccgcTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGAAATCGATCACGTGCCacac                   >  2:337830/1‑149 (MQ=255)
                                                                                                                             agcgacgaTGACCTGAAGGTCAACCTGCACCCGTACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGAAATCGATCACGTGCCACACCATCt              <  7:254148/142‑1 (MQ=255)
                                                                                                                              ggggCTATGACCTGCAGGTCAACCTGTACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGAAATCGATCACGTGCCAAACCATCtg             >  7:54179/1‑149 (MQ=255)
                                                                                                                                 gCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGAAATCGATCACGTGCCACACCATCtgtg           >  4:68212/1‑148 (MQ=255)
                                                                                                                                    tccgggtGCAGGTTGACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGAAATCGATCACGTGCCACACCATCTGTGGaa        >  4:288921/7‑148 (MQ=255)
                                                                                                                                     tGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGAAATCGATCACGTGCCACACCATCTGTGGaa        <  7:53779/147‑1 (MQ=255)
                                                                                                                                     tGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGAAATCGATCACGTGCCACACCATCTGTGGAaca      >  6:252424/1‑149 (MQ=255)
                                                                                                                                          tGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGAAATCGATCACGTGCCACACCATCTGTGGAACACCgg  >  2:165105/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                           
AGTTCACCCCGGCCAGCTGGGTCGATGGCCGCGTGGCCTATTGGCGTCAGGACGCCTCTGGCGAAGTGAGCCGTCGCCTGAACGACCCCAGCGGCGAGTCGGACAACGTCGGCGAAACACGCCGCTGGGGCTATGACCTGCAGGTCAACCTGCACCCGGACGAACGAACCGAGGTGTGGATGTCCTACTCCTGGCAGTACTCCAAGATCCTTGAGCCCAGCAGCGCGTTGCCAGGCAGCAAAGGCCAGGAAATCGATCACGTGCCACACCATCTGTGGAACACCGG  >  NC_002947/3778424‑3778709

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