Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,489,830 Δ1 bp 9.5% coding (264/597 nt) PP_4825 → MarC family membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,489,8300G. 9.5% 174.2 / 15.3 63coding (264/597 nt)PP_4825MarC family membrane protein
Reads supporting (aligned to +/- strand):  ref base G (26/27);  new base . (2/4);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.71e-01

AACTGGCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATGAAAAACTGCT  >  NC_002947/5489682‑5489979
                                                                                                                                                     |                                                                                                                                                      
aaCTGGCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATgg                                                                                                                                                       <  1:1035392/150‑1 (MQ=255)
   tGGCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTTACTTTAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGccc                                                                                                                                                    >  1:576543/1‑150 (MQ=255)
   tGGCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGccc                                                                                                                                                    >  1:963623/1‑150 (MQ=255)
     gCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAg                                                                                                                                                  <  1:1134417/150‑1 (MQ=255)
       cTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CCATGCCCAggg                                                                                                                                                <  1:1605339/150‑1 (MQ=255)
       cTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CCATGCCCAggg                                                                                                                                                <  1:1605700/150‑1 (MQ=255)
       cTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAggg                                                                                                                                                >  1:568814/1‑150 (MQ=255)
        tGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTATTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCCAT‑GCCCAGGGc                                                                                                                                               <  1:591867/150‑1 (MQ=255)
         ggAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTATTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCCAT‑GCCCAGGGcc                                                                                                                                              <  2:992018/150‑2 (MQ=255)
            aGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCCAT‑GCCCAGGGccgaa                                                                                                                                           <  1:323416/150‑5 (MQ=255)
            aGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CCATGCCCAGGGccgaa                                                                                                                                           <  1:323062/150‑5 (MQ=255)
                 gcgTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCCAGCcggcg                                                                                                                                      >  2:1041071/1‑150 (MQ=255)
                 gcgTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTACGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCcggcg                                                                                                                                      >  2:1833002/1‑150 (MQ=255)
                         aGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGc                                                                                                                             <  2:704875/150‑1 (MQ=255)
                           ccTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAg                                                                                                                           <  2:1831947/150‑1 (MQ=255)
                           ccTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAg                                                                                                                           <  2:1831411/150‑1 (MQ=255)
                            cTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGa                                                                                                                          >  1:1861942/1‑150 (MQ=255)
                               ggTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGTGAGCTTCGGAATTTTCGGCATCACCGCATATTTCTTGCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACaa                                                                                                                       <  1:329285/150‑1 (MQ=255)
                               ggTCGCCAGCGTGTTGCTGTACCTATTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACaa                                                                                                                       <  2:963623/150‑1 (MQ=255)
                                        cGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTcagca                                                                                                              <  1:666746/150‑1 (MQ=255)
                                              gCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGAtgt                                                                                                        <  1:662342/150‑1 (MQ=255)
                                                tGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGATGTGa                                                                                                      <  2:568814/150‑1 (MQ=255)
                                                tGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGATGTGa                                                                                                      >  1:788437/1‑150 (MQ=255)
                                                  tACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGATGTGAcc                                                                                                    >  2:507327/1‑150 (MQ=255)
                                                        ttCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCCAT‑GCCCAGGGccg                                                                                                                                             <  1:271714/104‑3 (MQ=255)
                                                        ttCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCCAT‑GCCCAGGGccg                                                                                                                                             >  2:271714/1‑102 (MQ=255)
                                                               gggcgATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGATGTGa                                                                                                      <  2:1321183/135‑1 (MQ=255)
                                                               gggcgATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGATGTGa                                                                                                      >  1:1321183/1‑135 (MQ=255)
                                                               gggcgATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGa                                                                                       <  2:607127/150‑1 (MQ=255)
                                                                 gcgATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGAcc                                                                                     >  1:935720/1‑150 (MQ=255)
                                                                   gATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCAt                                                                                   >  1:560012/1‑150 (MQ=255)
                                                                    aTCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCAtt                                                                                  >  2:184115/1‑150 (MQ=255)
                                                                     tCTTCGGAATTTTCGGCATAACCGCCGATGCGTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTc                                                                                 <  2:576543/150‑1 (MQ=255)
                                                                               tttCGGCATCACCGCCGATGCCTTCCGCATCGGTGCAGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGt                                                                       <  1:1310449/150‑1 (MQ=255)
                                                                                 tCGGCATCACCGCCGAGGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATTGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTg                                                                     <  1:493337/150‑1 (MQ=255)
                                                                                 tCGGCACCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTg                                                                     <  2:1847881/150‑1 (MQ=255)
                                                                                   ggCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGc                                                                   >  2:140826/1‑150 (MQ=255)
                                                                                    gCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGcc                                                                  >  1:762069/1‑150 (MQ=255)
                                                                                     cctCACCGCCGATGCCTTCCGATTTGGTGCTGGTCGCGTCCTGTTTTGTTCTGCCTTGGG‑CATGCCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTCCATCAGGATGTGACCCTCGTACCGCTGAGGATTGCGCTGACCGGTGGccc                                                                 <  2:500412/148‑1 (MQ=255)
                                                                                        cACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCTGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGa                                                                                                           <  2:731751/105‑1 (MQ=255)
                                                                                        cACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGa                                                                                                           >  1:731751/1‑105 (MQ=255)
                                                                                         aCCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGc                                                             >  2:436238/1‑150 (MQ=255)
                                                                                         aCCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGc                                                             >  2:221738/1‑150 (MQ=255)
                                                                                          ccgccgATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGAGATCACACG‑TCTCAACTTTCAGCAGGATGTTACCACCGCACCGCTGAGCATTCGGCTGACCGTTGGCCCCGGCa                                                            >  1:34537/1‑150 (MQ=255)
                                                                                           cgccgATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAgg                                                                                                                              >  1:132496/1‑83 (MQ=255)
                                                                                           cgccgATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAgg                                                                                                                              <  2:132496/83‑1 (MQ=255)
                                                                                                    cTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGc                                                  >  1:1105829/1‑150 (MQ=255)
                                                                                                     ttCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGcc                                                 >  2:1405642/1‑150 (MQ=255)
                                                                                                           aTCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACGATCGGTGCCctgctg                                           >  2:612737/1‑150 (MQ=255)
                                                                                                            tCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTgg                                          <  1:184115/150‑1 (MQ=255)
                                                                                                                 gCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATg                                     <  1:1405642/150‑1 (MQ=255)
                                                                                                                  cTGGCTGCGTGTTGTTGGTTTCGGCCCGGGG‑CATGGCCAAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATCCCGCTGACCGTTGGCCCCGGCAACATCGGTGCCCTGCTGGTGATgg                                    <  2:431615/150‑1 (MQ=255)
                                                                                                                   tGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATggg                                   <  2:305526/150‑1 (MQ=255)
                                                                                                                    ggCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATgggg                                  <  2:935720/150‑1 (MQ=255)
                                                                                                                      cAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGt                                <  2:1105829/150‑1 (MQ=255)
                                                                                                                        gCGTGCTGTTCCTTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGa                              <  1:1833002/150‑1 (MQ=255)
                                                                                                                        gCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCACGGCCAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGa                              <  1:1041071/150‑1 (MQ=255)
                                                                                                                          gTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGagc                            <  2:34537/150‑1 (MQ=255)
                                                                                                                            gCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGagcca                          <  2:1861942/150‑1 (MQ=255)
                                                                                                                                ttCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCg                      >  1:1330791/1‑150 (MQ=255)
                                                                                                                                   aTTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCAc                   >  2:1768635/1‑150 (MQ=255)
                                                                                                                                          cccTGGGCCAT‑GCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGAGTGAGCCAGCCGCACTGGGATg            >  2:459564/1‑150 (MQ=255)
                                                                                                                                           ccTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATGa           >  1:925816/1‑150 (MQ=255)
                                                                                                                                            cTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATGaa          >  2:1846922/1‑150 (MQ=255)
                                                                                                                                                     ggCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATGAAAAActgct  <  2:762069/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
AACTGGCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGG‑CATGGCCCAGGGCAAGCCGGCGGT‑ACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATGAAAAACTGCT  >  NC_002947/5489682‑5489979

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