Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,852,553 C→A 9.1% A78S (GCG→TCG)  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,5530CA 9.1% 255.0 / 12.1 77A78S (GCG→TCG) ilvBacetolactate synthase 2 large subunit
Reads supporting (aligned to +/- strand):  ref base C (36/34);  new base A (3/4);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-01

GTCCACTTCCTGGAAGGCGTCGGTGCCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCT  >  NC_000913/3852407‑3852699
                                                                                                                                                  |                                                                                                                                                  
gtCCACTTCCTGGAAGGCGTCGGTGCCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGAccgccg                                                                                                                                                  <  3:384919/149‑1 (MQ=255)
         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                                                                                                                                     <  1:499825/130‑1 (MQ=255)
                                tcGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGACGGTAAACCGGCGGt                                                                                                                                     >  2:499825/1‑130 (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                                                                       >  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: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                                                    <  1:263211/128‑1 (MQ=255)
                                                                                                                   gCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTaa                                                    >  2:263211/1‑128 (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                                                    <  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)
                                                                                                                                                  |                                                                                                                                                  
GTCCACTTCCTGGAAGGCGTCGGTGCCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCT  >  NC_000913/3852407‑3852699

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