Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,445,621 G→C 7.6% W857C (TGG→TGC 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,6210GC 7.6% 413.8 / 16.0 116W857C (TGG→TGCydbHputative membrane‑anchored protein, function unknown
Reads supporting (aligned to +/- strand):  ref base G (57/50);  new base C (4/5);  total (61/55)
Fisher's exact test for biased strand distribution p-value = 7.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 118 total aligned reads displayed.

GTTTATTAACCATGCAGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTGAAGTCGCTGCACCT  >  NC_000913/1445479‑1445764
                                                                                                                                              |                                                                                                                                               
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)
                    aaCATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCa                                                                                                                                         <  4:298032/131‑1 (MQ=255)
                    aaCATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGCCATGCCTGGa                                                                                                                                            <  3:13282/128‑1 (MQ=255)
                    aaCATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGCCATGCCTGGa                                                                                                                                            >  4:13282/1‑128 (MQ=255)
                        ttACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAAt                                                                                                                   >  2:1209691/1‑149 (MQ=255)
                        ttACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAAt                                                                                                                   <  1:512041/149‑1 (MQ=255)
                             ggTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGa                                                                                                                                            >  3:1154177/1‑119 (MQ=255)
                             ggTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGa                                                                                                                                            <  4:1154177/119‑1 (MQ=255)
                              gTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGa                                                                                                                                       >  4:1114376/1‑123 (MQ=255)
                              gTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGt                                                                                                             >  1:307218/1‑149 (MQ=255)
                               tACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTc                                                                                                            >  1:924505/1‑149 (MQ=255)
                                   aGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCgg                                                                                                                        >  3:720319/1‑133 (MQ=255)
                                    gTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGa                                                                                                                 <  4:1319629/139‑1 (MQ=255)
                                    gTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGa                                                                                                                 >  3:1319629/1‑139 (MQ=255)
                                          ttGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGaaaa                                                                                                 >  4:1156826/1‑149 (MQ=255)
                                            gATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAaga                                                                                               <  2:924505/149‑1 (MQ=255)
                                            gATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAaga                                                                                               <  2:1077904/149‑1 (MQ=255)
                                              tGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAgt                                                                                             >  1:320748/1‑149 (MQ=255)
                                               ggTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGAc                                                                                                                <  3:1464668/129‑1 (MQ=255)
                                               ggTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGAc                                                                                                                >  4:1464668/1‑129 (MQ=255)
                                                gTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGa                                                                                                                 <  1:300711/127‑1 (MQ=255)
                                                gTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGa                                                                                                                 >  2:300711/1‑127 (MQ=255)
                                                gTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACt                                                                                                               <  1:1302129/129‑1 (MQ=255)
                                                gTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACt                                                                                                               >  2:1302129/1‑129 (MQ=255)
                                                 tAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCaa                                                                                                       >  2:222047/1‑136 (MQ=255)
                                                 tAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCaa                                                                                                       <  1:222047/136‑1 (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                                                                                                         <  1:1388257/126‑1 (MQ=255)
                                                         gtacggtaAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGc                                                                                                         >  2:1388257/1‑126 (MQ=255)
                                                            cggtaAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACaa                                                                               <  1:588715/149‑1 (MQ=255)
                                                            cggtaAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACaa                                                                               <  3:898360/149‑1 (MQ=255)
                                                              gtaAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGa                                                                                         <  2:623009/137‑1 (MQ=255)
                                                              gtaAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGa                                                                                         >  1:623009/1‑137 (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)
                                                                            ccATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACt                                                                                                               >  1:340044/1‑101 (MQ=255)
                                                                            ccATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAaga                                                                                               >  1:775071/1‑117 (MQ=255)
                                                                            ccATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAaga                                                                                               <  2:775071/117‑1 (MQ=255)
                                                                               tcatGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGc                                                            >  3:1228357/1‑149 (MQ=255)
                                                                                   gAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCgg                                                                                                                        <  3:1459521/85‑1 (MQ=255)
                                                                                   gAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCgg                                                                                                                        >  4:1459521/1‑85 (MQ=255)
                                                                                   gAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTg                                                        >  1:361550/1‑149 (MQ=255)
                                                                                    aaGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGa                                                                                                                                            >  4:775646/1‑64 (MQ=255)
                                                                                    aaGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGa                                                                                                                                            <  3:775646/64‑1 (MQ=255)
                                                                                      gagaATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGa                                                                                         >  3:76423/1‑113 (MQ=255)
                                                                                      gagaATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGa                                                                                         <  4:76423/113‑1 (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)
                                                                                                    gCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCt                                       >  4:1322018/1‑149 (MQ=255)
                                                                                                       gtgGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGt                                    <  2:361550/149‑1 (MQ=255)
                                                                                                        tgGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGa                                                                                                                 <  3:693881/71‑1 (MQ=255)
                                                                                                        tgGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGa                                                                                                                 >  4:693881/1‑71 (MQ=255)
                                                                                                                tttACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTAGTTGGCTGTAc                           <  3:1156826/149‑1 (MQ=255)
                                                                                                                   aCGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTAt                                          <  3:121225/131‑1 (MQ=255)
                                                                                                                       tttGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAAt                                                                        >  1:1040367/1‑97 (MQ=255)
                                                                                                                       tttGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGc                    >  4:475714/1‑149 (MQ=255)
                                                                                                                       tttGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAACGAAAAt                                                                        <  2:1040367/97‑1 (MQ=255)
                                                                                                                           gCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTg                >  1:516822/1‑149 (MQ=255)
                                                                                                                           gCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTg                >  1:156101/1‑149 (MQ=255)
                                                                                                                             gATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCTAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCtt                                              <  3:1372268/117‑1 (MQ=255)
                                                                                                                             gATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAAt                                                                        <  2:1324450/91‑1 (MQ=255)
                                                                                                                             gATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGGGTGAGGAAAAACAATGAAAAt                                                                        >  1:1324450/1‑91 (MQ=255)
                                                                                                                                 tgctcCAGCCATGCCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTGaa              >  1:423839/3‑145 (MQ=255)
                                                                                                                                 tgatcCAGCCATGCCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTGaa              <  2:423839/142‑1 (MQ=255)
                                                                                                                                    tcCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAag                                                                                                >  4:131358/1‑60 (MQ=255)
                                                                                                                                    tcCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCt                                                              >  2:238521/1‑94 (MQ=255)
                                                                                                                                    tcCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCt                       <  3:309496/133‑1 (MQ=255)
                                                                                                                                    tcCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTGa               <  4:793748/141‑1 (MQ=255)
                                                                                                                                       aGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAgtg                                                                                            <  4:725606/61‑1 (MQ=255)
                                                                                                                                       aGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAgtg                                                                                            >  3:725606/1‑61 (MQ=255)
                                                                                                                                       aGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTa                                           >  2:948974/1‑110 (MQ=255)
                                                                                                                                       aGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTa                                           <  1:948974/110‑1 (MQ=255)
                                                                                                                                         gCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGt                             >  2:160005/1‑122 (MQ=255)
                                                                                                                                         gCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGt                             <  1:160005/122‑1 (MQ=255)
                                                                                                                                         gCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTGAAGTCGCTGCACCt  >  4:463533/1‑149 (MQ=255)
                                                                                                                                          cATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGt                                                    >  4:720982/1‑98 (MQ=255)
                                                                                                                                          cATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGt                                                    <  3:720982/98‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                               
GTTTATTAACCATGCAGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTGAAGTCGCTGCACCT  >  NC_000913/1445479‑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: