Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,445,621 G→C 6.5% 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 6.5% 446.2 / 17.5 123W857C (TGG→TGCydbHputative membrane‑anchored protein, function unknown
Reads supporting (aligned to +/- strand):  ref base G (61/54);  new base C (4/4);  total (65/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.97e-01

Only 100 of 123 total aligned reads displayed.

GTTTATTAACCATGCAGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTGAAGTCGCTGCAC  >  NC_000913/1445479‑1445762
                                                                                                                                              |                                                                                                                                             
     ttaACCATGCAGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTTGAGCAGaa                                                                                                                                    >  1:909570/1‑149 (MQ=255)
         ccATGCAGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGTATGGCTGGAGCAGacgtgc                                                                                                                                <  3:1569169/149‑6 (MQ=255)
         ccATGCAGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTg                                                                                                                                            >  2:565673/1‑137 (MQ=255)
         ccATGCAGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTg                                                                                                                                            <  1:565673/137‑1 (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)
                 gCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGcc                                                                                                                        >  2:809022/1‑149 (MQ=255)
                 gCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGcc                                                                                                                        >  3:1696564/1‑149 (MQ=255)
                  ccAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCg                                                                                                                       >  1:202014/1‑149 (MQ=255)
                   cAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCgg                                                                                                                      <  3:800386/149‑1 (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)
                                   aGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGaaa                                                                                                                  <  3:1147902/137‑1 (MQ=255)
                                   aGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGaaa                                                                                                                  >  4:1147902/1‑137 (MQ=255)
                                     tCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAg                                                                                                    >  4:159814/1‑149 (MQ=255)
                                     tCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAg                                                                                                    >  1:75080/1‑149 (MQ=255)
                                     tCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGCCATGCCTGGAgatta                                                                                                                                     >  3:1828475/1‑112 (MQ=255)
                                     tCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGCCATGCCTGGAgatta                                                                                                                                     <  4:1828475/116‑5 (MQ=255)
                                            gATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAaga                                                                                             <  1:76951/149‑1 (MQ=255)
                                            gATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAaga                                                                                             <  4:56983/149‑1 (MQ=255)
                                             aTGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTcc                                                                                                         <  4:1410712/136‑1 (MQ=255)
                                             aTGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTcc                                                                                                         >  3:1410712/1‑136 (MQ=255)
                                             aTGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAagag                                                                                            >  2:1439738/1‑149 (MQ=255)
                                              tGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTcc                                                                                                         >  3:321697/1‑135 (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                                                                                                    <  1:1584736/137‑1 (MQ=255)
                                                 tAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAg                                                                                                    >  2:1584736/1‑137 (MQ=255)
                                                  aaaaGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGc                                                                                                       <  4:825483/133‑1 (MQ=255)
                                                  aaaaGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGc                                                                                                       >  3:825483/1‑133 (MQ=255)
                                                   aaaGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAg                                                                                      <  2:75080/149‑1 (MQ=255)
                                                        ggtacggtaAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGaaaa                                                                                 >  3:819011/1‑149 (MQ=255)
                                                         gtacggtaAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGaaaaa                                                                                >  4:1031944/1‑149 (MQ=255)
                                                         gtacggtaAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGaaaaa                                                                                <  3:1046274/149‑1 (MQ=255)
                                                          tacggtaAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAAc                                                                               >  1:1099926/1‑149 (MQ=255)
                                                             ggtaAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGaaaaa                                                                                <  2:1725322/145‑1 (MQ=255)
                                                             ggtaAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGaaaaa                                                                                >  1:1725322/1‑145 (MQ=255)
                                                             ggtaAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAAt                                                                            <  4:1561228/149‑1 (MQ=255)
                                                              gtaAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAAc                                                                                                                                   >  4:1677405/1‑93 (MQ=255)
                                                              gtaAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAAc                                                                                                                                   <  3:1677405/93‑1 (MQ=255)
                                                                aacccTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATTaa                                                                         <  3:1255412/146‑1 (MQ=255)
                                                                    cTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACt                                                                                                             <  4:151690/109‑1 (MQ=255)
                                                                    cTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACt                                                                                                             >  3:151690/1‑109 (MQ=255)
                                                                     ttAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGAc                                                                                                              <  4:790859/107‑1 (MQ=255)
                                                                      tAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAAtttt                                                                   <  2:1099926/149‑1 (MQ=255)
                                                                       aaTTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGaaaa                                                                                               <  1:1084015/120‑1 (MQ=255)
                                                                       aaTTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGaaaa                                                                                               >  2:1084015/1‑120 (MQ=255)
                                                                          tACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTg                                                               <  4:83318/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                                                                                               <  3:929969/112‑1 (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                                                                                                                         >  3:1833503/1‑82 (MQ=255)
                                                                                   gAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGTAGAACGCACGTCTGc                                                                                                                         <  4:1833503/82‑1 (MQ=255)
                                                                                    aaGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGa                                                     >  1:1543141/1‑149 (MQ=255)
                                                                                     agagAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGAc                                                    <  1:684370/149‑1 (MQ=255)
                                                                                         aaTATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGa                                                                                       <  4:1082255/110‑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)
                                                                                              ttttACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCt                                                                                                                                             >  2:690029/1‑51 (MQ=255)
                                                                                                ttACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTa                                         <  1:903105/149‑1 (MQ=255)
                                                                                                          gcgcAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTgg                               >  1:298669/1‑149 (MQ=255)
                                                                                                          gcgcAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTgg                               >  3:983172/1‑149 (MQ=255)
                                                                                                          gcgcAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGATGGTTgg                               <  2:1032997/149‑1 (MQ=255)
                                                                                                           cgcAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGc                              <  3:1308945/149‑1 (MQ=255)
                                                                                                            gcAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGaaa                                                                                                <  2:1787243/82‑1 (MQ=255)
                                                                                                            gcAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGaaa                                                                                                >  1:1787243/1‑82 (MQ=255)
                                                                                                             cAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTg                            >  3:1666295/1‑149 (MQ=255)
                                                                                                                      cTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAaga                                                                                             <  1:70320/75‑1 (MQ=255)
                                                                                                                      cTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAaga                                                                                             >  2:70320/1‑75 (MQ=255)
                                                                                                                      cTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCg                   >  2:248487/1‑149 (MQ=255)
                                                                                                                        ttGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCa                 >  2:122538/1‑149 (MQ=255)
                                                                                                                         tGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCAt                <  2:298669/149‑1 (MQ=255)
                                                                                                                             gATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTa                                         <  3:370513/120‑1 (MQ=255)
                                                                                                                             gATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTa                                         >  4:370513/1‑120 (MQ=255)
                                                                                                                             gATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTGaa            <  2:1244453/149‑1 (MQ=255)
                                                                                                                             gATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTGaa            >  4:1463831/1‑149 (MQ=255)
                                                                                                                             gATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTGaa            >  1:22818/1‑149 (MQ=255)
                                                                                                                                aaTCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTGAAGTc         <  3:807058/149‑1 (MQ=255)
                                                                                                                                 tgctcCAGCCATGCCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGaaaa                                                                       >  4:1222506/3‑86 (MQ=255)
                                                                                                                                   ctcCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTGAAGTCGCt      >  3:1141966/1‑149 (MQ=255)
                                                                                                                                    tcCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAgtgt                                                                                         <  1:1475919/65‑1 (MQ=255)
                                                                                                                                    tcCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAgtgt                                                                                         >  2:1475919/1‑65 (MQ=255)
                                                                                                                                       aGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTGAAGTCGCTGCAc  >  3:883044/1‑149 (MQ=255)
                                                                                                                                       aGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTGAAGTCGCTGCAc  >  1:1586496/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                             
GTTTATTAACCATGCAGGCCAACATTACAGGTACCAGTCGCGTTGATGGTAAAAGCGGTACGGTAAACCTTAATTACCATCATGAAGAGAATATTTTTACGCTGTGGCGCAGTTTACGCTTTGGCGATAATCTCCAGGCATGGCTGGAGCAGAACGCACGTCTGCCGGGAAATGACTGTCCGCAAGGAAAAGAGTGTGAGGAAAAACAATGAAAATTTTACTGGCTGCGTTGACGTCATCTTTTATGCTGGTTGGCTGTACGCCTCGCATTGAAGTCGCTGCAC  >  NC_000913/1445479‑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: