Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,956,462 T→A 5.1% C378S (TGC→AGC)  ilvA → l‑threonine dehydratase, biosynthetic; also known as threonine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,956,4620TA 5.1% 263.5 / 2.9 78C378S (TGC→AGC) ilvAl‑threonine dehydratase, biosynthetic; also known as threonine deaminase
Reads supporting (aligned to +/- strand):  ref base T (37/37);  new base A (2/2);  total (39/39)
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.88e-01

GAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGGAAGAAAAAGGCAGCTTCCTCAAATTCTGCCAACTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTATATGGTCGGCG  >  NC_000913/3956315‑3956603
                                                                                                                                                   |                                                                                                                                             
gAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGGAAGAAAAAGGCAGCTTCCTCAAATTCTGCCAACTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCt                                                                                                                                               >  8:151218/1‑148 (MQ=255)
        cGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGGAAGAAAAAGGCAGCTTCCTCAAATTCTGCCAACTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCt                                                                                                                                         <  7:323261/146‑1 (MQ=255)
        cGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGGAAGAAAAAGGCAGCTTCCTCAAATTCTGCCAACTGCTTGGCGGGCGTTCGGGCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCtt                                                                                                                                        <  3:9610/147‑1 (MQ=255)
         gAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGGAAGAAAAAGGCAGCTTCCTCAAATTCTGCCAACTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTg                                                                                                                                      <  7:151218/148‑1 (MQ=255)
         gAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGGAAGAAAAAGGCAGCTTCCTCAAATTCTGCCAACTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTg                                                                                                                                      >  6:118160/1‑148 (MQ=255)
           aCTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGGAAGAAAAAGGCAGCTTCCTCAAATTCTGCCAACTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTc                                                                                                                                    <  5:118160/148‑1 (MQ=255)
                      aGCGTGAAGCGTTGTTGGCGGTGACCATTCCGGAAGAAAAAGGCAGCTTCCTCAAATTCTGCCAACTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCt                                                                                                                         >  3:356239/1‑148 (MQ=255)
                         gTGAAGCGTTGTTGGCGGTGACCATTCCGGAAGAAAAAGGCAGCTTCCTCAAATTCTGCCAACTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAg                                                                                                                      <  4:356239/148‑1 (MQ=255)
                           gAAGCGTTGTTGGCGGTGACCATTCCGGAAGAAAAAGGCAGCTTCCTCAAATTCTGCCAACTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGcc                                                                                                                    >  2:359553/1‑148 (MQ=255)
                            aaGCGTTGTTGGCGGTGACCATTCCGGAAGAAAAAGGCAGCTTCCTCAAATTCTGCCAACTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCcg                                                                                                                   >  6:236442/1‑148 (MQ=255)
                                gttgttGGCGGTGACCATTCCGGAAGAAAAAGGCAGCTTCCTCAAATTCTGCCAACTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGc                                                                                                               <  5:313515/148‑1 (MQ=255)
                                  tgttgGCGGTGACCATTCCGGAAGAAAAAGGCAGCTTCCTCAAATTCTGCCAACTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCt                                                                                                             >  6:60122/1‑148 (MQ=255)
                                           tGACCATTCCGGAAGAAAAAGGCAGCTTCCTCAAATTCTGCCAACTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAgcg                                                                                                    >  2:123627/1‑148 (MQ=255)
                                                aTTCCGGAAGAAAAAGGCAGCTTCCTCAAATTCTGCCAACTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCa                                                                                                  >  1:217432/1‑145 (MQ=255)
                                                  tCCGGAAGAAAAAGGCAGCTTCCTCAAATTCTGCCAACTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGa                                                                                              >  6:288566/1‑147 (MQ=255)
                                                          aaaaaGGCAGCTTCCTCAAATTCTGCCAACTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCa                                                                                     <  1:123627/148‑1 (MQ=255)
                                                           aaaaGGCAGCTTCCTCAAATTCTGCCAACTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAg                                                                                    <  5:60122/148‑1 (MQ=255)
                                                            aaaGGCAGCTTCCTCAAATTCTGCCAACTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCt                                                                                                                         <  6:127381/110‑1 (MQ=255)
                                                            aaaGGCAGCTTCCTCAAATTCTGCCAACTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCt                                                                                                                         >  5:127381/1‑110 (MQ=255)
                                                             aaGGCAGCTTCCTCAAATTCTGCCAACTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAAtttaacag                                                                                    >  7:378185/1‑141 (MQ=255)
                                                                gCAGCTTCCTCAAATTCTGCCAACTCCTTGGCGGGCGTGAGGTCACCGAGTTCAACTACCGGTTTCCCGATGCCAAAAACGCCTGCAGCTTTGTCGGGGCGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCt                                                                               <  4:83317/148‑1 (MQ=255)
                                                                  aGCTTCCTCAAATTCTGCCAACTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCa                                                                             <  5:236442/148‑1 (MQ=255)
                                                                  aGCTTCCTCAAATTCTGCCAACTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCa                                                                             <  6:297818/148‑1 (MQ=255)
                                                                    cTTCCTCAAATTCTGCCAACTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTc                                                                                                                                    >  8:193046/1‑91 (MQ=255)
                                                                    cTTCCTCAAATTCTGCCAACTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTc                                                                                                                                    <  7:193046/91‑1 (MQ=255)
                                                                    cTTCCTCAAATTCTGCCAACTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCAGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAg                                                                                    >  5:254557/1‑139 (MQ=255)
                                                                    cTTCCTCAAATTCTGCCAACTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCAGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAg                                                                                    >  5:249016/1‑139 (MQ=255)
                                                                    cTTCCTCAAATTCTGCCAACTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCAGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAg                                                                                    <  6:249016/139‑1 (MQ=255)
                                                                    cTTCCTCAAATTCTGCCAACTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCAGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAg                                                                                    <  6:254557/139‑1 (MQ=255)
                                                                        cTCAAATTCTGCCAACTGCTTGGTGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGc                                                                                      <  4:82614/133‑1 (MQ=255)
                                                                        cTCAAATTCTGCCAACTGCTTGGTGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGc                                                                                      >  3:82614/1‑133 (MQ=255)
                                                                        cTCAAATTCTGCCAACTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAacgacg                                                                       >  2:11749/1‑148 (MQ=255)
                                                                          cAAATTCTGCCAACTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGAcggc                                                                     >  4:224395/1‑148 (MQ=255)
                                                                           aaaTTCTGCCAACTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCTTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCCACGActgcg                                                                    >  6:294544/1‑148 (MQ=255)
                                                                                   ccAACTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTAc                                                               <  3:224395/145‑1 (MQ=255)
                                                                                   ccAACTGCTTGGCGGGCGTTCGGGCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGc                                                            <  5:288566/148‑1 (MQ=255)
                                                                                      aCTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTg                                                         >  8:14663/1‑148 (MQ=255)
                                                                                        tGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGt                                                       >  5:174745/1‑148 (MQ=255)
                                                                                            tGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGa                                                    >  6:304/1‑147 (MQ=255)
                                                                                             ggcgggcgTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGAtc                                                  <  6:174745/148‑1 (MQ=255)
                                                                                                      tCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGa                                                                                   >  3:257007/1‑106 (MQ=255)
                                                                                                      tCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGa                                                                                   <  4:257007/106‑1 (MQ=255)
                                                                                                      tCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCcgac                                          >  1:320964/1‑147 (MQ=255)
                                                                                                      tCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCcgac                                          >  4:228443/1‑147 (MQ=255)
                                                                                                      tCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCAGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGt                                                       >  4:196977/1‑134 (MQ=255)
                                                                                                      tCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCGGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGt                                                       <  3:196977/134‑1 (MQ=255)
                                                                                                        ggTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCcgacgac                                       <  1:359553/148‑1 (MQ=255)
                                                                                                          tCACCGAGTCCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCcgacgacga                                     <  5:269857/148‑1 (MQ=255)
                                                                                                           cACCGAGTCCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCcgacgacga                                     >  6:269857/1‑147 (MQ=255)
                                                                                                            aCCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGaaa                                   >  3:135033/1‑148 (MQ=255)
                                                                                                            aCCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGaa                                    <  2:217432/147‑1 (MQ=255)
                                                                                                            aCCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGaa                                    <  1:11749/147‑1 (MQ=255)
                                                                                                             ccGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCa                                                                                                  >  2:246580/1‑84 (MQ=255)
                                                                                                             ccGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCa                                                                                                  <  1:246580/84‑1 (MQ=255)
                                                                                                              cGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAAt                                  <  6:171865/147‑1 (MQ=255)
                                                                                                              cGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAAt                                  >  5:171865/1‑147 (MQ=255)
                                                                                                                 gTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAATGGCg                              <  5:304/148‑1 (MQ=255)
                                                                                                                    ccaCTACCGTTTTGCCGATGCCAAAAACGCCTGCATCGTTGTCGGTGTGCGCCTGAGCCGCGGCCGCGACGAGCGCAAAGAAATTTTGCAGATGCTCACCGACGGCGGCTACAGCGTGGGTGATCTCTCCGACGACGAAATGGCGaa                            <  5:294544/145‑1 (MQ=255)
                                                                                                                    cAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAATGGCGaa                            <  3:228443/147‑1 (MQ=255)
                                                                                                                      aCTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGAt                                                                                  <  7:130317/91‑1 (MQ=255)
                                                                                                                      aCTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGAt                                                                                  >  8:130317/1‑91 (MQ=255)
                                                                                                                       cTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTa                        >  3:2884/1‑148 (MQ=255)
                                                                                                                         aCCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGa                                                    <  5:200957/118‑1 (MQ=255)
                                                                                                                         aCCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGa                                                    >  6:200957/1‑118 (MQ=255)
                                                                                                                         aCCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGaa                                    <  6:147193/134‑1 (MQ=255)
                                                                                                                         aCCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGaa                                    >  5:147193/1‑134 (MQ=255)
                                                                                                                             ttttGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTg                  <  1:99555/148‑1 (MQ=255)
                                                                                                                              tttGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCcgacgacga                                     >  8:177746/1‑128 (MQ=255)
                                                                                                                              tttGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCcgacgacga                                     <  7:177746/128‑1 (MQ=255)
                                                                                                                              tttGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGaa                                    >  5:301116/1‑129 (MQ=255)
                                                                                                                              tttGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGaa                                    <  6:301116/129‑1 (MQ=255)
                                                                                                                                  ccGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAACAGCGCAATGAAATTTTGCAGATGCTCAACGACGGCTGCTACAGCGTGGTTGATCTCTCCGACGACGAAATGTCGAAGCTACACGTGCGCta             >  1:202715/1‑148 (MQ=255)
                                                                                                                                   cGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCtat            <  2:320964/148‑1 (MQ=255)
                                                                                                                                      tGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCt                                                                 <  6:292624/92‑1 (MQ=255)
                                                                                                                                      tGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCt                                                                 >  5:292624/1‑92 (MQ=255)
                                                                                                                                       gCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAg                                                             >  3:130035/1‑95 (MQ=255)
                                                                                                                                       gCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAg                                                             <  4:130035/95‑1 (MQ=255)
                                                                                                                                            aaaCGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTATATGGTcggc   >  5:190724/1‑148 (MQ=255)
                                                                                                                                              aCGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTATATGGTcggcg  >  6:91169/1‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
GAACGCTGCGAACTGGGCGAACAGCGTGAAGCGTTGTTGGCGGTGACCATTCCGGAAGAAAAAGGCAGCTTCCTCAAATTCTGCCAACTGCTTGGCGGGCGTTCGGTCACCGAGTTCAACTACCGTTTTGCCGATGCCAAAAACGCCTGCATCTTTGTCGGTGTGCGCCTGAGCCGCGGCCTCGAAGAGCGCAAAGAAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTATATGGTCGGCG  >  NC_000913/3956315‑3956603

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: