Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 423,990 C→G 5.6% N491K (AAC→AAG) ‡ malZ → maltodextrin glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913423,9900CG 5.6% 650.3 / 22.0 179N491K (AAC→AAG) ‡malZmaltodextrin glucosidase
Reads supporting (aligned to +/- strand):  ref base C (81/88);  new base G (5/5);  total (86/93)
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.84e-01

Only 100 of 180 total aligned reads displayed.

ATCAGCTCGACAGCCACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAAC  >  NC_000913/423842‑424137
                                                                                                                                                    |                                                                                                                                                   
   aGCTCGACAGCCACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGAt                                                                                                                                                  >  2:966405/1‑149 (MQ=255)
   aGCTCGACAGCCACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGAt                                                                                                                                                  <  1:926854/149‑1 (MQ=255)
    gCTCGACAGCCACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATc                                                                                                                                                 >  2:718454/1‑149 (MQ=255)
      tCGACAGCCACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCg                                                                                                                                               <  1:660265/149‑1 (MQ=255)
         aCAGCCACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATGTTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGttt                                                                                                                                            <  3:1252391/149‑1 (MQ=255)
         aCAGCCACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGttt                                                                                                                                            >  3:918908/1‑149 (MQ=255)
         aCAGCCACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGttt                                                                                                                                            <  1:966405/149‑1 (MQ=255)
                      tGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGt                                                                                                                               >  1:1287200/1‑149 (MQ=255)
                           gATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTcccc                                                                                                                          <  4:551232/149‑1 (MQ=255)
                           gATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTcccc                                                                                                                          <  4:1186201/149‑1 (MQ=255)
                           gATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTcccc                                                                                                                          <  2:804722/149‑1 (MQ=255)
                            aTTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCt                                                                                                                         <  2:1287200/149‑1 (MQ=255)
                              ttAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTgg                                                                                                                       <  4:1178138/149‑1 (MQ=255)
                              ttAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTgg                                                                                                                       >  3:823244/1‑149 (MQ=255)
                                aaaaCGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGt                                                                                                                               <  1:391932/139‑1 (MQ=255)
                                  aaCGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATc                                                                                                                                                 <  2:621042/119‑1 (MQ=255)
                                  aaCGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAgg                                                                                                                   >  3:611798/1‑149 (MQ=255)
                                  aaCGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAgg                                                                                                                   >  3:780826/1‑149 (MQ=255)
                                     gctgctCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAACGGATCGTTTTTGc                                                                                                                                         >  1:987212/1‑124 (MQ=255)
                                     gctgctCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAACGGATCGTTTTTGc                                                                                                                                         <  2:987212/124‑1 (MQ=255)
                                       tgctCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGa                                                                                                                                                   >  3:996172/1‑112 (MQ=255)
                                             gtcggGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAg                                                                                                        <  3:649229/149‑1 (MQ=255)
                                             gtcggGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAg                                                                                                        <  3:696696/149‑1 (MQ=255)
                                                  gATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAACGGATCGTTTTTGcc                                                                                                                                        <  2:840964/112‑1 (MQ=255)
                                                        gcgcgcCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGtt                                                                                             >  2:921321/1‑149 (MQ=255)
                                                          gcgcCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGaaa                                                                                                             >  2:298661/1‑131 (MQ=255)
                                                            gcCTGCCGCTGACGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCCAAAACGATCCGTTTTGCCGTaaa                                                                                                                                   >  3:963567/1‑107 (MQ=255)
                                                                    cTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGt                                                                                                                  <  1:324734/116‑1 (MQ=255)
                                                                    cTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGt                                                                                                                  >  2:324734/1‑116 (MQ=255)
                                                                     tGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGtt                                                                                                                                             <  1:552784/88‑1 (MQ=255)
                                                                     tGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGtt                                                                                                                                             >  2:552784/1‑88 (MQ=255)
                                                                      ggCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTAcc                                                                               >  1:100445/1‑149 (MQ=255)
                                                                       gCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCa                                                                              >  4:93883/1‑149 (MQ=255)
                                                                            ggtCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGaa                                                                         <  4:611798/149‑1 (MQ=255)
                                                                            ggtCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGaa                                                                         <  2:100445/149‑1 (MQ=255)
                                                                             gtCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAg                                                                             <  3:910461/144‑1 (MQ=255)
                                                                             gtCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAg                                                                             <  1:304402/144‑1 (MQ=255)
                                                                             gtCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAg                                                                             >  4:910461/1‑144 (MQ=255)
                                                                             gtCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAAt                                                                        <  1:921321/149‑1 (MQ=255)
                                                                                tGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCa                                                                                                         <  1:1332058/113‑1 (MQ=255)
                                                                                tGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCa                                                                                                         >  2:1332058/1‑113 (MQ=255)
                                                                                    tGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGaa                                                                         <  3:411481/141‑1 (MQ=255)
                                                                                     gTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTgcgc                                                                <  2:314369/149‑1 (MQ=255)
                                                                                        cACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGc                                                             <  2:683249/149‑1 (MQ=255)
                                                                                         aCCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCt                                                                                                                         <  1:1063815/88‑1 (MQ=255)
                                                                                         aCCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCt                                                                                                                         >  2:1063815/1‑88 (MQ=255)
                                                                                           ctggcctggTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGtttt                                                                                                                                           <  1:596676/68‑1 (MQ=255)
                                                                                           ctggcctggTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGtttt                                                                                                                                           >  2:596676/1‑68 (MQ=255)
                                                                                            tggcctggTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTaa                                                         <  4:686494/149‑1 (MQ=255)
                                                                                                ctggTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGaaa                                                     >  4:1461141/1‑149 (MQ=255)
                                                                                                 tggTGTACCGTGCATTTATTACGGTGATGATGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGt                                                                                                                  <  4:131141/87‑1 (MQ=255)
                                                                                                 tggTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGt                                                                                                                  >  3:131141/1‑87 (MQ=255)
                                                                                                 tggTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCa                                                                                                         >  1:1334357/1‑96 (MQ=255)
                                                                                                 tggTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCa                                                                                                         <  2:1334357/96‑1 (MQ=255)
                                                                                                 tggTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGaa                                                                         <  1:431985/128‑1 (MQ=255)
                                                                                                 tggTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGaaaa                                                    >  4:782120/1‑149 (MQ=255)
                                                                                                     gtACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCa                                                >  4:890859/1‑149 (MQ=255)
                                                                                                      tACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGaa                                                                         >  1:840015/1‑123 (MQ=255)
                                                                                                       aCCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCa                                                                                                         <  1:343015/90‑1 (MQ=255)
                                                                                                        ccGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGgc                                             <  4:780206/149‑1 (MQ=255)
                                                                                                          gTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTaa                                                         <  3:351976/135‑1 (MQ=255)
                                                                                                          gTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGgcgc                                           >  3:983747/1‑149 (MQ=255)
                                                                                                           tGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGaa                                                                         <  4:530209/118‑1 (MQ=255)
                                                                                                           tGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGaa                                                                         >  3:530209/1‑118 (MQ=255)
                                                                                                           tGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCt                                          <  3:573143/149‑1 (MQ=255)
                                                                                                           tGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCt                                          <  1:688101/149‑1 (MQ=255)
                                                                                                           tGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCt                                          >  4:573143/1‑149 (MQ=255)
                                                                                                            gCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTa                                         >  4:309566/1‑149 (MQ=255)
                                                                                                              aTTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGaaaa                                                    <  4:730060/136‑1 (MQ=255)
                                                                                                              aTTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGaaaa                                                    >  3:730060/1‑136 (MQ=255)
                                                                                                               tttATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGc                                                                                                                      >  4:499412/1‑69 (MQ=255)
                                                                                                                  attaCGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGaa                                                                         >  2:255443/1‑111 (MQ=255)
                                                                                                                  attaCGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGaa                                                                         <  1:255443/111‑1 (MQ=255)
                                                                                                                       ggTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGAtt                                                                    <  4:685065/111‑1 (MQ=255)
                                                                                                                           atgaAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAg                                                                                                           >  4:152208/1‑68 (MQ=255)
                                                                                                                             gaAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCa                        <  2:688944/149‑1 (MQ=255)
                                                                                                                               aGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAgg                      >  1:808089/1‑149 (MQ=255)
                                                                                                                                gTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCa                                                                                                         <  2:389826/65‑1 (MQ=255)
                                                                                                                                 tAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCt                            <  3:1468353/141‑1 (MQ=255)
                                                                                                                                 tAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCt                            >  4:1468353/1‑141 (MQ=255)
                                                                                                                                 tAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTc                         >  3:1426898/1‑144 (MQ=255)
                                                                                                                                 tAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTc                         <  4:1426898/144‑1 (MQ=255)
                                                                                                                                 tAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTg                    >  4:371498/1‑149 (MQ=255)
                                                                                                                                  aGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGc                   >  1:1036955/1‑149 (MQ=255)
                                                                                                                                    gACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGaaaa                                                    >  4:472832/1‑114 (MQ=255)
                                                                                                                                      cTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTa               <  3:497984/149‑1 (MQ=255)
                                                                                                                                         gATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCa                                                                                                         >  2:1124246/1‑56 (MQ=255)
                                                                                                                                         gATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCa                                                                                                         <  1:1124246/56‑1 (MQ=255)
                                                                                                                                         gATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGt                                                           <  3:1080032/102‑1 (MQ=255)
                                                                                                                                         gATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGt                                                           >  4:1080032/1‑102 (MQ=255)
                                                                                                                                            ggCAAAACGGATCGTTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTc                                     >  3:1240924/1‑121 (MQ=255)
                                                                                                                                             gCAAAACGGATCGTTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCt                            <  4:824361/129‑1 (MQ=255)
                                                                                                                                             gCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTAc                                        <  2:431004/117‑1 (MQ=255)
                                                                                                                                             gCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGaa        <  1:60207/149‑1 (MQ=255)
                                                                                                                                              cAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCa                                                >  4:985180/1‑108 (MQ=255)
                                                                                                                                               aaaaaCGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCa                                    <  3:683600/119‑1 (MQ=255)
                                                                                                                                               aaaaaCGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTc                         >  1:712717/1‑130 (MQ=255)
                                                                                                                                               aaaaaCGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGa      <  3:309566/149‑1 (MQ=255)
                                                                                                                                                  aaCGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATaa   <  3:553367/149‑1 (MQ=255)
                                                                                                                                                   aCGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAAc  >  1:1376087/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATCAGCTCGACAGCCACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAAC  >  NC_000913/423842‑424137

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