Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,852,559 T→A 5.3% 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 5.3% 428.5 / 8.8 114K76* (AAA→TAA) ‡ilvBacetolactate synthase 2 large subunit
Reads supporting (aligned to +/- strand):  ref base T (52/56);  new base A (3/3);  total (55/59)
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.85e-01

Only 100 of 114 total aligned reads displayed.

GGAAGGCGTCGGTGCCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGCCTG  >  NC_000913/3852418‑3852690
                                                                                                                                             |                                                                                                                                   
ggAAGGCGTCGGTGCCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTc                                                                                                                               >  3:880701/1‑148 (MQ=255)
  aaGGCGTCGGTGCCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCgg                                                                                                                             <  3:480952/148‑1 (MQ=255)
  aaGGCGTCGGTGCCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCgg                                                                                                                             <  1:707266/148‑1 (MQ=255)
    ggCGTCGGTGCCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTg                                                                                                                           >  2:49092/1‑148 (MQ=255)
         cGGTGCCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTgcgc                                                                                                                        >  3:994074/1‑146 (MQ=255)
               cGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTcc                                                                                                                <  3:130051/148‑1 (MQ=255)
                 atcatcGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCt                                                                                                              >  1:410563/1‑148 (MQ=255)
                 atcatcGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCt                                                                                                              >  3:715346/1‑148 (MQ=255)
                 atcatcGAGGCAGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTgcgcg                                                                                                                       <  2:1128558/139‑1 (MQ=255)
                 atcatcGAGGCAGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTgcgcg                                                                                                                       >  1:1128558/1‑139 (MQ=255)
                  tcatcGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCACGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTg                                                                                                             >  1:984882/1‑148 (MQ=255)
                   catcGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGa                                                                                                            <  1:49092/148‑1 (MQ=255)
                     tcGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGc                                                                                                          >  1:943165/1‑148 (MQ=255)
                      cGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCg                                                                                                         <  4:403769/148‑1 (MQ=255)
                       gAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGa                                                                                                        <  2:410563/148‑1 (MQ=255)
                       gAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGa                                                                                                        <  4:715346/148‑1 (MQ=255)
                        aGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCgg                                                                                                                             >  4:592393/1‑126 (MQ=255)
                        aGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCgg                                                                                                                             <  3:592393/126‑1 (MQ=255)
                            gggAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTc                                                                                                                               >  4:77743/1‑120 (MQ=255)
                            gggAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTc                                                                                                                               <  3:77743/120‑1 (MQ=255)
                             ggAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGc                                                                                                  >  4:79590/1‑148 (MQ=255)
                              gAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGcc                                                                                                 >  4:916000/1‑148 (MQ=255)
                              gAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGcc                                                                                                 <  3:744640/148‑1 (MQ=255)
                                      ccAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAg                                                                                                           <  3:320422/130‑1 (MQ=255)
                                       cAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTa                                                                                                                                    <  2:217051/104‑1 (MQ=255)
                                       cAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTa                                                                                                                                    >  1:217051/1‑104 (MQ=255)
                                       cAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTg                                                                                        <  2:984882/148‑1 (MQ=255)
                                            aTGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCAt                                                                                   <  3:79590/148‑1 (MQ=255)
                                                  aTCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACggg                                                                             <  4:49721/148‑1 (MQ=255)
                                                  aTCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACg                                                                               >  3:955770/1‑146 (MQ=255)
                                                  aTCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACg                                                                               <  4:955770/146‑1 (MQ=255)
                                                   tCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGTGATAAAGCCCGCGCCCTGTTCATGACGGGc                                                                            >  1:876573/1‑148 (MQ=255)
                                                       cGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACg                                                                               >  1:177820/1‑141 (MQ=255)
                                                       cGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACg                                                                               <  2:177820/141‑1 (MQ=255)
                                                           aTGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAatat                                                                    >  3:672685/1‑148 (MQ=255)
                                                                gTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTc                                                                                                                 >  2:895861/1‑98 (MQ=255)
                                                                gTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTc                                                                                                                 <  1:895861/98‑1 (MQ=255)
                                                                 tCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCa                                                                                    <  3:485763/126‑1 (MQ=255)
                                                                 tCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCa                                                                                    >  4:485763/1‑126 (MQ=255)
                                                                 tCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGc                                                                            >  3:1031370/1‑134 (MQ=255)
                                                                 tCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGc                                                                            <  4:1031370/134‑1 (MQ=255)
                                                                    aGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATg                                                                   <  1:358005/140‑1 (MQ=255)
                                                                    aGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATg                                                                   >  2:358005/1‑140 (MQ=255)
                                                                    aGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGAttt                                                           >  3:954367/1‑148 (MQ=255)
                                                                     gCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTa                                                                                                                                    <  4:194165/74‑1 (MQ=255)
                                                                     gCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTa                                                                                                                                    >  3:194165/1‑74 (MQ=255)
                                                                     gCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTg                                                          >  4:1038400/1‑148 (MQ=255)
                                                                     gCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTg                                                          <  2:74060/148‑1 (MQ=255)
                                                                      ccGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTg                                                                                                             >  4:926237/1‑96 (MQ=255)
                                                                       cgcgCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCg                                                        <  3:1038400/148‑1 (MQ=255)
                                                                           cATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTg                                                                                        <  1:150043/112‑1 (MQ=255)
                                                                           cATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTg                                                                                        >  2:150043/1‑112 (MQ=255)
                                                                             tCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCttt                                                  <  2:943165/148‑1 (MQ=255)
                                                                               ggCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTgg                                                <  4:954367/148‑1 (MQ=255)
                                                                                gCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGa                                                                        >  1:845142/1‑123 (MQ=255)
                                                                                gCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGa                                                                        <  2:845142/123‑1 (MQ=255)
                                                                                   aTGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGa                                                                        <  4:623646/120‑1 (MQ=255)
                                                                                   aTGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGa                                                                        >  3:623646/1‑120 (MQ=255)
                                                                                    tGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTaa                                           <  3:916000/148‑1 (MQ=255)
                                                                                       cGGTCACCAGGTTAGTCGCAGCCGGTCCGCTACAGGCCATATAGACCGACGGTAAACCGGCGGt                                                                                                                            <  4:945079/64‑1 (MQ=25)
                                                                                       cGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGAtt                                                            >  1:999823/1‑128 (MQ=255)
                                                                                       cGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGAtt                                                            <  2:999823/128‑1 (MQ=255)
                                                                                       cGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGACGGTAAACCGGCGGt                                                                                                                            >  3:945079/1‑64 (MQ=37)
                                                                                         gTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGa                                                                        >  1:1091930/1‑114 (MQ=255)
                                                                                         gTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGa                                                                        <  2:1091930/114‑1 (MQ=255)
                                                                                                  ttAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAAcag                             >  4:483848/1‑148 (MQ=255)
                                                                                                   tAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGCGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAAcagg                            >  2:1111176/1‑148 (MQ=255)
                                                                                                       cGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAAcaggcagg                        <  4:694694/148‑1 (MQ=255)
                                                                                                         cACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTgg                                                >  1:1008367/1‑122 (MQ=255)
                                                                                                           cccGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCgcgc                                                                                            >  1:568940/1‑76 (MQ=255)
                                                                                                           cccGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCgcgc                                                                                            <  2:568940/76‑1 (MQ=255)
                                                                                                           cccGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGa                       <  4:1026636/145‑1 (MQ=255)
                                                                                                            ccGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGa                   <  3:483848/148‑1 (MQ=255)
                                                                                                                 ccGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGccc                                                                                          >  4:1013401/1‑72 (MQ=255)
                                                                                                                 ccGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGccc                                                                                          <  3:1013401/72‑1 (MQ=255)
                                                                                                                 ccGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGc              >  4:858346/1‑148 (MQ=255)
                                                                                                                   gCTACAGGCCATAGAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCAt                                      >  3:82921/1‑122 (MQ=255)
                                                                                                                   gCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCAt                                      <  4:82921/122‑1 (MQ=255)
                                                                                                                   gCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGccc            <  4:672685/148‑1 (MQ=255)
                                                                                                                    cTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAAcaggc                           >  2:9206/1‑132 (MQ=255)
                                                                                                                    cTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAAcaggc                           <  1:9206/132‑1 (MQ=255)
                                                                                                                    cTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGa                   <  3:1011523/140‑1 (MQ=255)
                                                                                                                          gCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGc     >  1:225947/1‑148 (MQ=255)
                                                                                                                          gCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGc     >  3:968759/1‑148 (MQ=255)
                                                                                                                           ccATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTg                                                                                                             >  4:314636/1‑43 (MQ=255)
                                                                                                                           ccATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTg                                                                                                             <  3:314636/43‑1 (MQ=255)
                                                                                                                           ccATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCAt                                      >  3:1079830/1‑114 (MQ=255)
                                                                                                                           ccATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCAt                                      <  4:1079830/114‑1 (MQ=255)
                                                                                                                           ccATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAAcaggc                           <  1:315767/125‑1 (MQ=255)
                                                                                                                             aTACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGaa                  <  2:230152/132‑1 (MQ=255)
                                                                                                                             aTACAGACCGCCGGTTTACCGTCGGCGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGCCTg  <  1:1111176/148‑1 (MQ=255)
                                                                                                                              tACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACg                                                                               >  4:612135/1‑70 (MQ=255)
                                                                                                                              tACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACg                                                                               <  3:612135/70‑1 (MQ=255)
                                                                                                                                cgcaCCGACGGTAAACCGGCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTgg                                                >  2:363907/4‑99 (MQ=255)
                                                                                                                                cgcaCCGACGGTAAACCGGCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTgg                                                <  1:363907/96‑1 (MQ=255)
                                                                                                                                    ccgacgGTAAACCGGCGGTGCGCGCCATTCCCTGAg                                                                                                           <  3:945147/36‑1 (MQ=9)
                                                                                                                                    ccgacgGTAAACCGGCGGTGCGCGCCATTCCCTGAg                                                                                                           >  4:945147/1‑36 (MQ=12)
                                                                                                                                      gccgGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGc                                                     <  4:962918/88‑1 (MQ=255)
                                                                                                                                      gccgGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGc                                                     >  3:962918/1‑88 (MQ=255)
                                                                                                                                       ccgGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTTGCTTAAGGCa                                       <  3:669190/101‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                   
GGAAGGCGTCGGTGCCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGTCGCACCCGGTCCGCTACAGGCCATACAGACCGCCGGTTTACCGTCGGTGCGCGCCATTCCCTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGCCTG  >  NC_000913/3852418‑3852690

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