Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,988,049 C→T 8.4% intergenic (‑90/‑173) azuC ← / → yecR acid‑inducible small membrane‑associated protein/lipoprotein, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,988,0490CT 8.4% 249.3 / 9.8 83intergenic (‑90/‑173)azuC/yecRacid‑inducible small membrane‑associated protein/lipoprotein, function unknown
Reads supporting (aligned to +/- strand):  ref base C (39/37);  new base T (2/5);  total (41/42)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.43e-01

TTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAGTGC  >  NC_000913/1987903‑1988196
                                                                                                                                                   |                                                                                                                                                   
tttGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTAATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACa                                                                                                                                                    >  4:108186/1‑148 (MQ=255)
     aCGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATcac                                                                                                                                              <  3:3898/149‑1 (MQ=255)
      cGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATcaca                                                                                                                                             <  2:139141/149‑1 (MQ=255)
       gTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACAt                                                                                                                                            <  8:130100/149‑1 (MQ=255)
        tCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTTGCGCAGTTTCATCGCGGCACTCCGACAATTTGTTATACAGGAGTTATTGTCTTTGCGCTTATAATATGACAAACATCACAAAAAACAAACTTTATGTGATACAAATcaca                                                                                                                                             >  4:79187/1‑148 (MQ=255)
        tCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATa                                                                                                                                           <  7:119253/149‑1 (MQ=255)
           tGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTACGACCATTTGTTATACAGGTGTTATTGTATTTGCCCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTTATACAAATCACATaaa                                                                                                                                         >  1:186168/1‑148 (MQ=255)
            gCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATTTGTATCACATaaag                                                                                                                                        <  7:54235/148‑2 (MQ=255)
            gCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATTTGTATCACATaa                                                                                                                                          <  7:283812/146‑1 (MQ=255)
            gCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATaa                                                                                                                                          <  1:328622/146‑1 (MQ=255)
             cAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATTTGTATCACATaaagat                                                                                                                                      <  1:193536/149‑4 (MQ=255)
             cAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATa                                                                                                                                       <  3:135646/148‑1 (MQ=255)
                tAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGCCAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATAcccc                                                                                                                                   >  3:269809/1‑149 (MQ=255)
                   ttATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCtt                                                                                                                                 >  7:255518/1‑148 (MQ=255)
                     aTTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTaa                                                                                                                              >  7:22502/1‑149 (MQ=255)
                     aTTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTa                                                                                                                               <  5:189995/148‑1 (MQ=255)
                      ttGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTaa                                                                                                                              <  7:6718/148‑1 (MQ=255)
                      ttGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAAt                                                                                                                             <  1:77741/149‑1 (MQ=255)
                      ttGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAAt                                                                                                                             >  5:136285/1‑149 (MQ=255)
                        gAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTTTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTCTATGTGATACAAATCACATAAATACCCCTTTAATg                                                                                                                            <  6:102582/148‑1 (MQ=255)
                         aaCATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGt                                                                                                                           >  8:96027/1‑148 (MQ=255)
                         aaCATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGt                                                                                                                           >  7:20718/1‑148 (MQ=255)
                                ttttCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATaaaaa                                                                                                                   >  8:37016/1‑149 (MQ=255)
                                 tttCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAAt                                                                                                                  >  3:257088/1‑149 (MQ=255)
                                  ttgcggATGTTGCGCAGTGTCCTCGCGGCACTCCGACCTTTTGTTCTCCAGGTGTTAGTGGCTGTGCGCTTATAATATGACAACCATCACAAAAATCACTCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATaaaaa                                                                                                                   <  2:217407/143‑1 (MQ=255)
                                   tCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACACCCATCACAAAAATCAATCTTTATGTGATACAAATCACATAACTACCCCTTTAATGTTATAAAAATGa                                                                                                                >  8:131915/1‑149 (MQ=255)
                                     aGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGAt                                                                                                               <  7:96027/148‑1 (MQ=255)
                                     aGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATa                                                                                                              <  5:7048/149‑1 (MQ=255)
                                      ggATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATaa                                                                                                             <  8:22502/149‑1 (MQ=255)
                                      atgtttttGCGAAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGAt                                                                                                               <  1:102109/144‑1 (MQ=255)
                                          ttttGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATc                                                                                                           >  7:37513/1‑148 (MQ=255)
                                              gcgcAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCaaaa                                                                                                       <  7:9847/148‑1 (MQ=255)
                                                gcAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATTTGTATCACATaa                                                                                                                                          >  5:168058/1‑111 (MQ=255)
                                                gcAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATTTGTATCACATaa                                                                                                                                          <  6:168058/111‑1 (MQ=255)
                                                    tttCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAACATGATAATCAACAACCAGcc                                                                                                >  2:295311/1‑149 (MQ=255)
                                                      tCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCGATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGccca                                                                                              >  8:9009/1‑148 (MQ=255)
                                                      tCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGccc                                                                                               <  2:96076/148‑1 (MQ=255)
                                                           gcgGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAAGCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTAt                                                                                         <  6:336572/149‑1 (MQ=255)
                                                               cACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCAAATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCt                                                                                     <  3:108186/149‑1 (MQ=255)
                                                                 cTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGa                                                                                   <  2:79806/149‑1 (MQ=255)
                                                                 cTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGa                                                                                   >  6:224122/1‑149 (MQ=255)
                                                                  tCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGa                                                                                   >  8:212490/1‑148 (MQ=255)
                                                                         aTTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTa                                                                             >  4:91158/1‑147 (MQ=255)
                                                                         aTTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACa                                                                           <  8:20718/149‑1 (MQ=255)
                                                                          tttGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTa                                                                             >  6:334701/1‑146 (MQ=255)
                                                                          tttGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGAAACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTa                                                                             >  8:152471/1‑146 (MQ=255)
                                                                                ataCAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCATATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGc                                                                     <  1:55876/148‑1 (MQ=255)
                                                                                    aGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGaa                                                                <  8:37513/149‑1 (MQ=255)
                                                                                     ggTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGaaa                                                               >  2:310951/1‑149 (MQ=255)
                                                                                      gtgtTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGaaa                                                               >  4:32340/1‑148 (MQ=255)
                                                                                       tgtTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATa                                                             <  3:56123/149‑1 (MQ=255)
                                                                                        gtTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGa                                                                                   >  3:205670/1‑126 (MQ=255)
                                                                                        gtTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGa                                                                                   <  4:205670/126‑1 (MQ=255)
                                                                                         ttATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGccccc                                                                                             >  3:57705/1‑115 (MQ=255)
                                                                                         ttATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGccccc                                                                                             <  4:57705/115‑1 (MQ=255)
                                                                                         ttATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCGCTATTTCTGTCACCTACAGATGGCAAGAAATAgc                                                           >  8:126848/1‑149 (MQ=255)
                                                                                            ttGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAgcgac                                                        >  5:236018/1‑147 (MQ=255)
                                                                                                   tGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAg                                                  <  6:148068/148‑1 (MQ=255)
                                                                                                         tataATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCtt                                           >  5:322420/1‑149 (MQ=255)
                                                                                                                   cAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGCCAAGAAAAAGAGCCTGCCCGGCGTCTTTTCCGGCCAt                                 >  8:162997/1‑149 (MQ=255)
                                                                                                                    aaCCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCAtt                                <  4:332404/149‑1 (MQ=255)
                                                                                                                      ccATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTa                                                                             <  2:274470/102‑1 (MQ=255)
                                                                                                                      ccATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTa                                                                             >  1:274470/1‑102 (MQ=255)
                                                                                                                      ccATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTga                              >  7:198890/1‑148 (MQ=255)
                                                                                                                       gaTCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTc                             <  3:79187/148‑1 (MQ=255)
                                                                                                                       cATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTc                             >  8:53371/1‑149 (MQ=255)
                                                                                                                             aaaaaTCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCgcagca                       <  1:153588/149‑1 (MQ=255)
                                                                                                                                aatcaatcTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTg                    <  6:200855/149‑1 (MQ=255)
                                                                                                                                   caatcTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTAAAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCCGCACTGTaa                 >  5:51608/1‑149 (MQ=255)
                                                                                                                                    aatcTTTATGTTATACAAATCACATAAATACGACTTTAATGTTATAAAAATGATAATCAAAAAACATCCCACCTATTTCTGACACCTACAGATGGTAAGAAATAGCTCCTGCCAGGCGTCTTTTCATTCAATTGTCTCAGCACTGTAAc                <  8:86148/149‑1 (MQ=255)
                                                                                                                                    aatcTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTaa                 <  1:295311/148‑1 (MQ=255)
                                                                                                                                      tcTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCaa                                                                   >  7:316262/1‑96 (MQ=255)
                                                                                                                                      tcTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCaa                                                                   <  8:316262/96‑1 (MQ=255)
                                                                                                                                      tcTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGATATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAAggc              >  5:334002/1‑149 (MQ=255)
                                                                                                                                        tttATGTGATTTGTATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCgc                           <  2:105190/134‑1 (MQ=255)
                                                                                                                                        tttATGTGATTTGTATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCgc                           >  1:105190/1‑134 (MQ=255)
                                                                                                                                            tgtgATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTaaa         <  6:322420/148‑1 (MQ=255)
                                                                                                                                            tgtgATACAAATCACATAAATACCACTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTaaaa        <  4:305506/149‑1 (MQ=255)
                                                                                                                                             gtgATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTaaaa        >  1:25837/1‑148 (MQ=255)
                                                                                                                                               gATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAg     >  3:334166/1‑149 (MQ=255)
                                                                                                                                                aTACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAGt    >  1:226786/1‑149 (MQ=255)
                                                                                                                                                aTACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAGt    >  5:14096/1‑149 (MQ=255)
                                                                                                                                                  aCAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAGTGc  >  3:51067/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTATAATATGACAACCATCACAAAAATCAATCTTTATGTGATACAAATCACATAAATACCCCTTTAATGTTATAAAAATGATAATCAAAAAACAGCCCCCCTATTTCTGACACCTACAGATGGCAAGAAATAGCGCCTGCCAGGCGTCTTTTCCGGCCATTGTCGCAGCACTGTAACGCGTAAAATAGTGC  >  NC_000913/1987903‑1988196

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