Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,990,828 G→A intergenic (‑62/‑325) hemC ← / → cyaA hydroxymethylbilane synthase/adenylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,990,8280GA94.1% 281.5 / 5.0 85intergenic (‑62/‑325)hemC/cyaAhydroxymethylbilane synthase/adenylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (3/2);  new base A (39/41);  total (42/43)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTGTTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAAT  >  NC_000913/3990686‑3990975
                                                                                                                                              |                                                                                                                                                   
cAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTgtggt                                                                                                                                              <  2:454103/150‑1 (MQ=255)
   cTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGaa                                                                                                                                           >  2:1216604/1‑150 (MQ=255)
   cTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGaa                                                                                                                                           >  2:788659/1‑150 (MQ=255)
     tGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTGTTGTGGTGaaa                                                                                                                                          <  2:704165/149‑1 (MQ=255)
     tGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAAc                                                                                                                                         <  2:432563/150‑1 (MQ=255)
        cTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGt                                                                                                                                       >  2:576393/1‑149 (MQ=255)
          ttGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTGTTGTGGTGAAACGTAg                                                                                                                                     <  2:877360/149‑1 (MQ=255)
          ttGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTa                                                                                                                                      <  2:886763/148‑1 (MQ=255)
          ttGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAg                                                                                                                                     <  1:1383167/149‑1 (MQ=255)
                  gtgtGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTgcgc                                                                                                                           <  2:921624/151‑1 (MQ=255)
                  gtgtGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTgcgc                                                                                                                           >  2:286803/1‑151 (MQ=255)
                    gtgCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCaa                                                                                                                         <  2:145781/151‑1 (MQ=255)
                      gcctgccAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCaaa                                                                                                                        <  2:477234/150‑1 (MQ=255)
                            cAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTaaa                                                                                                                 <  1:725819/151‑1 (MQ=255)
                             agagTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTaaaa                                                                                                                <  1:1159566/151‑1 (MQ=255)
                                  gCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAgg                                                                                                           <  1:490023/151‑1 (MQ=255)
                                   cAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACACCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGt                                                                                                          <  1:1011152/151‑1 (MQ=255)
                                    tcgtgcccTTTGGCGTCTGGCAGTTCATAAGACTTTCGCTATCGTTCTTGTTGCCGGCTTTATCATCCGCGGTCCTTCCTATCATGCTTACCAGAGTGATGTCAGTGTTGTGCTGAAACGTAGACGCCTGCGCCAACCGTactataag                                                                                                            >  1:685104/8‑141 (MQ=255)
                                    aaGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTc                                                                                                         >  2:1183518/1‑151 (MQ=255)
                                      gTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTg                                                                                                       <  1:336612/151‑1 (MQ=255)
                                        gggCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGc                                                                                                     <  1:955206/151‑1 (MQ=255)
                                           cTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGt                                                                                                  >  2:472169/1‑151 (MQ=255)
                                            tttGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTg                                                                                                 >  2:16832/1‑151 (MQ=255)
                                             ttGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTgg                                                                                                >  2:581602/1‑151 (MQ=255)
                                             ttGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTg                                                                                                 <  1:141917/150‑1 (MQ=255)
                                             ttGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTg                                                                                                 <  2:13818/150‑1 (MQ=255)
                                              tGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTgg                                                                                                >  1:942114/1‑150 (MQ=255)
                                              tGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTgg                                                                                                <  2:762047/150‑1 (MQ=255)
                                                 cGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATcc                                                                                            >  2:239483/1‑151 (MQ=255)
                                                 cGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATcc                                                                                            >  2:267272/1‑151 (MQ=255)
                                                  gtgtGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCt                                                                                           <  1:1323670/151‑1 (MQ=255)
                                                   tgtgGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTg                                                                                          >  2:1336076/1‑151 (MQ=255)
                                                   tgtgGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTg                                                                                          >  2:1055439/1‑151 (MQ=255)
                                                     tgGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGa                                                                                         >  1:1011818/1‑150 (MQ=255)
                                                      ggCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACa                                                                                       <  1:223993/151‑1 (MQ=255)
                                                        cAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAgg                                                                                     >  1:437552/1‑151 (MQ=255)
                                                            tCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGtt                                                                                 >  2:1201516/1‑151 (MQ=255)
                                                                   aCATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGcc                                                                          <  2:1202535/151‑1 (MQ=255)
                                                                    cATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCg                                                                         >  1:693496/1‑151 (MQ=255)
                                                                     aTTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGt                                                                        >  1:473833/1‑151 (MQ=255)
                                                                       tGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGtt                                                                       <  2:321237/150‑1 (MQ=255)
                                                                        gTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGt                                                                     <  1:696206/151‑1 (MQ=255)
                                                                             aCATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGtaata                                                                 >  1:561789/1‑150 (MQ=255)
                                                                                 gCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGa                                                             >  1:243656/1‑150 (MQ=255)
                                                                                   ttGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTGTTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAAtt                                                          >  2:1299221/1‑151 (MQ=255)
                                                                                   ttGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGATAGGCGTTTCACGCCGTTGTAATAAGGAAtt                                                          >  2:782821/1‑151 (MQ=255)
                                                                                    tGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAAtt                                                          <  1:404030/150‑1 (MQ=255)
                                                                                     gTTACCGTCATTATCATCCGTGGTCCATCCTAACATACTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAAttt                                                         >  1:227593/1‑150 (MQ=255)
                                                                                         ccGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTAc                                                       >  1:957599/1‑148 (MQ=255)
                                                                                          cGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACag                                                     >  1:602202/1‑149 (MQ=255)
                                                                                           gTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAgaga                                                  <  2:1157103/151‑1 (MQ=255)
                                                                                             cATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAAt                                                >  1:1060055/1‑151 (MQ=255)
                                                                                             cATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAAt                                                <  1:859319/151‑1 (MQ=255)
                                                                                               ttATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATa                                               >  2:572586/1‑150 (MQ=255)
                                                                                                 atcatcCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATaaa                                             <  2:1152530/150‑1 (MQ=255)
                                                                                                 atcatcCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAAc                                            <  2:737957/151‑1 (MQ=255)
                                                                                                   catcCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATaaa                                             >  2:1413326/1‑148 (MQ=255)
                                                                                                   catcCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACgg                                          >  1:455196/1‑151 (MQ=255)
                                                                                                      ccGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTg                                        <  2:953523/150‑1 (MQ=255)
                                                                                                      ccGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGc                                       <  2:1090233/151‑1 (MQ=255)
                                                                                                      ccGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGc                                       <  1:414730/151‑1 (MQ=255)
                                                                                                      ccGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGc                                       >  1:112831/1‑151 (MQ=255)
                                                                                                       cGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCt                                      >  2:24052/1‑151 (MQ=255)
                                                                                                        gTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCt                                      >  1:946337/1‑150 (MQ=255)
                                                                                                        gTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTa                                     <  1:415896/151‑1 (MQ=255)
                                                                                                          ggTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTac                                    >  1:398572/1‑150 (MQ=255)
                                                                                                          ggTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTac                                    <  2:36683/150‑1 (MQ=255)
                                                                                                              cATCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTCATAAGGAATTTACAGAGAATAAACGGTCCTACACtt                                >  2:1158709/1‑150 (MQ=255)
                                                                                                               aTCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTtgt                              <  2:1384230/151‑1 (MQ=255)
                                                                                                                tCCTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTtgt                              >  2:1049543/1‑150 (MQ=255)
                                                                                                                 ccTAACAGCCTTGCCAGAGTGATGTCAGTATTGTGGTGGAACGTAGCCGCCTGCTCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTtgga                             <  1:915177/150‑3 (MQ=255)
                                                                                                                  cTAACATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTtgta                             <  2:916889/149‑1 (MQ=255)
                                                                                                                    aaCATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTtgtatg                           <  1:788622/149‑1 (MQ=255)
                                                                                                                     aCATCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAg                        <  2:318514/151‑1 (MQ=255)
                                                                                                                        tCCTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAgcgc                     >  2:552317/1‑151 (MQ=255)
                                                                                                                         ccTTGCCAGAGTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAgcgc                     >  1:506908/1‑150 (MQ=255)
                                                                                                                                agagTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCAtctttct             <  2:806632/151‑1 (MQ=255)
                                                                                                                                 gagTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCAtctttctt            <  2:885759/151‑1 (MQ=255)
                                                                                                                                  agTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTATCGCATctttcttt           <  2:669502/151‑1 (MQ=255)
                                                                                                                                  agTGATGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCAtctttct             <  2:1416454/149‑1 (MQ=255)
                                                                                                                                      aTGTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACgg       >  2:1278230/1‑151 (MQ=255)
                                                                                                                                        gTCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTc     <  1:1296689/151‑1 (MQ=255)
                                                                                                                                         tCAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTc     >  2:245420/1‑150 (MQ=255)
                                                                                                                                          cAGTATTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCa    >  1:964935/1‑150 (MQ=255)
                                                                                                                                           aGTGTTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATTAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAAt  >  2:1354152/1‑151 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
CAACTTGTCTTTGACATAGTGTGCCTGCCAGAGTGCAAGTGGGCTTTGGCGTGTGGCAATTCTTAAAACATTGTCTAACATGCTTGTTACCGTCATTATCATCCGTGGTCCATCCTAACATCCTTGCCAGAGTGATGTCAGTGTTGTGGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAAT  >  NC_000913/3990686‑3990975

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