Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,852,564 T→G 9.1% D74A (GAC→GCC)  ilvB ← acetolactate synthase 2 large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,852,5640TG 9.1% 273.1 / 17.4 89D74A (GAC→GCC) ilvBacetolactate synthase 2 large subunit
Reads supporting (aligned to +/- strand):  ref base T (40/40);  new base G (3/5);  total (44/45)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

CTGGAAGGCGTCGGTGCCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTAATGCCCTGC  >  NC_000913/3852416‑3852710
                                                                                                                                                    |                                                                                                                                                  
cTGGAAGGCGTCGGTGCCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGt                                                                                                                                                    <  2:267890/149‑1 (MQ=255)
 tGGAAGGCGTCGGTGCCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTc                                                                                                                                                   >  2:33504/1‑149 (MQ=255)
 tGGAAGGCGTCGGTGCCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTc                                                                                                                                                   >  1:232793/1‑149 (MQ=255)
 tGGAAGGCGTCGGTGCCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGAGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTc                                                                                                                                                   >  3:196433/1‑149 (MQ=255)
    aaGGCGTCGGTGCCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGACGGTAAACCGGCGGt                                                                                                                                                <  3:686993/149‑1 (MQ=255)
      ggCGTCGGTGCCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTgc                                                                                                                                              >  3:531237/1‑149 (MQ=255)
       gCGTCGGTGCCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTgcg                                                                                                                                             >  1:441863/1‑149 (MQ=255)
                ccGATCAACGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTcc                                                                                                                                    <  4:453485/149‑1 (MQ=255)
                 cGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTccc                                                                                                                                   <  1:435559/149‑1 (MQ=255)
                    tcatcGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGa                                                                                                                                <  2:525000/149‑1 (MQ=255)
                      atcGAGGCGGGAACCTTACCAGTGATGCAAATCGGCGGGAGGGAGTCCACCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGc                                                                                                                              >  2:37281/1‑149 (MQ=255)
                       tcGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGACGGTAAACCGGCGGt                                                                                                                                                >  2:499825/1‑130 (MQ=255)
                       tcGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGACGGTAAACCGGCGGt                                                                                                                                                <  1:499825/130‑1 (MQ=255)
                         gAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGGGCGCGCCATTCCCTGAGCGAt                                                                                                                           >  3:71003/1‑149 (MQ=255)
                         gAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGACGGTAAACCGGCGGt                                                                                                                                                >  4:761880/1‑128 (MQ=255)
                         gAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGACGGTAAACCGGCGGt                                                                                                                                                <  3:761880/128‑1 (MQ=255)
                              gggAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGc                                                                                                                      <  3:797771/149‑1 (MQ=255)
                                  aCCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATaaa                                                                                                                        <  1:174369/143‑1 (MQ=255)
                                  aCCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATaaa                                                                                                                        >  2:174369/1‑143 (MQ=255)
                                  aCCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCcgc                                                                                                                  >  3:388516/1‑149 (MQ=255)
                                       aCCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCt                                                                                                             >  4:680867/1‑149 (MQ=255)
                                                gCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCAt                                                                                                       >  1:388326/1‑146 (MQ=255)
                                                gCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCAt                                                                                                       <  2:388326/146‑1 (MQ=255)
                                                  aaaTCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACgg                                                                                                  <  2:441863/149‑1 (MQ=255)
                                                   aaTCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACggg                                                                                                 >  3:617053/1‑149 (MQ=255)
                                                          gggATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAat                                                                                          >  4:533026/1‑149 (MQ=255)
                                                            gATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAatat                                                                                        >  3:254058/1‑149 (MQ=255)
                                                             aTGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATg                                                                                       <  3:435589/149‑1 (MQ=255)
                                                               ggAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGc                                                                                     <  4:254058/149‑1 (MQ=255)
                                                                gggTCCAGCCGCGCATCGGCAATGGCGGCCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCg                                                                                    <  2:332981/147‑1 (MQ=255)
                                                                  gTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGa                                                                                  <  2:614186/149‑1 (MQ=255)
                                                                  gTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGa                                                                                  >  1:558362/1‑149 (MQ=255)
                                                                   tCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGAt                                                                                 >  4:626760/1‑149 (MQ=255)
                                                                   tCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGAt                                                                                 >  4:31580/1‑149 (MQ=255)
                                                                      aGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTg                                                                              >  4:431436/1‑149 (MQ=255)
                                                                       gCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGc                                                                             <  2:43019/149‑1 (MQ=255)
                                                                             cATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCTGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCtt                                                                       <  4:576067/149‑1 (MQ=255)
                                                                              aTCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCttt                                                                      >  1:257937/1‑149 (MQ=255)
                                                                              aTCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCttt                                                                      <  1:793440/149‑1 (MQ=255)
                                                                               tCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGATTACCGACGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTg                                                                     >  4:402408/1‑149 (MQ=255)
                                                                                  gCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGc                                                                                     >  4:573692/1‑130 (MQ=255)
                                                                                  gCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGc                                                                                     <  3:573692/130‑1 (MQ=255)
                                                                                    aaTGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTa                                                                >  3:9569/1‑149 (MQ=255)
                                                                                     aTGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTaa                                                               <  3:388467/149‑1 (MQ=255)
                                                                                     aTGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTaa                                                               <  4:531237/149‑1 (MQ=255)
                                                                                      tGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAg                                                              <  2:305133/149‑1 (MQ=255)
                                                                                           gTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATc                                                         <  4:186706/149‑1 (MQ=255)
                                                                                              aCCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAATGCATCGTa                                                      >  2:771130/1‑149 (MQ=255)
                                                                                               ccAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTaa                                                     >  1:503683/1‑149 (MQ=255)
                                                                                                cAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTaaa                                                    <  1:404368/149‑1 (MQ=255)
                                                                                                cAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTAACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTaac                                                    <  3:877684/149‑2 (MQ=255)
                                                                                                 aGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTCAAc                                                   <  1:771130/149‑1 (MQ=255)
                                                                                                 aGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAAc                                                   >  1:27950/1‑149 (MQ=255)
                                                                                                  ggTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAAca                                                  <  3:664876/149‑1 (MQ=255)
                                                                                                   gTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAAcag                                                 <  4:682278/149‑1 (MQ=255)
                                                                                                      aGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAAcaggca                                              >  3:43914/1‑149 (MQ=255)
                                                                                                        tCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAAcaggcagg                                            >  2:337128/1‑149 (MQ=255)
                                                                                                         cGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGa                                           <  3:287061/149‑1 (MQ=255)
                                                                                                          gCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTaa                                                               >  2:263211/1‑128 (MQ=255)
                                                                                                          gCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTaa                                                               <  1:263211/128‑1 (MQ=255)
                                                                                                           cACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATa                                         >  3:17922/1‑149 (MQ=255)
                                                                                                            aCCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCAt                                                                                                       >  1:725001/1‑86 (MQ=255)
                                                                                                            aCCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCAt                                                                                                       <  2:725001/86‑1 (MQ=255)
                                                                                                              ccGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGaa                                      <  4:68391/149‑1 (MQ=255)
                                                                                                                      cTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCgg                              >  1:863202/1‑149 (MQ=255)
                                                                                                                        aCAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGaa                            <  3:738417/149‑1 (MQ=255)
                                                                                                                         cAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAAt                           <  1:558363/149‑1 (MQ=255)
                                                                                                                          aGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGa                             >  2:235149/1‑146 (MQ=255)
                                                                                                                          aGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGa                             <  1:235149/146‑1 (MQ=255)
                                                                                                                            gCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGcc                        >  3:432733/1‑149 (MQ=255)
                                                                                                                             ccATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGCCt                       <  1:151044/149‑1 (MQ=255)
                                                                                                                                tACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTc                    >  2:646797/1‑149 (MQ=255)
                                                                                                                                 aCAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGAttt                                                                               >  3:6067/1‑89 (MQ=255)
                                                                                                                                 aCAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGAttt                                                                               <  4:6067/89‑1 (MQ=255)
                                                                                                                                    caCCGACGGTAAACCGGCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTaa                                                               <  3:807239/101‑1 (MQ=255)
                                                                                                                                    caCCGACGGTAAACCGGCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTaa                                                               >  4:807239/2‑102 (MQ=255)
                                                                                                                                     aCCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAAt               >  1:386175/1‑149 (MQ=255)
                                                                                                                                      ccgccgGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATc              <  2:712613/149‑1 (MQ=255)
                                                                                                                                       cgccgGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCt             >  4:809878/1‑149 (MQ=255)
                                                                                                                                          cgGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCCCAATCTTaa          >  1:78650/1‑149 (MQ=255)
                                                                                                                                           ggtaaaCCGGCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTAAt         <  4:614523/144‑1 (MQ=255)
                                                                                                                                           ggTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTAAt         <  2:558362/149‑1 (MQ=255)
                                                                                                                                           ggTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTAAt         <  3:533026/149‑1 (MQ=255)
                                                                                                                                           ggTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTAAt         <  4:388516/149‑1 (MQ=255)
                                                                                                                                              ttACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAAcag                                                 >  2:256809/1‑106 (MQ=255)
                                                                                                                                              ttACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAAcag                                                 <  1:256809/106‑1 (MQ=255)
                                                                                                                                              ttACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTAATGcc      <  4:43914/149‑1 (MQ=255)
                                                                                                                                                aCCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTAATGCCct    <  3:31580/149‑1 (MQ=255)
                                                                                                                                                  cGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTAATGCCctgc  >  4:284053/1‑149 (MQ=255)
                                                                                                                                                  cGTCGGTGCGCGCCATTCCCTGAGCAATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTAATGCCctgc  >  1:133557/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CTGGAAGGCGTCGGTGCCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTAATGCCCTGC  >  NC_000913/3852416‑3852710

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑