Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,247,144 T→A 5.2% D78V (GAT→GTT)  treA ← periplasmic trehalase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,247,1440TA 5.2% 294.8 / 2.4 88D78V (GAT→GTT) treAperiplasmic trehalase
Reads supporting (aligned to +/- strand):  ref base T (41/42);  new base A (3/2);  total (44/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCAGGCGGCTGTGGTGTTACCGGTGTTTCTTC  >  NC_000913/1246997‑1247286
                                                                                                                                                   |                                                                                                                                              
gTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGCGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAAt                                                                                                                                                <  4:166628/148‑1 (MQ=255)
  tAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAAt                                                                                                                                                >  1:137090/1‑146 (MQ=255)
      aCCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGc                                                                                                                                            >  1:169590/1‑146 (MQ=255)
           ccAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGatc                                                                                                                                     <  4:319781/148‑1 (MQ=255)
             aaaGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGatcat                                                                                                                                   >  4:103126/1‑148 (MQ=255)
               aGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAAACAGCAAGGatcat                                                                                                                                   >  4:386803/1‑146 (MQ=255)
                gTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAg                                                                                                                                <  8:6277/148‑1 (MQ=255)
                 tCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGc                                                                                                                               <  6:169413/148‑1 (MQ=255)
                      cAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCGAGGATCATCAGCGGATc                                                                                                                          >  7:403347/1‑148 (MQ=255)
                                gcgcAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTcggc                                                                                                                <  5:276318/148‑1 (MQ=255)
                                 cgcAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTcggca                                                                                                               <  5:271405/148‑1 (MQ=255)
                                     gTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTcggcacggc                                                                                                           <  2:151230/148‑1 (MQ=255)
                                      tGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTcggcacggca                                                                                                          >  8:274659/1‑148 (MQ=255)
                                          tGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGatc                                                                                                                                     >  3:392639/1‑117 (MQ=255)
                                          tGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGatc                                                                                                                                     <  4:392639/117‑1 (MQ=255)
                                           gCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCgg                                                                                                      >  8:284093/1‑147 (MQ=255)
                                              ccTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCaaa                                                                                                  <  7:415087/148‑1 (MQ=255)
                                                 ctGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGt                                                                                               >  1:422407/1‑148 (MQ=255)
                                                 ctGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAAACAGCAAGGATAATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGt                                                                                               >  6:49155/1‑148 (MQ=255)
                                                    gCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGt                                                                                               <  1:113262/145‑1 (MQ=255)
                                                     cGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATc                                                                                                        >  2:420935/1‑135 (MQ=255)
                                                     cGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATc                                                                                                        <  1:420935/135‑1 (MQ=255)
                                                      ggAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCGAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGtttttt                                                                                          <  8:403347/148‑1 (MQ=255)
                                                          cATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGCTCGGCACGGCATCGGCAAAGGTTTTTTGGTc                                                                                      >  7:415501/1‑148 (MQ=255)
                                                             atTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGCCAAAGGTTATTTGGGCCg                                                                                    >  8:329685/1‑147 (MQ=255)
                                                             atTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGATAACCAAATGGCGCAGATCAAATCCGCTCTGGTTATGCTGAATCCGATAATCAGAAAGGATCATCAGCGTATCGCTGATCGGCACGGCATCGGCAAAGGTTTTTTGGTCCg                                                                                    >  2:354/1‑147 (MQ=255)
                                                              tttCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCgg                                                                                   <  7:274659/147‑1 (MQ=255)
                                                               ttCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGa                                                                                  >  4:300318/1‑147 (MQ=255)
                                                                   cGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAg                                                                             >  5:949/1‑148 (MQ=255)
                                                                    gCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATTGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGaaaa                                                                               <  2:169590/145‑1 (MQ=255)
                                                                    gCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGaaaa                                                                               <  2:422407/145‑1 (MQ=255)
                                                                    gCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAg                                                                             <  6:949/147‑1 (MQ=255)
                                                                    gCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAAACAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGaaaa                                                                               <  5:49155/145‑1 (MQ=255)
                                                                      cttcttTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCg                                                                                    <  3:48340/138‑1 (MQ=255)
                                                                      cttcttTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCg                                                                                    >  4:48340/1‑138 (MQ=255)
                                                                       ttcttTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTg                                                                         >  3:72731/1‑148 (MQ=255)
                                                                        tcttTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGttttt                                                                                           <  4:290355/129‑1 (MQ=255)
                                                                        tcttTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGttttt                                                                                           >  3:290355/1‑129 (MQ=255)
                                                                        acttACGGCAGGGTGAAATAGACGTAAACTAAATGGCCCAGATCAAATCCGCTCTCGTTTTGCTGCATCCGATAAACAGCAAGCATCATCAGCGGATCGCTCTACGGCACGGAATCGGCAAAGGTTTTTTGCTCCCGAAAAAGATTgg                                                                        >  2:362301/2‑148 (MQ=255)
                                                                         cttTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTg                                                                         <  2:137090/146‑1 (MQ=255)
                                                                            tCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAAACAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGt                                                                     <  3:386803/147‑1 (MQ=255)
                                                                                 aGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGca                                                               <  4:72731/148‑1 (MQ=255)
                                                                                    gTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTg                                                                 >  6:402831/1‑143 (MQ=255)
                                                                                    gTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTg                                                                 <  5:402831/143‑1 (MQ=255)
                                                                                            gACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTATTGCACATCATTaaa                                                     >  3:306615/1‑147 (MQ=255)
                                                                                                ttAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCg                                                >  7:440362/1‑148 (MQ=255)
                                                                                                 tAACGAAATGGCGCAGATCAAATCCGCACTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGTCATCGGCAAAGGTTTTTAGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCgg                                               >  5:11634/1‑148 (MQ=255)
                                                                                                  aaCGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCgg                                               <  4:409733/147‑1 (MQ=255)
                                                                                                   aCGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGc                                              <  7:284093/147‑1 (MQ=255)
                                                                                                     gAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTaaa                                                     <  7:258602/138‑1 (MQ=255)
                                                                                                     gAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTaaa                                                     >  8:258602/1‑138 (MQ=255)
                                                                                                     gAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGccc                                            >  8:92952/1‑147 (MQ=255)
                                                                                                      aaaTGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATGAGCGGCTCGCTGTTCGGCAAGGCATCGGCAAAGGTTTTTTGGTCCGGAAATAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCta                                          <  1:354/148‑1 (MQ=255)
                                                                                                      aaaTGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTc                                                                                                                   <  7:376166/75‑1 (MQ=255)
                                                                                                      aaaTGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTc                                                                                                                   >  8:376166/1‑75 (MQ=255)
                                                                                                       aaTGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGt                                                                                               >  7:88987/1‑94 (MQ=255)
                                                                                                       aaTGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGt                                                                                               <  8:88987/94‑1 (MQ=255)
                                                                                                            cgcAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATaaaa                                    >  8:398825/1‑148 (MQ=255)
                                                                                                              cAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAata                                  >  6:82882/1‑148 (MQ=255)
                                                                                                                gATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATc                                <  5:82882/148‑1 (MQ=255)
                                                                                                                 aTCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGc                                                                >  5:207569/1‑115 (MQ=255)
                                                                                                                 aTCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGc                                                                <  6:207569/115‑1 (MQ=255)
                                                                                                                 aTCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAatat                                 >  5:261443/1‑146 (MQ=255)
                                                                                                                 aTCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAatat                                 >  8:242102/1‑146 (MQ=255)
                                                                                                                  tCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTaa                                                      <  2:64513/124‑1 (MQ=255)
                                                                                                                  tCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTaa                                                      >  1:64513/1‑124 (MQ=255)
                                                                                                                  tCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCAg                              >  8:58930/1‑148 (MQ=255)
                                                                                                                  tCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCAg                              >  4:408765/1‑148 (MQ=255)
                                                                                                                  tCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCAg                              >  8:13158/1‑148 (MQ=255)
                                                                                                                  tCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGACATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCAg                              >  5:163189/1‑148 (MQ=255)
                                                                                                                  tCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGAACGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGaaaa                                                                               >  1:30820/1‑99 (MQ=255)
                                                                                                                  tCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGAACGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGaaaa                                                                               <  2:30820/99‑1 (MQ=255)
                                                                                                                      aTCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTaaa                                                     >  3:276981/1‑121 (MQ=255)
                                                                                                                      aTCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTaaa                                                     <  4:276981/121‑1 (MQ=255)
                                                                                                                       tCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCAggcggc                         <  7:13158/148‑1 (MQ=255)
                                                                                                                           ctctGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCCCGGCATCGGCAAAGGTTTTTTGGTCCGGCATAAGGTTGGCGTTTTGCACAGCATTACCCAGCGGCCCTAATAAAATATCAGGCGGCtgtg                     <  1:362301/148‑1 (MQ=255)
                                                                                                                           ctctGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCAGGCGGCtgtg                     <  3:300318/148‑1 (MQ=255)
                                                                                                                              tGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGCTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCAGGCGGCTgtgg                    <  8:415501/146‑1 (MQ=255)
                                                                                                                               ggTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCAGGCGGCTGTGgtgt                 >  7:161649/1‑148 (MQ=255)
                                                                                                                                gTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCAGGCGGCTGTGgtgt                 <  8:161649/147‑1 (MQ=255)
                                                                                                                                   ttGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCAGGCGGCTGTGGTGTTAc              <  7:242102/147‑1 (MQ=255)
                                                                                                                                    tgctgcATCCGATAATCAGCAAGGATCTTCAGCGGATCGCTGTTCGGCACGCCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCGCTAATAAAATATCAGGCGGCTGTGGTGTTACCg            >  6:331063/1‑148 (MQ=255)
                                                                                                                                     gctgcATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTcaa                                                     <  2:307990/106‑1 (MQ=255)
                                                                                                                                     gctgcATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTcaa                                                     >  1:307990/1‑106 (MQ=255)
                                                                                                                                         cATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCAGGCGGCTGTGGTGTTACCGgtg         <  8:440362/146‑1 (MQ=255)
                                                                                                                                         cATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCAGGCGGCTGTGGTGTTACCGGTGtt       <  7:126624/148‑1 (MQ=255)
                                                                                                                                            ccGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAat                                   <  3:399884/117‑1 (MQ=255)
                                                                                                                                            ccGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAat                                   >  4:399884/1‑117 (MQ=255)
                                                                                                                                             cGATAATCAGCAAGGATCATCAGCGGATCGGTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCAGGCGGCTGTGGTGTTACCGGTGTttctt   <  8:97636/148‑1 (MQ=255)
                                                                                                                                              gATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCAGGCGGCTGTGGTGTTACCGGTGTttcttc  >  3:112238/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
GTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCAGGCGGCTGTGGTGTTACCGGTGTTTCTTC  >  NC_000913/1246997‑1247286

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