Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,363,063 G→T 7.8% L235I (CTC→ATC)  glmM ← phosphoglucosamine mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,363,0630GT 7.8% 197.5 / 4.5 62L235I (CTC→ATC) glmMphosphoglucosamine mutase
Reads supporting (aligned to +/- strand):  ref base G (26/31);  new base T (3/2);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

GTACCGACCACCCCGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGT‑GTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCACGGAATACGCTGGGCGCAACCTTGTA  >  NC_002947/5362921‑5363210
                                                                                                                                                |                                                                                                                                                   
gTAACGACCACCCCGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGT‑GTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGc                                                                                                                                               <  8:34667/149‑1 (MQ=255)
 tACCGACCACCCCGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGT‑GTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCa                                                                                                                                              <  7:49814/149‑1 (MQ=255)
   ccGACCACCCCGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGT‑GTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATg                                                                                                                                            >  5:120686/1‑149 (MQ=255)
      accaccCCGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGT‑GTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTcc                                                                                                                                         >  6:102776/1‑149 (MQ=255)
       cccccccGCCGTGCCGCTTGCCCGTTTCCTGCAGGTCCCGTTCAATGATGAACAGCCTCTCGTCAACGTCGCCGATGGCGCCCGT‑GTGGTCGAACTTCAGCACGCGAGCGCCGTCACCGTCACAT‑GTCATGCCGAGATCTGCATGTCCg                                                                                                                                        <  8:233515/146‑1 (MQ=255)
       ccaccCCGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGT‑GTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGATCTCGGCAttgcct                                                                                                                                        <  3:364262/149‑6 (MQ=255)
        caccCCGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGT‑GTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGa                                                                                                                                       <  4:232824/149‑1 (MQ=255)
         accCCGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGGTGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGT‑GTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGAc                                                                                                                                      >  7:24172/1‑149 (MQ=255)
         accCCGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGT‑GTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGAc                                                                                                                                      <  4:181852/149‑1 (MQ=255)
           cccGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAATAGCAGCTCGTCACCGTCGACGATGGCGCCCGT‑GTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGATATCGGCATGTACgcaca                                                                                                                                    >  8:328923/1‑145 (MQ=255)
           cccGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGT‑GTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGAc                                                                                                                                      >  3:304848/1‑147 (MQ=255)
                   gCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGT‑GTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGATCTCGGCAttgc                                                                                                                                          <  7:23708/135‑4 (MQ=255)
                   gCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGT‑GTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGATCTCGGCAttgc                                                                                                                                          >  8:23708/1‑132 (MQ=255)
                    cAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGT‑GTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCg                                                                                                                            <  1:220043/148‑1 (MQ=255)
                          gCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGT‑GTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTg                                                                                                                       >  7:99927/1‑147 (MQ=255)
                          gCCATGGTCCTGCAGGTCGCGAGCAATGATGAATAGCAGCTCGTCACCGTCGACGATGGCGCCCGT‑GTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGc                                                                                                                      <  7:328923/148‑1 (MQ=255)
                           ccATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGT‑GTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAg                                                                                                                    >  4:344158/1‑149 (MQ=255)
                           ccATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGT‑GTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAg                                                                                                                    <  4:304848/149‑1 (MQ=255)
                            cATGTTCCTGCAGGTCGCGAGCAATGATGAACATAAGCTCGTCCCCGACGACGATCGCGCCCGT‑GTGCTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGc                                                                                                                   >  1:235197/1‑149 (MQ=255)
                            cATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGT‑GTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGc                                                                                                                   >  8:105921/1‑149 (MQ=255)
                             aTGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGT‑GTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGc                                                                                                                   >  2:344069/1‑148 (MQ=255)
                                      cAGGTCGCGAGCAATGGTGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGT‑GTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGAtgt                                                                                                         <  8:24172/149‑1 (MQ=255)
                                      cAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGT‑GTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGAtgt                                                                                                         <  2:66353/149‑1 (MQ=255)
                                      cAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGT‑GTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGAtgt                                                                                                         <  7:105921/149‑1 (MQ=255)
                                      cAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGT‑GTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGAtgt                                                                                                         >  7:75553/1‑149 (MQ=255)
                                      cAGGTCGCGAGAAATGATGAACATCAGCTCGTCACCGTCGACGATGGCGCCCGT‑GTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGAtgt                                                                                                         <  8:57845/149‑1 (MQ=255)
                                        ggTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGT‑GTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGAtgtgt                                                                                                       <  6:62840/149‑1 (MQ=255)
                                        ggTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGT‑GTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGAtgtgt                                                                                                       <  3:266971/149‑1 (MQ=255)
                                        ggTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGT‑GTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGAtgtgt                                                                                                       <  8:188861/149‑1 (MQ=255)
                                              gAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTc                                                                                                    >  1:82693/1‑147 (MQ=255)
                                                         aaCAGCAGCTCGTCACCGTCGACGATGGCGCCCGT‑GTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGAtt                                                                                                               <  2:52541/124‑1 (MQ=255)
                                                         aaCAGCAGCTCGTCACCGTCGACGATGGCGCCCGT‑GTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGAtt                                                                                                               >  1:52541/1‑124 (MQ=255)
                                                           cagcagCGCGTCACCGTCGACGATCGCGCCCGT‑GTGGGCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTGCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCga                                                                                    <  2:235197/149‑2 (MQ=255)
                                                              cagcTCGTCACCGTCGACGATGGCGCCCGT‑GTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGa                                                                                 >  7:315345/1‑149 (MQ=255)
                                                              cagcTCGTCACCGTCGACGATGGCGCCCGT‑GTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGa                                                                                 >  1:271183/1‑149 (MQ=255)
                                                                gcTCGTCACCGTCGACGATGGCGCCCGT‑GTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATg                                                                               <  3:344158/149‑1 (MQ=255)
                                                                 cTCGTCACCGTCGACGATGGCGCCCGT‑GTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGt                                                                              <  5:102776/149‑1 (MQ=255)
                                                                    gTCACCGTCGACGATGGCGCCCGT‑GTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTc                                                                            <  8:75553/148‑1 (MQ=255)
                                                                     tCACCGTCGACGATGGCGCCCGT‑GTGGTCGACCATCAGGACGCGATCGCCGTCACCGTCAAAG‑GCTATGGCGAGATCTTCATGTCCGACCAGTACAGCGGCCTGCAGGGTTTCGATGTGTGTCGAGCCGGAACCTTCGTTGTTGTTCAg                                                                          <  6:140624/149‑1 (MQ=255)
                                                                                      gcCCGT‑GTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCAt                                                         <  8:99927/149‑1 (MQ=255)
                                                                                       cccGT‑GTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATg                                                        <  6:120686/149‑1 (MQ=255)
                                                                                       cccGT‑GTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATg                                                        <  1:344069/149‑1 (MQ=255)
                                                                                               ggTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGAGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACg                                                 <  4:35600/149‑1 (MQ=255)
                                                                                                gTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGtt                                                                                   >  6:46331/1‑114 (MQ=255)
                                                                                                gTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGtt                                                                                   <  5:46331/114‑1 (MQ=255)
                                                                                                gTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACgg                                                >  2:259221/1‑149 (MQ=255)
                                                                                                gTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACgg                                                >  5:252876/1‑149 (MQ=255)
                                                                                                       aTCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTc                                         <  2:52283/149‑1 (MQ=255)
                                                                                                        tCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCg                                        >  7:52799/1‑149 (MQ=255)
                                                                                                           gCACGCGATCGCCGTCACCGTCAAAGCGCAATGCCGAGATCGGCATGTCCGAc                                                                                                                                      <  1:191055/53‑1 (MQ=255)
                                                                                                           gCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGAc                                                                                                                                      >  2:191055/1‑52 (MQ=255)
                                                                                                           gCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCa                                     >  5:43114/1‑149 (MQ=255)
                                                                                                                  aTCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCt                              <  2:271183/149‑1 (MQ=255)
                                                                                                                   tCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTc                             >  3:266491/1‑149 (MQ=255)
                                                                                                                       cgtcaccgtcaAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCACg                          >  2:87985/1‑148 (MQ=255)
                                                                                                                        gtcaccgtcaAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGcgtattcc                        >  6:241604/1‑142 (MQ=255)
                                                                                                                        gtcaccgtcaAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCACGGa                        <  8:296855/149‑1 (MQ=255)
                                                                                                                            ccgtcaAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCACGGAATAc                    <  6:252876/149‑1 (MQ=255)
                                                                                                                               tcaAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCACGGAATACGCt                 >  2:256546/1‑149 (MQ=255)
                                                                                                                                   aG‑GCAATGCCGAGATCGGCATGTCAGACCAGTACAGCGGCGTGCAGCGATTCGATGTGTGTCGCGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTGAc                           >  7:307008/1‑135 (MQ=255)
                                                                                                                                       acaTGCCGATCTCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCACGGAATACGCTGGgcgc           >  1:184069/3‑148 (MQ=255)
                                                                                                                                             cGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCACGGAATACGCTGGGCGCAACCTTg    >  4:150671/1‑149 (MQ=255)
                                                                                                                                               agaTCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCACGGAATACGCTGGGCGCAACCTTgta  <  2:82693/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GTACCGACCACCCCGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGT‑GTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAG‑GCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCACGGAATACGCTGGGCGCAACCTTGTA  >  NC_002947/5362921‑5363210

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: ‑