Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,852,559 T→A 9.6% K76* (AAA→TAA) ‡ 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,5590TA 9.6% 268.7 / 14.5 83K76* (AAA→TAA) ‡ilvBacetolactate synthase 2 large subunit
Reads supporting (aligned to +/- strand):  ref base T (38/37);  new base A (3/5);  total (41/42)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CTGGAAGGCGTCGGTGCCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTAATGCC  >  NC_000913/3852416‑3852706
                                                                                                                                               |                                                                                                                                                   
cTGGAAGGCGTCGGTGCCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGt                                                                                                                                                <  2:267890/149‑1 (MQ=255)
 tGGAAGGCGTCGGTGCCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTc                                                                                                                                               >  1:232793/1‑149 (MQ=255)
 tGGAAGGCGTCGGTGCCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTc                                                                                                                                               >  2:33504/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                                                                                                                                            <  3:761880/128‑1 (MQ=255)
                         gAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGACGGTAAACCGGCGGt                                                                                                                                            >  4:761880/1‑128 (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                                                                              >  1:558362/1‑149 (MQ=255)
                                                                  gTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGa                                                                              <  2:614186/149‑1 (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:793440/149‑1 (MQ=255)
                                                                              aTCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCttt                                                                  >  1:257937/1‑149 (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                                                                           <  4:6067/89‑1 (MQ=255)
                                                                                                                                 aCAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGAttt                                                                           >  3:6067/1‑89 (MQ=255)
                                                                                                                                    caCCGACGGTAAACCGGCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTaa                                                           >  4:807239/2‑102 (MQ=255)
                                                                                                                                    caCCGACGGTAAACCGGCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTaa                                                           <  3:807239/101‑1 (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)
                                                                                                                                           ggTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTAAt     <  2:558362/149‑1 (MQ=255)
                                                                                                                                           ggTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTAAt     <  4:388516/149‑1 (MQ=255)
                                                                                                                                           ggTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTAAt     <  3:533026/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)
                                                                                                                                               |                                                                                                                                                   
CTGGAAGGCGTCGGTGCCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTAATGCC  >  NC_000913/3852416‑3852706

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