Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,101,014 A→G D85D (GAT→GAC gnd ← 6‑phosphogluconate dehydrogenase, decarboxylating

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,101,0140AG98.7% 269.4 / ‑3.9 76D85D (GAT→GACgnd6‑phosphogluconate dehydrogenase, decarboxylating
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (45/30);  total (46/30)
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.51e-01

CCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAATCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGACAAACTCTTTCACCGTATAGTAAGGAACCAGTTTCTTGCCTGGATTTTCGGCAATCACTTCTTCCGTCTTCT  >  NC_000913/2100869‑2101159
                                                                                                                                                 |                                                                                                                                                 
cccccTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCaa                                                                                                                                               <  2:826896/150‑1 (MQ=255)
           ccGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAgcagca                                                                                                                                       >  2:1061568/1‑147 (MQ=255)
                aGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTGCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCt                                                                                                                              <  1:796117/151‑1 (MQ=255)
                aGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCt                                                                                                                              <  1:558193/151‑1 (MQ=255)
                  aaaCACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGc                                                                                                                            >  1:1081968/1‑151 (MQ=255)
                  aaaCACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGc                                                                                                                            >  2:596603/1‑151 (MQ=255)
                   aaCACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCa                                                                                                                           <  2:217397/151‑1 (MQ=255)
                     cacCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCAcc                                                                                                                         <  2:1153235/151‑1 (MQ=255)
                              cGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATTATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCAc                                                                                                                >  2:731224/1‑151 (MQ=255)
                                aTGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCAcc                                                                                                               >  1:583335/1‑150 (MQ=255)
                                aTGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCa                                                                                                              <  2:175696/151‑1 (MQ=255)
                                 tGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCAt                                                                                                             >  2:237903/1‑151 (MQ=255)
                                     gTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAAc                                                                                                         <  1:181914/151‑1 (MQ=255)
                                        aaaGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGGCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGAGCCTGCACCTGCTTTCACCATTAACAgg                                                                                                      >  2:820470/1‑151 (MQ=255)
                                         aaGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAgg                                                                                                      >  1:730728/1‑150 (MQ=255)
                                         aaGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAAAGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAgg                                                                                                      >  1:479814/1‑150 (MQ=255)
                                           gCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATg                                                                                                   <  1:145448/151‑1 (MQ=255)
                                              ctctGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGcg                                                                                                 >  1:796421/1‑150 (MQ=255)
                                              ctctGCTGAAAGCCCATGATTACGAGGAGTAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGCAAATAGCAGCATCCGCGCCTGCACCTGCTTTCACCATTAAAAGGAAGc                                                                                                  >  1:465196/1‑149 (MQ=255)
                                               tctGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGcga                                                                                                >  2:216419/1‑150 (MQ=255)
                                                ctgctgAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGcga                                                                                                <  2:1212138/149‑1 (MQ=255)
                                                  gctgAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACATGCTTTCACCATTAACAGGATGcgac                                                                                               <  1:1116898/148‑1 (MQ=255)
                                                   ctgAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAgg                                                                                           >  1:138748/1‑151 (MQ=255)
                                                     gAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCg                                                                                         >  2:707832/1‑151 (MQ=255)
                                                     gAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCg                                                                                         >  2:925042/1‑151 (MQ=255)
                                                       aaGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGtt                                                                                       <  1:340644/151‑1 (MQ=255)
                                                       aaGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGtt                                                                                       >  1:217750/1‑151 (MQ=255)
                                                        aGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGttt                                                                                      >  1:359214/1‑151 (MQ=255)
                                                         gCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTc                                                                                     >  1:449027/1‑151 (MQ=255)
                                                           tCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTcc                                                                                    >  2:1080141/1‑150 (MQ=255)
                                                             aCGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCaga                                                                                 <  2:609038/151‑1 (MQ=255)
                                                                aTTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGAt                                                                              <  1:472777/151‑1 (MQ=255)
                                                                  tACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTc                                                                            >  2:1234472/1‑151 (MQ=255)
                                                                  tACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTc                                                                            >  1:650638/1‑151 (MQ=255)
                                                                   acgtcgAATAGTCTCCTGGAAGAATGTGTTACCGCCATCAATCATGACGTCTCCTTTATCGAGATATGGTTTGATTGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCagaaagtc                                                                            <  1:650437/150‑6 (MQ=255)
                                                                   acgacgAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTc                                                                            <  1:1158662/150‑1 (MQ=255)
                                                                   acgacgAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTc                                                                            <  2:256213/150‑1 (MQ=255)
                                                                    cgacgaATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCg                                                                           >  2:1239631/1‑150 (MQ=255)
                                                                      acgaATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGAc                                                                         >  1:1072033/1‑150 (MQ=255)
                                                                          aTAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGACAAAct                                                                    >  2:1190266/1‑151 (MQ=255)
                                                                           tAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGACAAActc                                                                   >  1:295273/1‑151 (MQ=255)
                                                                            aGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGACAAActct                                                                  >  1:142259/1‑151 (MQ=255)
                                                                             gtgtCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGACAAACTCtt                                                                 <  1:1092533/151‑1 (MQ=255)
                                                                                 ccTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGACAAACTCTTTCAc                                                             >  2:203644/1‑151 (MQ=255)
                                                                                 ccTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGACAAACTCTTTCAc                                                             <  1:975906/151‑1 (MQ=255)
                                                                                     gaagaaGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGACAAACTCTTTCACCGta                                                         <  2:989935/151‑1 (MQ=255)
                                                                                     gaagaaGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGACAAACTCTTTCACCGta                                                         <  1:207559/151‑1 (MQ=255)
                                                                                             tgtTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGACAAACACTTTCACCGTATAGTaac                                                  <  1:690511/150‑2 (MQ=255)
                                                                                                  ccaccaTCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGACAAACTCTTTCACCGTATAGTAAGGAAcc                                             >  2:269892/1‑150 (MQ=255)
                                                                                                  ccaccaTCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGACAAACTCTTTCACCGTATAGTAAGGAACCa                                            <  1:538759/151‑1 (MQ=255)
                                                                                                   caccaTCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGACAAACTCTTTCACCGTATAGTAAGGAACCAg                                           <  1:463585/151‑1 (MQ=255)
                                                                                                    accaTCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGACAAACTCTTTCACCGTATAGTAAGGAACCAGt                                          <  2:295289/151‑1 (MQ=255)
                                                                                                    accaTCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGACAAACTCTTTCACCGTATAGTAAGGAACCAGt                                          >  2:657625/1‑151 (MQ=255)
                                                                                                    accaTCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGACAAACTCTTTCACCGTATAGTAAGGAACCAGt                                          >  1:853052/1‑151 (MQ=255)
                                                                                                     ccaTCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGACAAACTCTTTCACCGTATAGTAAGGAACCAGtt                                         >  2:234746/1‑151 (MQ=255)
                                                                                                      caTCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGACAAACTCTTTCACCGTATAGTAAGGAACCAGttt                                        >  1:934015/1‑151 (MQ=255)
                                                                                                      caTCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGACAAACTCTTTCACCGTATAGTAAGGAACCAGttt                                        <  2:653206/151‑1 (MQ=255)
                                                                                                       aTCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGACAAACTCTTTCACCGTATAGTAAGGAACCAGTTTc                                       >  1:790972/1‑151 (MQ=255)
                                                                                                          aaTGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGACAAACTCTTTCACCGTATAGTAAGGAACCAGTTTCTTg                                    >  1:315878/1‑151 (MQ=255)
                                                                                                                   gTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGACAAACTCTTTCACCGTATAGTAAGGAACCAGTTTCTTGCCTGGAttt                           <  2:106066/151‑1 (MQ=255)
                                                                                                                   gTCTCCTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGACAAACTCTTTCACCGTATAGTAAGGAACCAGTTTCTTGCCTGGAtt                            <  2:783993/150‑1 (MQ=255)
                                                                                                                       ccTTTATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGACAAACTCTTTCACCGTATAGTAAGGAACCAGTTTCTTGCCTGGAtttt                          <  2:70899/148‑1 (MQ=255)
                                                                                                                         tttATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGACAAACTCTTTCACCGTATAGTAAGGAACCAGTTTCTTGCCTGGATTTTCGGc                      >  1:843146/1‑150 (MQ=255)
                                                                                                                         tttATCGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGACAAACTCTTTCACCGTATAGTAAGGAACCAGTTTCTTGCCTGGATTTTCGGc                      >  2:1234002/1‑150 (MQ=255)
                                                                                                                              cGAGATATGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGACAAACTCTTTCACCGTATAGTAAGGAACCAGTTTCTTGCCTGGATTTTCGGCAATCAc                >  2:382848/1‑151 (MQ=255)
                                                                                                                                  atatGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGACAAACTCTTTCACCGTATAGTAAGGAACCAGTTTCTTGCCTGGATTTTCGGCAATCActtct            >  1:832963/1‑151 (MQ=255)
                                                                                                                                  atatGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGACAAACTCTTTCACCGTATAGTAAGGAACCAGTTTCTTGCCTGGATTTTCGGCAATCActtct            >  1:96188/1‑151 (MQ=255)
                                                                                                                                  atatGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCCGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGACAAACTCTTTCACCGTATAGTAAGGAACCAGTTTCTTGCCTGGATTTTCGGCAATCActtct            >  2:208167/1‑151 (MQ=255)
                                                                                                                                    atGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGACAAACTCTTTCACCGTATAGTAAGGAACCAGTTTCTTGCCTGGATTTTCGGCAATCACttcttc          >  1:840727/1‑151 (MQ=255)
                                                                                                                                     tGGTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGACAAACTCTTTCACCGTATAGTAAGGAACCAGTTTCTTGCCTGGATTTTCGGCAATCACTTCTTcc         >  1:384334/1‑151 (MQ=255)
                                                                                                                                       gTTTGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGACAAACTCTTTCACCGTATAGTAAGGAACCAGTTTCTTGCCTGGATTTTCGGCAATCACTTCTTCCGt       <  2:664401/151‑1 (MQ=255)
                                                                                                                                         ttGAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGACAAACTCTTTCACCGTATAGTAAGGAACCAGTTTCTTGCCTGGATTTTCGGCAATCACTTCTTCCGtct     >  2:320727/1‑151 (MQ=255)
                                                                                                                                         ttGAGGGAATCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGACAAACTCTTTCACCGTATAGTAAGGAACCAGTTTCTTGCCTGGATTTTCGGCAATCACTTCTTCCGtct     >  2:1139278/1‑151 (MQ=255)
                                                                                                                                           gAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGACAAACTCTTTCACCGTATAGTAAGGAACCAGTTTCTTGCCTGGATTTTCGGCAATCACTTCTTCCGtcttc   <  2:219518/151‑1 (MQ=255)
                                                                                                                                           gAGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGACAAACTCTTTCACCGTATAGTAAGGAACCAGTTTCTTGCCTGGATTTTCGGCAATCACTTCTTCCGtcttc   >  2:107638/1‑151 (MQ=255)
                                                                                                                                            aGGGAGTCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGACAAACTCTTTCACCGTATAGTAAGGAACCAGTTTCTTGCCTGGATTTTCGGCAATCACTTCTTCCGtcttct  >  1:28607/1‑151 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CCCCCTCTTCACCGCCAGAAACACCGGTACCGATGAAGTTAAAGCCCTCTGCTGAAAGCTCACGATTACGACGAATAGTGTCCTGGAAGAAGGTGTTACCACCATCAATGATGATGTCTCCTTTATCGAGATATGGTTTGAGGGAATCAATAGCAGCATCCGTGCCTGCACCTGCTTTCACCATTAACAGGATGCGACGAGGCGTTTCCAGAGATTCGACAAACTCTTTCACCGTATAGTAAGGAACCAGTTTCTTGCCTGGATTTTCGGCAATCACTTCTTCCGTCTTCT  >  NC_000913/2100869‑2101159

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