Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,616,606 G→C 8.2% A171A (GCG→GCC)  bepA → OM protein maintenance and assembly metalloprotease and chaperone, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,616,6060GC 8.2% 340.7 / 17.7 98A171A (GCG→GCC) bepAOM protein maintenance and assembly metalloprotease and chaperone, periplasmic
Reads supporting (aligned to +/- strand):  ref base G (44/45);  new base C (4/4);  total (48/50)
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.83e-01

Only 100 of 101 total aligned reads displayed.

GTTATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCC  >  NC_000913/2616458‑2616749
                                                                                                                                                    |                                                                                                                                               
gTTATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCg                                                                                                                                                 >  2:154978/1‑149 (MQ=255)
 ttATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGa                                                                                                                                                <  3:327244/149‑1 (MQ=255)
 ttATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGa                                                                                                                                                >  3:40065/1‑149 (MQ=255)
   aTTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATg                                                                                                                                              <  4:33209/149‑1 (MQ=255)
              gAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGc                                                                                                                                            <  3:342782/140‑1 (MQ=255)
              gAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGc                                                                                                                                            >  4:342782/1‑140 (MQ=255)
              gAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCCATCGCCAGtaaaa                                                                                                                                    <  2:1131428/148‑5 (MQ=255)
              gAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCCATCGCCAGtaaaa                                                                                                                                    >  1:1131428/1‑144 (MQ=255)
                 aGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGc                                                                                                                                <  3:618205/149‑1 (MQ=255)
                  gTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCg                                                                                                                               <  2:183467/149‑1 (MQ=255)
                     aaCTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCgggg                                                                                                                            >  3:156886/1‑149 (MQ=255)
                      aCTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCCATCGCCAGtaaa                                                                                                                                     <  1:680040/139‑4 (MQ=255)
                      aCTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCCATCGCCAGtaaa                                                                                                                                     >  2:680040/1‑136 (MQ=255)
                        cggCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCgg                                                                                                                              <  2:798801/143‑1 (MQ=255)
                        cggCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCgg                                                                                                                              >  1:798801/2‑144 (MQ=255)
                                 ttATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGTGGCGCTg                                                                                                                >  3:937808/1‑149 (MQ=255)
                                     ggCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCg                                                                                                            <  4:156886/149‑1 (MQ=255)
                                          aCGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTaca                                                                                                       >  3:238119/1‑149 (MQ=255)
                                          aCGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTaca                                                                                                       >  1:1417466/1‑149 (MQ=255)
                                              aaTCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTgg                                                                                                   >  1:632434/1‑149 (MQ=255)
                                                  tcCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGc                                                                                                                                            >  1:844695/1‑104 (MQ=255)
                                                  tcCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGc                                                                                                                                            <  2:844695/104‑1 (MQ=255)
                                                  tcCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCggg                                                                                               >  1:372912/1‑149 (MQ=255)
                                                  tcCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCCATCGCCAGtaaaa                                                                                                                                    <  3:519203/112‑5 (MQ=255)
                                                  tcCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCCATCGCCAGtaaaa                                                                                                                                    >  4:519203/1‑108 (MQ=255)
                                                       cGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAAcgc                                                                                          <  4:40065/149‑1 (MQ=255)
                                                       cGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAAcgc                                                                                          <  2:632434/149‑1 (MQ=255)
                                                            cccAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAg                                                                                     >  1:261734/1‑149 (MQ=255)
                                                             ccAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAgg                                                                                    >  3:1398642/1‑149 (MQ=255)
                                                             ccAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAgg                                                                                    >  1:1459524/1‑149 (MQ=255)
                                                               aaCGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAgggg                                                                                  >  3:475856/1‑149 (MQ=255)
                                                                aCGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGga                                                                                 >  1:1439539/1‑149 (MQ=255)
                                                                 cGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGgat                                                                                <  2:1417466/149‑1 (MQ=255)
                                                                    cACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGgatgat                                                                             <  1:154978/149‑1 (MQ=255)
                                                                    cACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGgatgat                                                                             >  1:469773/1‑149 (MQ=255)
                                                                      ccTGGCGCGAGCGATGGAAGATCAGCAGCGCAGGGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCa                                                                           <  2:249880/149‑1 (MQ=255)
                                                                       gtGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGt                                                                                        >  1:962273/2‑135 (MQ=255)
                                                                       cTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGt                                                                                        <  2:962273/135‑1 (MQ=255)
                                                                         ggCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGtt                                                                        <  4:1398642/149‑1 (MQ=255)
                                                                         ggCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGtt                                                                        <  4:475856/149‑1 (MQ=255)
                                                                                gcgaTGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGAt                                                                                                                          >  2:1259168/1‑92 (MQ=255)
                                                                                gcgaTGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGCGAt                                                                                                                          <  1:1259168/92‑1 (MQ=255)
                                                                                    tGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAgg                                                                                    <  2:1002984/126‑1 (MQ=255)
                                                                                    tGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAgg                                                                                    >  1:1002984/1‑126 (MQ=255)
                                                                                      gAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCgg                                                                                                >  2:1005493/1‑112 (MQ=255)
                                                                                      gAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCgg                                                                                                <  1:1005493/112‑1 (MQ=255)
                                                                                         ggtCAGCGGCGCGCTGCGCTGCTGATCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGa                                                        >  1:483521/3‑149 (MQ=255)
                                                                                          aTCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGaa                                                       <  2:1459524/149‑1 (MQ=255)
                                                                                           tCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCa                                                                           >  1:1049652/1‑128 (MQ=255)
                                                                                           tCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCa                                                                           <  2:1049652/128‑1 (MQ=255)
                                                                                            cagcagCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACa                                                     <  2:469773/149‑1 (MQ=255)
                                                                                              gcagcgcagcGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAgg                                                   <  2:261734/149‑1 (MQ=255)
                                                                                              gcagcgcagcGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAgg                                                   >  3:774373/1‑149 (MQ=255)
                                                                                                 gcgcagcgcGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAg                                                <  2:1439539/149‑1 (MQ=255)
                                                                                                  cgcagcgcGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGc                                               <  4:238119/149‑1 (MQ=255)
                                                                                                            cgcTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCAtt                                     >  1:18952/1‑149 (MQ=255)
                                                                                                              cTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTgg                                   <  2:372912/149‑1 (MQ=255)
                                                                                                                   cTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAg                                                >  3:454692/1‑131 (MQ=255)
                                                                                                                   cTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAg                                                <  4:454692/131‑1 (MQ=255)
                                                                                                                   cTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTc                              >  2:813096/1‑149 (MQ=255)
                                                                                                                   cTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTc                              >  1:1205855/1‑149 (MQ=255)
                                                                                                                   cTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTc                              <  1:813096/149‑1 (MQ=255)
                                                                                                                   cTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTc                              >  1:1124185/1‑149 (MQ=255)
                                                                                                                   cTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGACCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTc                              >  3:662572/1‑149 (MQ=255)
                                                                                                                      ggTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCgg                                                                                                                              >  2:623667/1‑50 (MQ=255)
                                                                                                                      ggTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCgg                                                                                                                              <  1:623667/50‑1 (MQ=255)
                                                                                                                      ggTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAgg                           >  1:338220/1‑149 (MQ=255)
                                                                                                                       gTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGCAAGCGCACCGCATTGGTATTCAGGt                          <  2:483521/149‑1 (MQ=255)
                                                                                                                          ggCGCGTTAGGTTCTATTTTACTGGCTATGGCCAGTCCGCAGGCGGGGATCGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAg                            <  1:416643/144‑1 (MQ=255)
                                                                                                                          ggCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAg                            >  2:416643/1‑144 (MQ=255)
                                                                                                                              cgTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCaa                   <  2:1205855/149‑1 (MQ=255)
                                                                                                                              cgTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCaa                   <  4:662572/149‑1 (MQ=255)
                                                                                                                               gTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGgatgat                                                                             <  1:1236449/90‑1 (MQ=255)
                                                                                                                               gTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGgatgat                                                                             >  2:1236449/1‑90 (MQ=255)
                                                                                                                                    gTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCg             >  3:1469078/1‑149 (MQ=255)
                                                                                                                                    gTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCg             >  2:1304449/1‑149 (MQ=255)
                                                                                                                                     ttCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGAcc                                          <  4:812389/119‑1 (MQ=255)
                                                                                                                                     ttCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGAcc                                          >  3:812389/1‑119 (MQ=255)
                                                                                                                                     ttCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCgg            >  3:919417/1‑149 (MQ=255)
                                                                                                                                     ttCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCgg            <  1:199774/149‑1 (MQ=255)
                                                                                                                                     ttCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCgg            >  1:856637/1‑149 (MQ=255)
                                                                                                                                     ttCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCgg            >  2:199774/1‑149 (MQ=255)
                                                                                                                                          ttttACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTc       >  3:786469/1‑149 (MQ=255)
                                                                                                                                          ttttACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTc       <  2:856637/149‑1 (MQ=255)
                                                                                                                                          ttttACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTc       <  4:774373/149‑1 (MQ=255)
                                                                                                                                            ttACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGa     <  1:1304449/149‑1 (MQ=255)
                                                                                                                                             tACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACg                 <  1:1020347/136‑1 (MQ=255)
                                                                                                                                             tACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACg                 >  2:1020347/1‑136 (MQ=255)
                                                                                                                                             tACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGAt    <  2:1124185/149‑1 (MQ=255)
                                                                                                                                             gaCTGGCCATCGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAAt                                                          >  4:336733/2‑95 (MQ=255)
                                                                                                                                             gaCTGGCCATCGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAAt                                                          <  3:336733/94‑1 (MQ=255)
                                                                                                                                              aCTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATc   <  2:338220/149‑1 (MQ=255)
                                                                                                                                                tGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGAtt        <  3:174928/142‑1 (MQ=255)
                                                                                                                                                tGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGAtt        >  4:174928/1‑142 (MQ=255)
                                                                                                                                                tGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATcc  <  4:159420/148‑1 (MQ=255)
                                                                                                                                                tGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATcc  >  3:159420/1‑148 (MQ=255)
                                                                                                                                                 ggCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGa                                            <  3:874581/105‑1 (MQ=255)
                                                                                                                                                 ggCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGa                                            >  4:874581/1‑105 (MQ=255)
                                                                                                                                                    gATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAg                                                                          <  4:1389501/72‑1 (MQ=255)
                                                                                                                                                    gATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAg                                                                          >  3:1389501/1‑72 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GTTATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCC  >  NC_000913/2616458‑2616749

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