Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,616,609 G→C 8.0% M172I (ATG→ATC 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,6090GC 8.0% 345.9 / 15.8 100M172I (ATG→ATCbepAOM protein maintenance and assembly metalloprotease and chaperone, periplasmic
Reads supporting (aligned to +/- strand):  ref base G (45/47);  new base C (4/4);  total (49/51)
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.84e-01

Only 100 of 103 total aligned reads displayed.

ATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGC  >  NC_000913/2616461‑2616755
                                                                                                                                                    |                                                                                                                                                  
aTTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATg                                                                                                                                                    <  4:33209/149‑1 (MQ=255)
           gAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGc                                                                                                                                                  <  3:342782/140‑1 (MQ=255)
           gAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGc                                                                                                                                                  >  4:342782/1‑140 (MQ=255)
           gAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCCATCGCCAGtaaaa                                                                                                                                          >  1:1131428/1‑144 (MQ=255)
           gAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCCATCGCCAGtaaaa                                                                                                                                          <  2:1131428/148‑5 (MQ=255)
              aGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGc                                                                                                                                      <  3:618205/149‑1 (MQ=255)
               gTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCg                                                                                                                                     <  2:183467/149‑1 (MQ=255)
                  aaCTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCgggg                                                                                                                                  >  3:156886/1‑149 (MQ=255)
                   aCTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCCATCGCCAGtaaa                                                                                                                                           >  2:680040/1‑136 (MQ=255)
                   aCTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCCATCGCCAGtaaa                                                                                                                                           <  1:680040/139‑4 (MQ=255)
                     cggCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCgg                                                                                                                                    >  1:798801/2‑144 (MQ=255)
                     cggCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCgg                                                                                                                                    <  2:798801/143‑1 (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                                                                                                                                                  <  2:844695/104‑1 (MQ=255)
                                               tcCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGc                                                                                                                                                  >  1:844695/1‑104 (MQ=255)
                                               tcCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCggg                                                                                                     >  1:372912/1‑149 (MQ=255)
                                               tcCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCCATCGCCAGtaaaa                                                                                                                                          <  3:519203/112‑5 (MQ=255)
                                               tcCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCCATCGCCAGtaaaa                                                                                                                                          >  4:519203/1‑108 (MQ=255)
                                                    cGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAAcgc                                                                                                <  2:632434/149‑1 (MQ=255)
                                                    cGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAAcgc                                                                                                <  4:40065/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:475856/149‑1 (MQ=255)
                                                                      ggCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGtt                                                                              <  4:1398642/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                                                                                                      <  1:1005493/112‑1 (MQ=255)
                                                                                   gAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCgg                                                                                                      >  2:1005493/1‑112 (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                                                         >  3:774373/1‑149 (MQ=255)
                                                                                           gcagcgcagcGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAgg                                                         <  2:261734/149‑1 (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                                    >  1:1124185/1‑149 (MQ=255)
                                                                                                                cTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTc                                    <  1:813096/149‑1 (MQ=255)
                                                                                                                cTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTc                                    >  1:1205855/1‑149 (MQ=255)
                                                                                                                cTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTc                                    >  2:813096/1‑149 (MQ=255)
                                                                                                                cTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGACCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTc                                    >  3:662572/1‑149 (MQ=255)
                                                                                                                   ggTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCgg                                                                                                                                    <  1:623667/50‑1 (MQ=255)
                                                                                                                   ggTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCgg                                                                                                                                    >  2:623667/1‑50 (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                                                                                   >  2:1236449/1‑90 (MQ=255)
                                                                                                                            gTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGgatgat                                                                                   <  1:1236449/90‑1 (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                  >  1:856637/1‑149 (MQ=255)
                                                                                                                                  ttCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCgg                  <  1:199774/149‑1 (MQ=255)
                                                                                                                                  ttCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCgg                  >  3:919417/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)
                                                                                                                                          gaCTGGCCATCGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAAt                                                                >  4:336733/2‑95 (MQ=255)
                                                                                                                                           aCTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATc         <  2:338220/149‑1 (MQ=255)
                                                                                                                                             tGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGAtt              >  4:174928/1‑142 (MQ=255)
                                                                                                                                             tGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGAtt              <  3:174928/142‑1 (MQ=255)
                                                                                                                                             tGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATcc        >  3:159420/1‑148 (MQ=255)
                                                                                                                                             tGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATcc        <  4:159420/148‑1 (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)
                                                                                                                                                  aTGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGc  <  4:1469078/149‑1 (MQ=255)
                                                                                                                                                   tGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGa                                                        >  3:707240/1‑94 (MQ=255)
                                                                                                                                                   tGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGa                                                        <  4:707240/94‑1 (MQ=255)
                                                                                                                                                    ggCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTc                    <  4:101289/129‑1 (MQ=255)
                                                                                                                                                    ggCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTc                    >  3:101289/1‑129 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ATTCCGATAACGAAAGTCAACTGGCTTCAGTTATGGCGCACGAAATCTCCCACGTCACCCAACGTCACCTGGCGCGAGCGATGGAAGATCAGCAGCGCAGCGCGCCGCTGACCTGGGTCGGCGCGTTAGGTTCTATTTTACTGGCGATGGCCAGTCCGCAGGCGGGGATGGCGGCGCTGACCGGTACACTGGCGGGAACGCGTCAGGGGATGATCAGTTTCACCCAGCAAAATGAACAGGAAGCGGACCGCATTGGTATTCAGGTGCTGCAACGCTCGGGATTCGATCCGCAGGC  >  NC_000913/2616461‑2616755

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: