Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,445,616 G→C 7.9% A856P (GCA→CCA)  ydbH → putative membrane‑anchored protein, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,445,6160GC 7.9% 400.2 / 19.5 114A856P (GCA→CCA) ydbHputative membrane‑anchored protein, function unknown
Reads supporting (aligned to +/- strand):  ref base G (55/50);  new base C (4/5);  total (59/55)
Fisher's exact test for biased strand distribution p-value = 7.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01

Only 100 of 117 total aligned reads displayed.

ATAATCTCGGTTTATTAACCATGCAGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTGAAGTCGCTGCACCT  >  NC_000913/1445470‑1445764
                                                                                                                                                  |                                                                                                                                                    
aTAATCTCGGTTTATTAACCATGCAGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGGAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCa                                                                                                                                                    <  2:903436/149‑1 (MQ=255)
         gTTTATTAACCATGCAGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATgg                                                                                                                                                 >  2:947400/1‑143 (MQ=255)
         gTTTATTAACCATGCAGGCCAACACTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATATTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATgg                                                                                                                                                 <  1:947400/143‑1 (MQ=255)
          tttATTAACCATGCAGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGCCATGCCTGGAga                                                                                                                                          <  3:1210852/149‑2 (MQ=255)
              ttaACCATGCAGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGaa                                                                                                                                      >  2:512041/1‑149 (MQ=255)
                aaCCATGCAGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGc                                                                                                                                          <  4:138482/143‑1 (MQ=255)
                aaCCATGCAGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGc                                                                                                                                          >  3:138482/1‑143 (MQ=255)
                        aGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTg                                                                                                                            >  1:1077904/1‑149 (MQ=255)
                         ggCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGc                                                                                                                           >  4:1307903/1‑149 (MQ=255)
                            cAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGCCATGCCTGGAg                                                                                                                                           >  4:373526/1‑130 (MQ=255)
                            cAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGCCATGCCTGGAg                                                                                                                                           <  3:373526/130‑1 (MQ=255)
                             aaCATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCa                                                                                                                                         >  3:298032/1‑131 (MQ=255)
                             aaCATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGCCATGCCTGGa                                                                                                                                            >  4:13282/1‑128 (MQ=255)
                                 ttACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAAt                                                                                                                   <  1:512041/149‑1 (MQ=255)
                                      ggTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGa                                                                                                                                            >  3:1154177/1‑119 (MQ=255)
                                       gTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGa                                                                                                                                       <  3:1114376/123‑1 (MQ=255)
                                       gTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGa                                                                                                                                       >  4:1114376/1‑123 (MQ=255)
                                       gTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGt                                                                                                             >  1:307218/1‑149 (MQ=255)
                                            aGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCgg                                                                                                                        <  4:720319/133‑1 (MQ=255)
                                            aGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCgg                                                                                                                        >  3:720319/1‑133 (MQ=255)
                                             gTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGa                                                                                                                 >  3:1319629/1‑139 (MQ=255)
                                             gTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGa                                                                                                                 <  4:1319629/139‑1 (MQ=255)
                                             gTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCaa                                                                                                       <  1:81624/149‑1 (MQ=255)
                                                   ttGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGaaaa                                                                                                 >  4:1156826/1‑149 (MQ=255)
                                                     gATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAaga                                                                                               <  2:1077904/149‑1 (MQ=255)
                                                       tGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAgt                                                                                             >  1:320748/1‑149 (MQ=255)
                                                        ggTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGAc                                                                                                                >  4:1464668/1‑129 (MQ=255)
                                                        ggTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGAc                                                                                                                <  3:1464668/129‑1 (MQ=255)
                                                         gTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGa                                                                                                                 >  2:300711/1‑127 (MQ=255)
                                                         gTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGa                                                                                                                 <  1:300711/127‑1 (MQ=255)
                                                         gTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACt                                                                                                               >  2:1302129/1‑129 (MQ=255)
                                                         gTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACt                                                                                                               <  1:1302129/129‑1 (MQ=255)
                                                          tAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCaa                                                                                                       <  1:222047/136‑1 (MQ=255)
                                                          tAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCaa                                                                                                       >  2:222047/1‑136 (MQ=255)
                                                          tAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACACACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGtgtg                                                                                          >  3:475772/1‑149 (MQ=255)
                                                           aaaaGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACACACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGa                                                                                         <  4:475772/149‑1 (MQ=255)
                                                            aaaGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAg                                                                                        <  2:307218/149‑1 (MQ=255)
                                                            aaaGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAg                                                                                        >  2:588715/1‑149 (MQ=255)
                                                               gCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGaa                                                                                     <  3:1307903/149‑1 (MQ=255)
                                                                  gtacggtaAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGc                                                                                                         >  2:1388257/1‑126 (MQ=255)
                                                                  gtacggtaAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGc                                                                                                         <  1:1388257/126‑1 (MQ=255)
                                                                     cggtaAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACaa                                                                               <  3:898360/149‑1 (MQ=255)
                                                                     cggtaAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACaa                                                                               <  1:588715/149‑1 (MQ=255)
                                                                       gtaAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGa                                                                                         >  1:623009/1‑137 (MQ=255)
                                                                       gtaAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGa                                                                                         <  2:623009/137‑1 (MQ=255)
                                                                       gtaAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATg                                                                             >  4:227876/1‑149 (MQ=255)
                                                                            ccTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAAt                                                                        <  2:320748/149‑1 (MQ=255)
                                                                            ccTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAAt                                                                        <  1:1209691/149‑1 (MQ=255)
                                                                             cTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAAtt                                                                       >  2:1138349/1‑149 (MQ=255)
                                                                              ttAATTACCATTATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGa                                                                                                                 >  4:1004836/1‑106 (MQ=255)
                                                                              ttAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGa                                                                                                                 <  3:1004836/106‑1 (MQ=255)
                                                                                aaTTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTa                                                                    <  1:1138349/149‑1 (MQ=255)
                                                                                  ttACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAaga                                                                                               <  4:469648/120‑1 (MQ=255)
                                                                                  ttACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAaga                                                                                               >  3:469648/1‑120 (MQ=255)
                                                                                  ttACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGCCATGCCTGGAgatta                                                                                                                                       >  1:207149/1‑76 (MQ=39)
                                                                                  ttACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGCCATGCCTGGAgatta                                                                                                                                       <  2:207149/80‑5 (MQ=39)
                                                                                     ccATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACt                                                                                                               >  1:340044/1‑101 (MQ=255)
                                                                                     ccATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACt                                                                                                               <  2:340044/101‑1 (MQ=255)
                                                                                     ccATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAaga                                                                                               <  2:775071/117‑1 (MQ=255)
                                                                                     ccATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAaga                                                                                               >  1:775071/1‑117 (MQ=255)
                                                                                        tcatGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGc                                                            >  3:1228357/1‑149 (MQ=255)
                                                                                            gAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCgg                                                                                                                        >  4:1459521/1‑85 (MQ=255)
                                                                                            gAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCgg                                                                                                                        <  3:1459521/85‑1 (MQ=255)
                                                                                            gAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTg                                                        >  1:361550/1‑149 (MQ=255)
                                                                                            gAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTg                                                        <  3:227876/149‑1 (MQ=255)
                                                                                             aaGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGa                                                                                                                                            <  3:775646/64‑1 (MQ=255)
                                                                                             aaGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGa                                                                                                                                            >  4:775646/1‑64 (MQ=255)
                                                                                               gagaATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGa                                                                                         <  4:76423/113‑1 (MQ=255)
                                                                                               gagaATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGa                                                                                         >  3:76423/1‑113 (MQ=255)
                                                                                                  aaTATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGAAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCa                                                  >  1:1123210/1‑149 (MQ=255)
                                                                                                   atatTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCAt                                                 >  2:415720/1‑149 (MQ=255)
                                                                                                         ttACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTa                                           >  3:42604/1‑149 (MQ=255)
                                                                                                         ttACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTa                                           <  4:1228357/149‑1 (MQ=255)
                                                                                                                gtgGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGt                                    <  2:361550/149‑1 (MQ=255)
                                                                                                                 tgGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGa                                                                                                                 <  3:693881/71‑1 (MQ=255)
                                                                                                                            aCGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTAt                                          >  4:121225/1‑131 (MQ=255)
                                                                                                                            aCGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTAt                                          <  3:121225/131‑1 (MQ=255)
                                                                                                                                tttGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAAt                                                                        >  1:1040367/1‑97 (MQ=255)
                                                                                                                                tttGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGc                    >  4:475714/1‑149 (MQ=255)
                                                                                                                                tttGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGc                    >  2:752368/1‑149 (MQ=255)
                                                                                                                                tttGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAACGAAAAt                                                                        <  2:1040367/97‑1 (MQ=255)
                                                                                                                                   ggCGATAATCTCCAGGCATGGCTGGAGCAGATCGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCTTCTTTTATGCTGGTTGGCTGTACGCCTAGCAtt                 <  1:415720/149‑1 (MQ=255)
                                                                                                                                    gCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTg                >  1:516822/1‑149 (MQ=255)
                                                                                                                                    gCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTg                >  1:156101/1‑149 (MQ=255)
                                                                                                                                      gATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCTAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCtt                                              <  3:1372268/117‑1 (MQ=255)
                                                                                                                                      gATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCTAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCtt                                              >  4:1372268/1‑117 (MQ=255)
                                                                                                                                      gATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAAt                                                                        <  2:1324450/91‑1 (MQ=255)
                                                                                                                                      gATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTGaa              >  2:1137156/1‑149 (MQ=255)
                                                                                                                                      gATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGGGTGAGGAAAAACAATGAAAAt                                                                        >  1:1324450/1‑91 (MQ=255)
                                                                                                                                             tcCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAag                                                                                                <  3:131358/60‑1 (MQ=255)
                                                                                                                                             tcCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAag                                                                                                >  4:131358/1‑60 (MQ=255)
                                                                                                                                             tcCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCt                       >  4:309496/1‑133 (MQ=255)
                                                                                                                                             tcCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTGa               >  3:793748/1‑141 (MQ=255)
                                                                                                                                             tcCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTGa               <  4:793748/141‑1 (MQ=255)
                                                                                                                                                aGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAgtg                                                                                            <  4:725606/61‑1 (MQ=255)
                                                                                                                                                aGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTa                                           <  1:948974/110‑1 (MQ=255)
                                                                                                                                                aGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTa                                           >  2:948974/1‑110 (MQ=255)
                                                                                                                                                  gCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGt                             <  1:160005/122‑1 (MQ=255)
                                                                                                                                                  gCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGt                             >  2:160005/1‑122 (MQ=255)
                                                                                                                                                  gCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTGAAGTCGCTGCACCt  >  4:463533/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
ATAATCTCGGTTTATTAACCATGCAGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTGAAGTCGCTGCACCT  >  NC_000913/1445470‑1445764

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