Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,228,183 C→T 5.1% S35S (AGC→AGT ycgM → putative isomerase/hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,228,1830CT 5.1% 315.1 / 4.4 99S35S (AGC→AGTycgMputative isomerase/hydrolase
Reads supporting (aligned to +/- strand):  ref base C (45/48);  new base T (3/2);  total (49/50)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

GTCGGCGCGGGTTTTCATATCCACGATAAGGTGAGGGGAACGTTATGTATCAACATCACAACTGGCAAGGTGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTCGGTTCAGTTCATCATGAAGTCGAACTGGCGGTGTTGATTGGCGCGACACTGCGTCAGG  >  NC_000913/1228035‑1228328
                                                                                                                                                    |                                                                                                                                                 
gTCGGCGCGGGTTTTCATATCCACGATAAGGTGAGGGGAACGTTATGTATCAACATCACAACTGGCAAGGTGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAgc                                                                                                                                                   >  6:403081/1‑149 (MQ=255)
     cgcgGGTTTTCATATCCACGATAAGGTGAGGGGAACGTTATGTATCAACATCACAACTGGCAAGGTGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGt                                                                                                                                              <  7:67315/149‑1 (MQ=255)
     cgcgGGTTTTCATATCCACGATAAGGTGAGGGGAACGTTATGTATCAACATCACAACTGGCAAGGTGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGt                                                                                                                                              <  3:132450/149‑1 (MQ=255)
        gggTTTTCATATCCACGATAAGGTGAGGGGAACGTTATGTATCAACATCACAACTGGCAAGGTGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGcc                                                                                                                                           <  6:62057/149‑1 (MQ=255)
                atatCCACGATAAGGTGAGGGGAACGTTATGTATCAACATCACAACTGGCAAGGTGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGc                                                                                                                                   >  7:199306/1‑149 (MQ=255)
                 tatCCACGATAAGGTGAGGGGAACGTTATGTATCAACATCACAACTGGCAAGGTGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGc                                                                                                                                   >  6:161111/1‑148 (MQ=255)
                 tatCCACGATAAGGTGAGGGGAACGTTATGTATCAACATCACAACTGGCAAGGTGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGACAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGc                                                                                                                                   >  6:456595/1‑148 (MQ=255)
                  atCCACGATAAGGTGAGGGGAACGTTATGTATCAACATCACAACTGGCAAGGTGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGc                                                                                                                                   >  8:125290/1‑147 (MQ=255)
                   tCCACGATAAGGTGAGGGGAACGTTATGTATCAACATCACAACTGGCAAGGTGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGAGCAGTGCCCGAAGAGCCAg                                                                                                                                >  3:52690/1‑149 (MQ=255)
                    ccACGATAAGGTGAGGGGAACGTTATGTATCAACATCACAACTGGCAAGGTGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGt                                                                                                                               >  1:109643/1‑149 (MQ=255)
                      aCGATAAGGTGAGGGGAACGTTATGTATCAACATCACAACTGGCAAGGTGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGc                                                                                                                             >  3:126829/1‑149 (MQ=255)
                        gATAAGGTGAGGGGAACGTTATGTATCAACATCACAACTGGCAAGGTGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTg                                                                                                                           >  2:268819/1‑149 (MQ=255)
                         aTAAGGTGAGGGGAACGTTATGTATCAACATCACAACTGGCAAGGTGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGt                                                                                                                          <  5:110013/149‑1 (MQ=255)
                            aGGTGAGGGGAACGTTATGTATCAACATCACAACTGGCAAGGTGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTtta                                                                                                                       <  5:161111/149‑1 (MQ=255)
                              gTGAGGGGAACGTTATGTATCAACATCACAACTGGCAAGGTGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTttat                                                                                                                      >  5:377376/1‑148 (MQ=255)
                               tGAGGGGAACGTTATGTATCAACATCACAACTGGCAAGGTGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGc                                                                                                                                   >  6:370026/1‑134 (MQ=255)
                               tGAGGGGAACGTTATGTATCAACATCACAACTGGCAAGGTGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGc                                                                                                                                   <  5:370026/134‑1 (MQ=255)
                                 aGGGGAACGTTATGTATCAACATCACAACTGGCAAGGTGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTaaa                                                                                                                  <  8:199306/149‑1 (MQ=255)
                                 aGGGGAACGTTATGTATCAACATCACAACTGGCAAGGTGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTaaa                                                                                                                  >  2:100924/1‑149 (MQ=255)
                                     gAACGGTATGTATCAACATCACAACTGGCAAGGTGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAg                                                                                                              <  8:223206/149‑1 (MQ=255)
                                          ttATGTATCAACATCACAACTGGCAAGGTGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGaa                                                                                                                                       >  1:466015/1‑119 (MQ=255)
                                          ttATGTATCAACATCACAACTGGCAAGGTGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGaa                                                                                                                                       <  2:466015/119‑1 (MQ=255)
                                             tGTATCAACATCACAACTGGCAAGGTGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCa                                                                                                      >  3:47364/1‑149 (MQ=255)
                                                aTCAACATCACAACTGGCAAGGTGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACtg                                                                                                   >  6:388914/1‑149 (MQ=255)
                                                   aaCATCACAAATGGCAAGGTGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCATTGCCCGAATATCCAGTGCTGTTGATTAAACCAGAAACGGCAATGTGc                                                                                                >  8:56421/1‑149 (MQ=255)
                                                       tCACAACTGGCAAGGTGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATc                                                                                            <  7:125290/149‑1 (MQ=255)
                                                         acaACTGGCAAGGTGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTg                                                                                          >  6:129594/1‑149 (MQ=255)
                                                          caACTGGCAAGGTGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGc                                                                                         <  5:388914/149‑1 (MQ=255)
                                                           aaCTGGCAAGGTGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGt                                                                                                                          >  8:418594/1‑115 (MQ=255)
                                                           aaCTGGCAAGGTGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGt                                                                                                                          <  7:418594/115‑1 (MQ=255)
                                                           aaCTGGCAAGGTGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCg                                                                                        >  7:403104/1‑149 (MQ=255)
                                                           aaCTGGCAAGGTGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCg                                                                                        >  1:370804/1‑149 (MQ=255)
                                                             cTGGCAAGGTGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACg                                                                                                         <  8:124650/130‑1 (MQ=255)
                                                             cTGGCAAGGTGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACg                                                                                                         >  7:124650/1‑130 (MQ=255)
                                                                gCAAGGTGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCa                                                                                     <  2:109643/147‑1 (MQ=255)
                                                                   aGGTGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCAt                                                                                <  5:129594/149‑1 (MQ=255)
                                                                     gTGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTg                                                                              <  4:126829/149‑1 (MQ=255)
                                                                      tGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTgg                                                                             >  6:273404/1‑149 (MQ=255)
                                                                           tgctgGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATc                                                                        >  7:246848/1‑149 (MQ=255)
                                                                             ctgGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATcc                                                                       <  2:175852/148‑1 (MQ=255)
                                                                             ctgGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATcc                                                                       >  7:119166/1‑148 (MQ=255)
                                                                             ctgGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATcc                                                                       >  1:175852/1‑148 (MQ=255)
                                                                             ctgGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATc                                                                        <  1:100924/147‑1 (MQ=255)
                                                                              tgGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGc                                                                                                                                   <  2:93375/87‑1 (MQ=255)
                                                                              tgGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGc                                                                                                                                   >  1:93375/1‑87 (MQ=255)
                                                                               ggATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCAt                                                                    >  6:379866/1‑149 (MQ=255)
                                                                                gATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATc                                                                   >  8:50183/1‑149 (MQ=255)
                                                                                gATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATc                                                                   >  7:219645/1‑149 (MQ=255)
                                                                                 aTTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGGAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATcc                                                                  <  4:3880/149‑1 (MQ=255)
                                                                                    aTCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGTGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGAt                                                               >  5:181341/1‑149 (MQ=255)
                                                                                    aTCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGTGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGAt                                                               >  5:187901/1‑149 (MQ=255)
                                                                                     tCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGAtt                                                              >  2:78706/1‑149 (MQ=255)
                                                                                       cGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGAttt                                                             <  5:273404/148‑1 (MQ=255)
                                                                                       cGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTc                                                            >  1:307039/1‑149 (MQ=255)
                                                                                        ggTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTCg                                                           <  1:268819/149‑1 (MQ=255)
                                                                                        ggTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTCg                                                           <  5:403081/149‑1 (MQ=255)
                                                                                         gTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTCgg                                                          >  2:284/1‑149 (MQ=255)
                                                                                             gTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTCGgtt                                                        >  8:188928/1‑147 (MQ=255)
                                                                                              tAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTATAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTCGgttca                                                      <  4:47364/148‑1 (MQ=255)
                                                                                              tAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGTGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTCGgttcag                                                     <  6:187901/149‑1 (MQ=255)
                                                                                              tAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTg                                                                              >  2:66485/1‑124 (MQ=255)
                                                                                              tAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTg                                                                              <  1:66485/124‑1 (MQ=255)
                                                                                              tAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTCGgttcag                                                     >  2:194097/1‑149 (MQ=255)
                                                                                               aaaGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGTGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTCGgttcag                                                     <  6:181341/148‑1 (MQ=255)
                                                                                                     gtCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTCGGTTCAGTtcatc                                               <  2:370804/148‑1 (MQ=255)
                                                                                                      tCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTCGGTTCAGTTcatcat                                             <  1:78706/149‑1 (MQ=255)
                                                                                                        tgtgTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTCGGTTCAGTtcatca                                              <  1:37151/146‑1 (MQ=255)
                                                                                                        tgtgTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTCGGTTCAGTtcatca                                              >  2:37151/1‑146 (MQ=255)
                                                                                                        tgtgTTGGCAGTAACTATGCCAAAAATATTAAAGAGATGGGCAGTGCACTACACGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGAACTGCGGCAGCCATTATCGATCCCATCCGATTTCGGTTCAGTTAATCATg                                            >  8:441834/1‑148 (MQ=255)
                                                                                                          tgtTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTCGGTTCAGTTCATCATGaa                                          <  8:219645/148‑1 (MQ=255)
                                                                                                          tgtTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTCGGTTCAGTTCATCATGaa                                          <  7:320131/148‑1 (MQ=255)
                                                                                                            ttGGCAGTAACTATTCCAAACATATTAAAGCTATGGGCATCGCAGTTCCCGAAGAGCCATTGCTGTTTATTAAACCAGAAAAGGCACTGTGCTATCTTCGGCAGCCATTGGCGATCCCATCCGATTTCGTTTCAGTTCATCATGAAGt                                        <  7:56421/148‑1 (MQ=255)
                                                                                                               gCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTCGGTTCAGTTCATCATGAAGTCGaa                                    <  2:307039/149‑1 (MQ=255)
                                                                                                                cAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGAttt                                                             <  6:402101/123‑1 (MQ=255)
                                                                                                                  gTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTCGGTTCAGTTCATCATGAAGTCGAACTg                                 >  1:334881/1‑149 (MQ=255)
                                                                                                                   tAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTCGGTTCAGTTCATCATGAAGTCGAACTgg                                >  6:185510/1‑149 (MQ=255)
                                                                                                                   tAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTCGGTTCAGTTCATCATGAAGTCGAACTgg                                <  1:284/149‑1 (MQ=255)
                                                                                                                         tGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTCGGTTCAGTTCATCATGAAGTCGAACTGGCGgtgt                          <  8:119166/149‑1 (MQ=255)
                                                                                                                           ccAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATc                                                                   <  4:146217/106‑1 (MQ=255)
                                                                                                                            cAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTCGGTTCAGTTCATCATGAAGTCGAACTGGCg                              <  7:39344/142‑1 (MQ=255)
                                                                                                                            cAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTCGGTTCAGTTCATCATGAAGTCGAACTGGCg                              >  8:39344/1‑142 (MQ=255)
                                                                                                                            cAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTCGGTTCAGTTCATCATGAAGTCGAACTGGCGGTGTTGa                       >  2:234458/1‑149 (MQ=255)
                                                                                                                            cAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTCGGTTCAGTTCATCATGAAGTCGAACTAGCGGTGTTGa                       <  5:456595/149‑1 (MQ=255)
                                                                                                                             aaaCATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAACCGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTGCGGTTCAGGTCATCATGAAGTCGAACTGGCGGTGTTgct                      >  6:7520/1‑147 (MQ=255)
                                                                                                                             aaaCATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTCGGTTCAGTTCATCATGAAGTCGAACTGGCGGTGTTGAt                      <  5:185510/149‑1 (MQ=255)
                                                                                                                             aaaCATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTCGGTTCAGTTCATCATGAAGTCGAACTGGCGGTGTTGAt                      <  7:50183/149‑1 (MQ=255)
                                                                                                                              aaCATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTCGGTTCAGTTcatcat                                             <  2:102099/125‑1 (MQ=255)
                                                                                                                              aaCATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTCGGTTCAGTTcatcat                                             >  1:102099/1‑125 (MQ=255)
                                                                                                                               aCATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGGGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTCGGTTCAGTTCATCATGAAGTCGAACTGGCGGTGTTGATTg                    <  4:52690/149‑1 (MQ=255)
                                                                                                                                   atTAAAGAGATGGGCAGCGGTGTGCGCGAAGAGCCTCTGATGTTTATTTAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTTTCGATCACATCCGATTTCGgttcagttct                                                 <  6:412678/116‑2 (MQ=255)
                                                                                                                                   atTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTCGGTTCAGTTCATCATGAAGTCGAACTGGCGGTGTTGATTGgcgc                >  5:39058/1‑149 (MQ=255)
                                                                                                                                   atTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGAAGCCATTGTCGATCCCATCCGATTTCGgttcagttca                                                 >  5:412678/1‑116 (MQ=255)
                                                                                                                                      aaaGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTCGGTTCAGTtcatc                                               <  3:72923/115‑1 (MQ=255)
                                                                                                                                      aaaGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTCGGTTCAGTtcatc                                               >  4:72923/1‑115 (MQ=255)
                                                                                                                                       aaGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTCGGTTCAGTTCATCATGAAGTCGAACTGGCGGTGTTGATTGGCGCGaca            >  2:260874/1‑149 (MQ=255)
                                                                                                                                          agaTGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTCGGTTCAGTTCATCATGAAGTCGAACTGGCGGTGTTGATTGGCGCGACACTg         >  8:47160/1‑149 (MQ=255)
                                                                                                                                           gaTGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTCGGTTCAGTTCATCATGAAGTCGAACTGGCGGTGTTGATTGGCGCGACACTGc        <  8:246848/149‑1 (MQ=255)
                                                                                                                                               ggCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTCGGTTCAGTTCATCATGAAGTCGAACTGGCGGTGTTGATTGGCGCGACACTGCGTCa    <  2:334881/149‑1 (MQ=255)
                                                                                                                                                gCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTCGGTTCAGTTCATCATGAAGTCGAACTGGCGGTGTTGATTGGCGCGACACTGCGTCAg   <  7:188928/149‑1 (MQ=255)
                                                                                                                                                 cAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTCGGTTCAGTTCATCATGAAGTCGAACTGGCGGTGTTGATTGGCGCGACACTGCGTCAgg  <  1:194097/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GTCGGCGCGGGTTTTCATATCCACGATAAGGTGAGGGGAACGTTATGTATCAACATCACAACTGGCAAGGTGCGCTGCTGGATTATCCGGTGAGTAAAGTAGTCTGTGTTGGCAGTAACTATGCCAAACATATTAAAGAGATGGGCAGCGCAGTGCCCGAAGAGCCAGTGCTGTTTATTAAACCAGAAACGGCACTGTGCGATCTGCGGCAGCCATTGGCGATCCCATCCGATTTCGGTTCAGTTCATCATGAAGTCGAACTGGCGGTGTTGATTGGCGCGACACTGCGTCAGG  >  NC_000913/1228035‑1228328

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: