Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 423,995 C→G 5.4% P493R (CCG→CGG) ‡ malZ → maltodextrin glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913423,9950CG 5.4% 663.2 / 19.5 182P493R (CCG→CGG) ‡malZmaltodextrin glucosidase
Reads supporting (aligned to +/- strand):  ref base C (83/89);  new base G (5/5);  total (88/94)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

Only 100 of 183 total aligned reads displayed.

TCGACAGCCACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTG  >  NC_000913/423848‑424143
                                                                                                                                                   |                                                                                                                                                    
   aCAGCCACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGttt                                                                                                                                                  <  1:966405/149‑1 (MQ=255)
     aGCCACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGGCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTg                                                                                                                                                >  1:942208/1‑149 (MQ=255)
      gCCACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAACGGATCGTTTTTGc                                                                                                                                               <  2:762612/149‑1 (MQ=255)
                tGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGt                                                                                                                                     >  1:1287200/1‑149 (MQ=255)
                  cgcgATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTc                                                                                                                                   <  2:260565/149‑1 (MQ=255)
                     gATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTcccc                                                                                                                                <  4:1186201/149‑1 (MQ=255)
                     gATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTcccc                                                                                                                                <  4:551232/149‑1 (MQ=255)
                      aTTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCt                                                                                                                               <  4:1451356/149‑1 (MQ=255)
                      aTTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCt                                                                                                                               <  2:1287200/149‑1 (MQ=255)
                        ttAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTgg                                                                                                                             <  4:1178138/149‑1 (MQ=255)
                          aaaaCGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGt                                                                                                                                     <  1:391932/139‑1 (MQ=255)
                          aaaaCGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGt                                                                                                                                     >  2:391932/1‑139 (MQ=255)
                            aaCGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAgg                                                                                                                         >  3:611798/1‑149 (MQ=255)
                               gctgctCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAACGGATCGTTTTTGc                                                                                                                                               >  1:987212/1‑124 (MQ=255)
                               gctgctCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAACGGATCGTTTTTGc                                                                                                                                               <  2:987212/124‑1 (MQ=255)
                                      ggtcggGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCa                                                                                                               <  2:942208/149‑1 (MQ=255)
                                       gtcggGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAg                                                                                                              >  2:1042591/1‑149 (MQ=255)
                                       gtcggGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAg                                                                                                              <  3:696696/149‑1 (MQ=255)
                                       gtcggGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAg                                                                                                              <  4:581963/149‑1 (MQ=255)
                                       gtcggGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAg                                                                                                              <  4:773326/149‑1 (MQ=255)
                                       gtcggGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAg                                                                                                              <  3:649229/149‑1 (MQ=255)
                                                  gcgcgcCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGtt                                                                                                   >  2:921321/1‑149 (MQ=255)
                                                      gcCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTaaa                                                                                                                                         <  4:963567/107‑1 (MQ=255)
                                                          gCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGACGTTATTCgcgc                                                                                           <  4:918908/149‑1 (MQ=255)
                                                              cTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGt                                                                                                                        >  2:324734/1‑116 (MQ=255)
                                                              cTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGt                                                                                                                        <  1:324734/116‑1 (MQ=255)
                                                               tGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGtt                                                                                                                                                   >  2:552784/1‑88 (MQ=255)
                                                               tGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTAc                                                                                      >  1:683249/1‑149 (MQ=255)
                                                                ggCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTAcc                                                                                     >  1:100445/1‑149 (MQ=255)
                                                                   ggtggtCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGc                                                                                  <  4:780826/149‑1 (MQ=255)
                                                                       gtCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAg                                                                                   <  3:910461/144‑1 (MQ=255)
                                                                       gtCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAg                                                                                   >  2:304402/1‑144 (MQ=255)
                                                                       gtCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAg                                                                                   <  1:304402/144‑1 (MQ=255)
                                                                       gtCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAAt                                                                              <  1:921321/149‑1 (MQ=255)
                                                                          tGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCa                                                                                                               >  2:1332058/1‑113 (MQ=255)
                                                                          tGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGAt                                                                           >  4:553367/1‑149 (MQ=255)
                                                                           ggCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGAtt                                                                          <  4:1335393/149‑1 (MQ=255)
                                                                            gCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGAt                                                                           <  1:1314551/147‑1 (MQ=255)
                                                                            gCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTg                                                                         >  2:688101/1‑149 (MQ=255)
                                                                              tGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGaa                                                                               >  4:411481/1‑141 (MQ=255)
                                                                              tGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTgc                                                                        <  3:546467/148‑1 (MQ=255)
                                                                                ttCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCt                                                                     >  3:780206/1‑149 (MQ=255)
                                                                                  cACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGc                                                                   <  2:683249/149‑1 (MQ=255)
                                                                                   aCCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCt                                                                                                                               <  1:1063815/88‑1 (MQ=255)
                                                                                     ctggcctggTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGtttt                                                                                                                                                 >  2:596676/1‑68 (MQ=255)
                                                                                      tggcctggTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTaa                                                               <  4:686494/149‑1 (MQ=255)
                                                                                          ctggTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGaaa                                                           >  4:1461141/1‑149 (MQ=255)
                                                                                           tggTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGt                                                                                                                        >  3:131141/1‑87 (MQ=255)
                                                                                           tggTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCa                                                                                                               >  1:1334357/1‑96 (MQ=255)
                                                                                           tggTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGaa                                                                               <  1:431985/128‑1 (MQ=255)
                                                                                           tggTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGaa                                                                               >  2:431985/1‑128 (MQ=255)
                                                                                           tggTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGaaaa                                                          <  1:1174128/149‑1 (MQ=255)
                                                                                           tggTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGaaaa                                                          >  4:782120/1‑149 (MQ=255)
                                                                                            ggTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAg                                                         <  3:1461141/149‑1 (MQ=255)
                                                                                              tgtACCGTGCATTTATTATGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGt                                                                                                                        <  3:902468/84‑1 (MQ=255)
                                                                                              tgtACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCATAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGt                                                                                                                        >  4:902468/1‑84 (MQ=255)
                                                                                              tgtACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTc                                                       <  3:93883/149‑1 (MQ=255)
                                                                                                tACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGaa                                                                               <  2:840015/123‑1 (MQ=255)
                                                                                                tACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGaa                                                                               >  1:840015/1‑123 (MQ=255)
                                                                                                 aCCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCa                                                                                                               >  2:343015/1‑90 (MQ=255)
                                                                                                 aCCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCa                                                                                                               <  1:343015/90‑1 (MQ=255)
                                                                                                    gTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTaa                                                               >  4:351976/1‑135 (MQ=255)
                                                                                                     tGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGaa                                                                               <  4:530209/118‑1 (MQ=255)
                                                                                                     tGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCt                                                <  3:573143/149‑1 (MQ=255)
                                                                                                     tGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCt                                                <  1:688101/149‑1 (MQ=255)
                                                                                                      gCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTa                                               >  4:309566/1‑149 (MQ=255)
                                                                                                         tttATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGc                                                                                                                            <  3:499412/69‑1 (MQ=255)
                                                                                                            attaCGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGaa                                                                               <  1:255443/111‑1 (MQ=255)
                                                                                                                     atgaAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAg                                                                                                                 <  3:152208/68‑1 (MQ=255)
                                                                                                                          gTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGaaaa                                                          >  2:135420/1‑118 (MQ=255)
                                                                                                                           tAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCt                                  >  4:1468353/1‑141 (MQ=255)
                                                                                                                           tAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTc                               >  3:1426898/1‑144 (MQ=255)
                                                                                                                            aGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGc                         >  1:1036955/1‑149 (MQ=255)
                                                                                                                              gACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGaaaa                                                          <  3:472832/114‑1 (MQ=255)
                                                                                                                              gACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAg                             <  2:1107744/143‑1 (MQ=255)
                                                                                                                              gACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTg                       >  4:497984/1‑149 (MQ=255)
                                                                                                                                   gATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCa                                                                                                               <  1:1124246/56‑1 (MQ=255)
                                                                                                                                   gATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGt                                                                 <  3:1080032/102‑1 (MQ=255)
                                                                                                                                   gATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGt                                                                 >  4:1080032/1‑102 (MQ=255)
                                                                                                                                   gATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGc                  <  1:579829/149‑1 (MQ=255)
                                                                                                                                   gATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGc                  <  3:371498/149‑1 (MQ=255)
                                                                                                                                   gATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGc                  >  4:288660/1‑149 (MQ=255)
                                                                                                                                    aTGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCg                 >  3:975785/1‑149 (MQ=255)
                                                                                                                                      ggCAAAACGGATCGTTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTc                                           >  3:1240924/1‑121 (MQ=255)
                                                                                                                                      ggCAAAACGGATCGTTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGa               >  3:1272351/1‑149 (MQ=255)
                                                                                                                                       gCAAAACGGATCGTTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCt                                  <  4:824361/129‑1 (MQ=255)
                                                                                                                                       gCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGaa              <  3:782120/149‑1 (MQ=255)
                                                                                                                                       gCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGaa              <  1:60207/149‑1 (MQ=255)
                                                                                                                                        cAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCa                                                      <  3:985180/108‑1 (MQ=255)
                                                                                                                                         aaaaaCGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCa                                          >  4:683600/1‑119 (MQ=255)
                                                                                                                                         aaaaaCGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCa                                          <  3:683600/119‑1 (MQ=255)
                                                                                                                                         aaaaaCGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTc                               <  2:712717/130‑1 (MQ=255)
                                                                                                                                         aaaaaCGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTc                               >  1:712717/1‑130 (MQ=255)
                                                                                                                                             aCGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAAc        >  1:1376087/1‑149 (MQ=255)
                                                                                                                                               gATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACgt      >  3:377040/1‑149 (MQ=255)
                                                                                                                                                aTCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCt                                  <  3:347372/120‑1 (MQ=255)
                                                                                                                                                aTCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGaa              >  2:1212595/1‑140 (MQ=255)
                                                                                                                                                 tCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGc                                                                   >  2:1130300/1‑86 (MQ=255)
                                                                                                                                                 tCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGt                      <  1:747778/131‑1 (MQ=255)
                                                                                                                                                   cGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACgtggtg  >  2:88512/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TCGACAGCCACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTG  >  NC_000913/423848‑424143

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