Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,445,616 G→C 6.3% 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 6.3% 452.6 / 17.2 127A856P (GCA→CCA) ydbHputative membrane‑anchored protein, function unknown
Reads supporting (aligned to +/- strand):  ref base G (63/56);  new base C (4/4);  total (67/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.26e-01

Only 100 of 127 total aligned reads displayed.

GATAATCTCGGTTTATTAACCATGCAGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTGAAGTCGCTGCAC  >  NC_000913/1445469‑1445762
                                                                                                                                                   |                                                                                                                                                  
gATAATCTCGGTTTATTAACCATGCAGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGc                                                                                                                                                   >  3:894552/1‑149 (MQ=255)
 aTAATCTCGGTTTATTAACCATGCAGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCa                                                                                                                                                  >  4:820376/1‑149 (MQ=255)
  tAATCTCGGTTTATTAACCATGCAGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCAt                                                                                                                                                 <  4:894552/149‑1 (MQ=255)
          gTTTATTAACCATGCAGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGCCATGCCTGGAg                                                                                                                                         >  2:722747/1‑149 (MQ=255)
           tttATTAACCATGCAGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGCCATGCCTGGAga                                                                                                                                        <  1:722747/149‑2 (MQ=255)
               ttaACCATGCAGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTTGAGCAGaa                                                                                                                                    >  1:909570/1‑149 (MQ=255)
                   ccATGCAGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTg                                                                                                                                            <  1:565673/137‑1 (MQ=255)
                   ccATGCAGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTg                                                                                                                                            >  2:565673/1‑137 (MQ=255)
                   ccATGCAGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCa                                                                                                                                >  1:645217/1‑149 (MQ=255)
                     aTGCAGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGCCATGCCTGGAga                                                                                                                                        >  1:102831/1‑138 (MQ=255)
                     aTGCAGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGCCATGCCTGGAga                                                                                                                                        <  2:102831/139‑2 (MQ=255)
                      tGCAGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGt                                                                                                                             >  1:1032997/1‑149 (MQ=255)
                            ccAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCg                                                                                                                       >  1:202014/1‑149 (MQ=255)
                             cAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCgg                                                                                                                      <  3:800386/149‑1 (MQ=255)
                              aaCATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCggg                                                                                                                     >  3:83318/1‑149 (MQ=255)
                                cATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGaa                                                                                                                   >  2:76951/1‑149 (MQ=255)
                                  ttACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAAt                                                                                                                 <  1:809022/149‑1 (MQ=255)
                                     cAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGAc                                                                                                              <  2:202014/149‑1 (MQ=255)
                                      aGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACt                                                                                                             >  4:1255412/1‑149 (MQ=255)
                                        gTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGt                                                                                                           >  3:1561228/1‑149 (MQ=255)
                                             aGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGaaa                                                                                                                  >  4:1147902/1‑137 (MQ=255)
                                             aGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGaaa                                                                                                                  <  3:1147902/137‑1 (MQ=255)
                                               tCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAg                                                                                                    >  4:159814/1‑149 (MQ=255)
                                               tCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAg                                                                                                    >  1:75080/1‑149 (MQ=255)
                                               tCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGCCATGCCTGGAgatta                                                                                                                                     >  3:1828475/1‑112 (MQ=255)
                                                      gATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAaga                                                                                             <  4:56983/149‑1 (MQ=255)
                                                       aTGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTcc                                                                                                         >  3:1410712/1‑136 (MQ=255)
                                                       aTGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTcc                                                                                                         <  4:1410712/136‑1 (MQ=255)
                                                       aTGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAagag                                                                                            >  2:1439738/1‑149 (MQ=255)
                                                        tGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTcc                                                                                                         >  3:321697/1‑135 (MQ=255)
                                                        tGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTcc                                                                                                         <  4:321697/135‑1 (MQ=255)
                                                        tGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAgt                                                                                           <  1:1439738/149‑1 (MQ=255)
                                                        tGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAgt                                                                                           <  2:645217/149‑1 (MQ=255)
                                                         ggTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAgtg                                                                                          <  4:1696564/149‑1 (MQ=255)
                                                          gTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAgtgt                                                                                         >  2:684370/1‑149 (MQ=255)
                                                           tAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAg                                                                                                    >  2:1584736/1‑137 (MQ=255)
                                                           tAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAg                                                                                                    <  1:1584736/137‑1 (MQ=255)
                                                            aaaaGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGc                                                                                                       >  3:825483/1‑133 (MQ=255)
                                                            aaaaGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGc                                                                                                       <  4:825483/133‑1 (MQ=255)
                                                             aaaGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAg                                                                                      <  2:75080/149‑1 (MQ=255)
                                                                  ggtacggtaAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGaaaa                                                                                 >  3:819011/1‑149 (MQ=255)
                                                                   gtacggtaAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGaaaaa                                                                                >  4:1031944/1‑149 (MQ=255)
                                                                      cggtaAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACaa                                                                             <  1:1108659/149‑1 (MQ=255)
                                                                       ggtaAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGaaaaa                                                                                <  2:1725322/145‑1 (MQ=255)
                                                                       ggtaAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAAt                                                                            <  4:1561228/149‑1 (MQ=255)
                                                                        gtaAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAAc                                                                                                                                   >  4:1677405/1‑93 (MQ=255)
                                                                        gtaAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAAc                                                                                                                                   <  3:1677405/93‑1 (MQ=255)
                                                                              cTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACt                                                                                                             <  4:151690/109‑1 (MQ=255)
                                                                              cTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACt                                                                                                             >  3:151690/1‑109 (MQ=255)
                                                                               ttAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGAc                                                                                                              >  3:790859/1‑107 (MQ=255)
                                                                                tAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAAtttt                                                                   <  2:1099926/149‑1 (MQ=255)
                                                                                 aaTTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGaaaa                                                                                               <  1:1084015/120‑1 (MQ=255)
                                                                                 aaTTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGaaaa                                                                                               >  2:1084015/1‑120 (MQ=255)
                                                                                 aaTTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTa                                                                  <  3:1031944/149‑1 (MQ=255)
                                                                                     aCCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGaaaa                                                                                               >  3:94049/1‑116 (MQ=255)
                                                                                     aCCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGaaaa                                                                                               <  4:94049/116‑1 (MQ=255)
                                                                                     aCCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTgg                                                              <  2:909570/149‑1 (MQ=255)
                                                                                        atcatGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCaa                                                                                                     <  1:1612608/107‑1 (MQ=255)
                                                                                        atcatGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCaa                                                                                                     >  2:1612608/1‑107 (MQ=255)
                                                                                         tcatGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACg                                                                                                                                  >  1:101778/1‑77 (MQ=255)
                                                                                         tcatGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGaaaa                                                                                               >  4:929969/1‑112 (MQ=255)
                                                                                         tcatGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGc                                                          >  2:903105/1‑149 (MQ=255)
                                                                                           atGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGt                                                        <  3:159814/149‑1 (MQ=255)
                                                                                            tGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGtt                                                       >  1:1244453/1‑149 (MQ=255)
                                                                                             gAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGTAGAACGCACGTCTGc                                                                                                                         <  4:1833503/82‑1 (MQ=255)
                                                                                             gAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGTAGAACGCACGTCTGc                                                                                                                         >  3:1833503/1‑82 (MQ=255)
                                                                                             gAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGc                                                                                                                         >  1:179988/1‑82 (MQ=255)
                                                                                             gAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGc                                                                                                                         <  2:179988/82‑1 (MQ=255)
                                                                                                   aaTATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGa                                                                                       <  4:1082255/110‑1 (MQ=255)
                                                                                                   aaTATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGtt                                                       <  1:1555289/142‑1 (MQ=255)
                                                                                                   aaTATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGtt                                                       >  2:1555289/1‑142 (MQ=255)
                                                                                                   aaTATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCa                                                >  4:807058/1‑149 (MQ=255)
                                                                                                   aaTATTTTTACACTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGa                                                                                       >  3:1082255/1‑110 (MQ=255)
                                                                                                        ttttACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCt                                                                                                                                             <  1:690029/51‑1 (MQ=255)
                                                                                                         tttACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCtttt                                          >  1:588546/1‑149 (MQ=255)
                                                                                                          ttACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTa                                         <  1:903105/149‑1 (MQ=255)
                                                                                                                    gcgcAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTgg                               >  3:983172/1‑149 (MQ=255)
                                                                                                                    gcgcAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTgg                               >  1:298669/1‑149 (MQ=255)
                                                                                                                    gcgcAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGATGGTTgg                               <  2:1032997/149‑1 (MQ=255)
                                                                                                                     cgcAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGc                              <  3:1308945/149‑1 (MQ=255)
                                                                                                                      gcAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGaaa                                                                                                >  1:1787243/1‑82 (MQ=255)
                                                                                                                       cAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTg                            >  3:1666295/1‑149 (MQ=255)
                                                                                                                                cTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAaga                                                                                             >  2:70320/1‑75 (MQ=255)
                                                                                                                                cTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAaga                                                                                             <  1:70320/75‑1 (MQ=255)
                                                                                                                                cTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCg                   >  2:248487/1‑149 (MQ=255)
                                                                                                                                  ttGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCa                 >  2:122538/1‑149 (MQ=255)
                                                                                                                                  ttGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCa                 >  3:1729367/1‑149 (MQ=255)
                                                                                                                                       gATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTa                                         >  4:370513/1‑120 (MQ=255)
                                                                                                                                       gATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTGaa            >  4:1463831/1‑149 (MQ=255)
                                                                                                                                       gATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTGaa            <  2:1244453/149‑1 (MQ=255)
                                                                                                                                          aaTCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTGAAGTc         <  3:807058/149‑1 (MQ=255)
                                                                                                                                           tgctcCAGCCATGCCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGaaaa                                                                       >  4:1222506/3‑86 (MQ=255)
                                                                                                                                           tgctcCAGCCATGCCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGaaaa                                                                       <  3:1222506/84‑1 (MQ=255)
                                                                                                                                             ctcCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTGAAGTCGCt      >  3:1141966/1‑149 (MQ=255)
                                                                                                                                              tcCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAgtgt                                                                                         >  2:1475919/1‑65 (MQ=255)
                                                                                                                                              tcCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAgtgt                                                                                         <  1:1475919/65‑1 (MQ=255)
                                                                                                                                                cAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTGAAGTCGCTGCa   <  2:588546/149‑1 (MQ=255)
                                                                                                                                                 aGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTGAAGTCGCTGCAc  >  1:1586496/1‑149 (MQ=255)
                                                                                                                                                 aGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTGAAGTCGCTGCAc  >  3:883044/1‑149 (MQ=255)
                                                                                                                                                 aGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTGAAGTCGCTGCAc  <  4:819011/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GATAATCTCGGTTTATTAACCATGCAGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTGAAGTCGCTGCAC  >  NC_000913/1445469‑1445762

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: