Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,224,601 A→C 6.4% M26R (ATG→AGG)  gltK ← glutamate/aspartate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,224,6010AC 6.4% 213.5 / 4.6 73M26R (ATG→AGG) gltKglutamate/aspartate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (32/35);  new base C (3/2);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

AGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAG  >  NC_002947/1224455‑1224743
                                                                                                                                                  |                                                                                                                                              
aGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGCCTCCCATGACCATc                                                                                                                                              >  8:78250/1‑149 (MQ=255)
   aCCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAg                                                                                                                                            <  2:35598/148‑1 (MQ=255)
    ccAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGtt                                                                                                                                          >  5:64428/1‑149 (MQ=255)
      aGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGttt                                                                                                                                         >  3:74584/1‑148 (MQ=255)
        gTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAg                                                                                                                                      <  6:55927/149‑1 (MQ=255)
         tGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAGTTGAAGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGCCAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCAGATGACTCCCCTGACGATCAGTTTGAgg                                                                                                                                     >  8:249929/1‑149 (MQ=255)
                cagcagcagCGGAATCGAGCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCAt                                                                                                                                >  4:362962/1‑147 (MQ=255)
                cagcagcagCGGAATCGAGCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCAt                                                                                                                                <  3:362962/147‑1 (MQ=255)
                cagcagcagCGGAATCGAGCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGa                                                                                                                              >  8:31747/1‑149 (MQ=255)
                  gcagcagcGGAATCGAGCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGt                                                                                                                                   <  5:28507/142‑1 (MQ=255)
                  gcagcagcGGAATCGAGCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGt                                                                                                                                   >  6:28507/1‑142 (MQ=255)
                  gcagcagcGGAATCGAGCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGAcc                                                                                                                            <  2:34097/149‑1 (MQ=255)
                   cagcagcGGAATCGAGCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCa                                                                                                                           >  4:203697/1‑149 (MQ=255)
                    tgcatcGGAATGGAGCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACATGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACGATCAGTTTGAgg                                                                                                                                     <  8:90611/137‑1 (MQ=255)
                    agcagcGGAATCGAGCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAgg                                                                                                                                     >  7:90611/1‑138 (MQ=255)
                    agcagcGGAATCGAGCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGGTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGGCATGACCAt                                                                                                                          <  5:37285/149‑1 (MQ=255)
                                    ggAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCTGGGAGggca                                                                                                          <  3:94277/149‑1 (MQ=255)
                                      aaGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGggcagg                                                                                                        >  7:188940/1‑149 (MQ=255)
                                       agtagtTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCTGGGAGggcatg                                                                                                        <  2:183828/148‑3 (MQ=255)
                                       agtagtTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGggcagg                                                                                                        >  7:352065/1‑148 (MQ=255)
                                           gtTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAggg                                                                                                   >  6:270068/1‑149 (MQ=255)
                                           gtTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACAATCAGTTTGAGCGTCATGACCATGCCCTCCCAGAGGGCAGTCAggg                                                                                                   >  8:120537/1‑149 (MQ=255)
                                              gACGTAGGCGCCGGCCAGGTTAGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCgg                                                                                                >  6:71136/1‑149 (MQ=255)
                                               aCGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCggg                                                                                               <  1:342606/149‑1 (MQ=255)
                                                 gTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGga                                                                                             <  3:203697/149‑1 (MQ=255)
                                                    ggCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCagag                                                                                                              >  7:109185/1‑129 (MQ=255)
                                                    ggCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCagag                                                                                                              <  8:109185/129‑1 (MQ=255)
                                                    ggCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGgatga                                                                                          >  6:249928/1‑149 (MQ=255)
                                                    ggCGCCGGCAAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCGAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGgatga                                                                                          <  7:156865/149‑1 (MQ=255)
                                                               ggTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGggcaggca                                                                                                      >  5:105827/1‑126 (MQ=255)
                                                               ggTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGggcaggca                                                                                                      <  6:105827/126‑1 (MQ=255)
                                                               ggTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGaaa                                                                               <  6:64428/149‑1 (MQ=255)
                                                                 ttCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAAt                                                                              >  3:86303/1‑148 (MQ=255)
                                                                      cagcagTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAAtccattt                                                                        <  8:352065/149‑1 (MQ=255)
                                                                       agcagTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAAtccatttc                                                                       <  3:102874/149‑1 (MQ=255)
                                                                           gTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAAtccatttccatt                                                                   <  4:74584/149‑1 (MQ=255)
                                                                                cTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACtt                                                              <  7:227922/149‑1 (MQ=255)
                                                                                 ttGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTg                                                             >  1:229384/1‑149 (MQ=255)
                                                                                    acgacAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGccc                                                          >  2:142707/1‑149 (MQ=255)
                                                                                       acAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCAt                                                                                                                                >  3:80262/1‑76 (MQ=255)
                                                                                       acAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCAt                                                                                                                                <  4:80262/76‑1 (MQ=255)
                                                                                             gcATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGa                                                  >  3:235458/1‑148 (MQ=255)
                                                                                               aTCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAAt                                               <  4:86303/149‑1 (MQ=255)
                                                                                                  aGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATca                                             <  5:249928/148‑1 (MQ=255)
                                                                                                  aGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATca                                             <  4:162765/148‑1 (MQ=255)
                                                                                                   gTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTAGCTGAAATCCATTTCCATTTACTTGCCCCCAACGGAAATcagg                                           >  7:316509/1‑149 (MQ=255)
                                                                                                     gCCATGAAGGTGCCCAGCACGATACCGCCGATGACTCCCATCACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATcaggc                                          <  5:10863/148‑1 (MQ=255)
                                                                                                      ccAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAAtcatc                                                                          >  8:343899/1‑112 (MQ=255)
                                                                                                      ccAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAAtcatc                                                                          <  7:343899/115‑4 (MQ=255)
                                                                                                      ccAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATcaggcca                                        <  4:175937/149‑1 (MQ=255)
                                                                                                        aGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATcaggccagg                                      >  8:261213/1‑149 (MQ=255)
                                                                                                        aGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATcaggccagg                                      >  6:193545/1‑149 (MQ=255)
                                                                                                             ggTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGACCAGGCACgg                                 >  6:56655/1‑149 (MQ=255)
                                                                                                              gttcccAGCACGATACCGCCGATGACTCCCATGACCTTCAGTTTGAGGGTCATGACCATGCCCTCCCAGACGGCAGGCAGGGCGTGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACgg                                 <  7:249929/145‑1 (MQ=255)
                                                                                                                  ccAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACt                            <  8:188940/149‑1 (MQ=255)
                                                                                                                  ccAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACt                            <  7:31747/149‑1 (MQ=255)
                                                                                                                    aGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACttt                          <  1:285046/149‑1 (MQ=255)
                                                                                                                       aCGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTttctt                       <  7:78250/149‑1 (MQ=255)
                                                                                                                          aTACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCAc                                   >  8:246241/1‑134 (MQ=255)
                                                                                                                          aTACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCAc                                   <  7:246241/134‑1 (MQ=255)
                                                                                                                              cgccgATGACTCCCATGACCATCAGTTTCAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGAcc                <  6:317795/149‑1 (MQ=255)
                                                                                                                                  gATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATgc            <  7:261213/149‑1 (MQ=255)
                                                                                                                                       tggtcATGACCCTCAAACTGATGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACt                                                               <  5:262873/89‑1 (MQ=37)
                                                                                                                                       tggtcATGACCCTCAAACTGATGGTAATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACt                                                               >  6:262873/5‑93 (MQ=37)
                                                                                                                                       cTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTttctt                       >  7:302721/1‑133 (MQ=255)
                                                                                                                                       cTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTttctt                       <  8:302721/133‑1 (MQ=255)
                                                                                                                                          tcATGACCCTCAAACTGATGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGag     >  1:296617/2‑148 (MQ=255)
                                                                                                                                          ccATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGag     >  7:40312/1‑148 (MQ=255)
                                                                                                                                           cATGACCCTCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGagca   <  8:17137/149‑1 (MQ=255)
                                                                                                                                           cATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACttt                          <  5:27770/126‑1 (MQ=255)
                                                                                                                                           cATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACttt                          >  6:27770/1‑126 (MQ=255)
                                                                                                                                            aTGACCCTCAAACTGATGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCTATTTACTTGCCCCCAACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGagcag  >  6:96458/1‑149 (MQ=255)
                                                                                                                                            aTGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGagcag  >  5:239360/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
AGAACCAGGTGATCACCAGCAGCAGCGGAATCGAGCGGAAGTAGTTGACGTAGGCGCCGGCCAGGTTCGACAGCAGTTTGCTTGACGACAGGCGCATCAGTGCCAGGATGGTGCCCAGCACGATACCGCCGATGACTCCCATGACCATCAGTTTGAGGGTCATGACCATGCCCTCCCAGAGGGCAGGCAGGGCGGGGATGATTTCGCTGAAATCCATTTCCATTTACTTGCCCCCCACGGAAATCAGGCCAGGCACGGCGACTTTCTTCTCGACCATGCGCATGAGCAG  >  NC_002947/1224455‑1224743

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