Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,852,061 T→C 6.4% M242V (ATG→GTG)  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,0610TC 6.4% 189.9 / 3.8 62M242V (ATG→GTG) ilvBacetolactate synthase 2 large subunit
Reads supporting (aligned to +/- strand):  ref base T (29/29);  new base C (2/2);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

GGTCATCAAAACGCGCACCGAGCACTATCAACAAATCCGCCTCCTGCAAAATATAGTTGGTGCTGCGCACGCCGTGCATCCCCAGCATACCCAGCGACAACGGATGCGCTTTTGGCAACATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATCACACCGCCGCCCAGATAAAGCACCGGGCGTTTGGCAGCGTTAATCATCGCCGCTGCGTCACGAATGCTTTCTTCGCTAAAGGCGGGGGC  >  NC_000913/3851920‑3852208
                                                                                                                                             |                                                                                                                                                   
ggTCATCAAAACGCGCACCGAGCACTATCAACAAATCCGCCTCCTGCAAAATATAGTTGGTGCTGCGCACGCCGTGCATCCCCAGCATACCCAGCGACAACGGATGCGCTTTTGGCAACATGCCCAGCGCCATTAAAGTCATggtggt                                                                                                                                               <  7:174923/148‑1 (MQ=255)
     taaaaaCTCTCACCGATCACTATCAACAAATCCGCCTCCTTAAAAATTTAGTTGGTGCTGCGCACGCCGTGCATCCCCAGCATACCCAGCGACAACGGATGCGCTTTTGGGAAGATGCGCAGCGCCATTAAAGTCATGGTGGTAGGCa                                                                                                                                          <  4:347679/146‑1 (MQ=255)
      caAAACGCGCACCGAGCACTATCAACAAATCCGCCTCCTGCAAAATATAGTTGGTGCTGCGCACGCCGTGCATCCCCAGCATACCCAGCGACAACGGATGCGCTTTTGGCAACATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAg                                                                                                                                         >  3:288771/1‑148 (MQ=255)
       aaaaCGCGCACCGAGCACTATCAACAAATCCGCCTCCTGCAAAATATAGTTGGTGCTGCGCACGCCGTGCATCCCCCGAATCCCCAGCGACAACGGATGCGCTTTTGGCAACATGCCCAGCGCCATTAAAGTCATGTTGGTAGGCAGt                                                                                                                                        >  6:204413/1‑148 (MQ=255)
                   gAGCACTATCAACAAATCCGCCTCCTGCAAAATATAGTTGGTGCTGCGCACGCCGTGCATCCCCAGCATACCCAGCGACAACGGATGCGCTTTTGGCAACATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTcc                                                                                                                            >  6:15387/1‑148 (MQ=255)
                      cACTATCAACAAATCCGCCTCCTGCAAAATATAGTTGGTGCTGCGCACGCCGTGCATCCCCAGCATACCCAGCGACAACGGATGCGCTTTTGGCAACATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGcc                                                                                                                         <  2:334771/148‑1 (MQ=255)
                          aTCAACAAATCCGCCTCCTGCAAAATATAGTTGGTGCTGCGCACGCCGTGCATCCCCAGCATACCCAGCGACAACGGATGCGCTTTTGGCAACATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGtt                                                                                                                     >  3:51377/1‑148 (MQ=255)
                           tCAACAAATCCGCCGACTGCAAAATATAGTTGGGGCTGCGCACGCCGTGCATCCCCAGCATACCCACCGACAACGGATGCGCTTTTGGCAACATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGtt                                                                                                                     <  5:38663/147‑1 (MQ=255)
                            caacaaATCCGCCTCCTGCAAAATATAGTTGGTGCTGCGCACGCCGTGCATCCCCAGCATACCCAGCGACAACGGATGCGCTTTTGGCAACATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCa                                                                                                                   >  1:400543/1‑148 (MQ=255)
                                    ccGCCTCCTGCAAAATATAGTTGGTGCTGCGCACGCCGTGCATCCCCAGCATACCCAGCGACAACGGATGCGCTTTTGGCAACATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCg                                                                                                           >  4:335802/1‑148 (MQ=255)
                                       cctcctGCAAAATATAGTTGGTGCTGCGCACGCCGTGCATCCCCAGCATACCCAGCGACAACGGATGCGCTTTTGGCAACATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGc                                                                                                        >  1:275517/1‑148 (MQ=255)
                                           ctGCAAAATATAGTTGGTGCTGCGCACGCCGTGCATCCCCAGCATACCCAGCGACAACGGATGCGCTTTTGGCAACATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCggg                                                                                                     <  5:89277/147‑1 (MQ=255)
                                             gCAAAATATAGTTGGTGCTGCGCACGCCGTGCATCCCCAGCATACCCAGCGACAACGGATGCGCTTTTGGCAACATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGgcg                                                                                                   >  8:252369/1‑147 (MQ=255)
                                              cAAAATATAGTTGGTGCTGCGCACGCCGTGCATCCCCAGCATACCCAGCGACAACGGATGCGCTTTTGGCAACATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGgcgc                                                                                                  <  6:272041/147‑1 (MQ=255)
                                               aaaaTATAGTTGGTGCTGCGCACGCCGTGCATCCCCAGCATACCCAGCGACAACGGATGCGCTTTTGGCAACATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCAt                                                                                                <  4:288771/148‑1 (MQ=255)
                                               aaaaTATAGTTGGTGCTGCGCACGCCGTGCATCCCCAGCATACCCAGCGACAACGGATGCGCTTTTGGCAACATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCAt                                                                                                <  2:275517/148‑1 (MQ=255)
                                                  atatAGTTGGTGCTGCGCACGCCGTGCATCCCCAGCATACCCAGCGACAACGGATGCGCTTTTGGCAACATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGa                                                                                             <  2:85043/148‑1 (MQ=255)
                                                   tataGTTGGTGCTGCGCACGCCGTGCATCCCCAGCATACCCAGCGACAACGGATGCGCTTTTGGCAACATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGAt                                                                                            <  5:204978/148‑1 (MQ=255)
                                                     taGTTGGTGCTGCGCACGCCGTGCATCCCCAGCATACCCAGCGACAACGGATGCGCTTTTGGCAACATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATc                                                                                           <  5:15387/147‑1 (MQ=255)
                                                     taGTTGGTGCTGCGCACGCCGTGCATCCCCAGCATACCCAGCGACAACGGATGCGCTTTTGGCAACATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATc                                                                                           <  5:66351/147‑1 (MQ=255)
                                                     taGTTGGTGCTGCGCACGCCGTGCATCCCCAGCATACCCAGCGACAACGGATGCGCTTTTGGCAACATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATc                                                                                           <  2:400543/147‑1 (MQ=255)
                                                                          tGCATCCCCAGCATACCCAGCGACAACGGATGCGCTTTTGGCAACATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATCACACCGCCGCCCAGATAAAGCa                                                                     <  7:99642/148‑1 (MQ=255)
                                                                            cATCCCCAGCATACCCAGCGACAACGGATGCGCTTTTGGCAACATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATCACACCGCCGCCCAGATAAAGCAcc                                                                   >  7:136292/1‑148 (MQ=255)
                                                                              tCCCCAGCATACCCAGCGACAACGGATGCGCTTTTGGCAACATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCAACCGTGCGGGCGCATTGATCACACCGCCGCCCAGATAAAGCACCgg                                                                 >  8:366818/1‑148 (MQ=255)
                                                                               ccccAGCATACCCAGCGACAACGGATGCGCTTTTGGCAACATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATCACACCGCCGCCCAGATAAAGCACCggg                                                                >  2:388777/1‑148 (MQ=255)
                                                                                cccAGCATACCCAGCGACAACGGATGCGCTTTTGGCAACATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATCACACCGCCGCCCAGATAAAGCACCGGGc                                                               >  1:177728/1‑148 (MQ=255)
                                                                                cccAGCATACCCAGCGACAACGGATGCGCTTTTGGCAACATGCCCAGCGCCATTAAAGTCACGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATCACACCGCCGCCCAGATAAAGCACCGGGc                                                               >  1:155128/1‑148 (MQ=255)
                                                                                  cAGCATACCCAGCGACAACGGATGCGCTTTTGGCAACATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATCACACCGCCGCCCAGATAAAGCACCGGGCGt                                                             <  1:89638/148‑1 (MQ=255)
                                                                                        acacAGCGACAACGGATGCGCTTTTGGCAACATGCCCAGAGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATCACACCGCCGCCCAGATAAAGAACCGGGCGTTTgg                                                         >  3:219666/4‑146 (MQ=255)
                                                                                        aCCCAGCGACAACGGATGCGCTTTTGGCAACATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATCACACCGCCGCCCAGATAAAGCa                                                                     >  5:380572/1‑134 (MQ=255)
                                                                                        aCCCAGCGACAACGGATGCGCTTTTGGCAACATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATCACACCGCCGCCCAGATAAAGCa                                                                     <  6:380572/134‑1 (MQ=255)
                                                                                           cAGCGACAACGGATGCGCTTTTGGCAACATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATCACACCGCCGCCCAGATAAAGCACCGGGCGTTTGGCAGCg                                                    >  6:65916/1‑148 (MQ=255)
                                                                                              cGACAACGGATGCGCTTTTGGCAACATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATcaca                                                                                        <  5:260401/109‑1 (MQ=255)
                                                                                              cGACAACGGATGCGCTTTTGGCAACATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATcaca                                                                                        >  6:260401/1‑109 (MQ=255)
                                                                                                      gATGCGCTTTTGGCAACATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATCACACCGCCGCCCAGATAAAGCACCGGGCGTTTGGCAGCGTTAATCATcgc                                         >  4:54089/1‑148 (MQ=255)
                                                                                                      gATGCGCTTTTGGCAACATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATCACACCGCCGCCCAGATAAAGCACCGGGCGTTTGGCAGCGTTAATCATcgc                                         <  2:234278/148‑1 (MQ=255)
                                                                                                        tGCGCTTTTGGCAACATGCCCAGCGCCAGTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATCACACCGCCGCCCAGATACAGCACCGGGCGTTTGGCAGCGTTAATCATcgccg                                       <  6:200186/148‑1 (MQ=255)
                                                                                                           gcTTTTGGCAACATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATCACACCGCCGCCCAGATAAAGCACCGGGCGTTTGGCAGCGTTAATCATCGCCGCTg                                    <  5:144587/148‑1 (MQ=255)
                                                                                                            cTTTTGGCAACATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATCACACCGCCGCCCAGATAAAGCACCGGGCGtt                                                            >  4:313383/1‑123 (MQ=255)
                                                                                                            cTTTTGGCAACATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATCACACCGCCGCCCAGATAAAGCACCGGGCGtt                                                            <  3:313383/123‑1 (MQ=255)
                                                                                                             ttttGGCAACATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATCACACCGCCGCCCAGATAAAGCACCGGGCGTTTGGCAGCGTTAATCATCGCCGCTGCg                                  >  8:388770/1‑148 (MQ=255)
                                                                                                                   cAACATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATCACACCGCCGCCCAGATAAAGCACCGGGCGTTTGGCAGCGTTAATCATCGCCGCTGCGTCACGa                            <  1:388777/148‑1 (MQ=255)
                                                                                                                   cAACATGCCC‑GCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGGTCACACCGCCGCCCAGATAAAGCACCGGGCGTTTGGCAGCGTTAATCATCGCCGCTGCGTCACGaa                           >  8:203431/1‑148 (MQ=255)
                                                                                                                     aCATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATCACACCGCCGCCCAGATAAAGCACCGCGCGTTTGGCAGCGTTAATCATCGCCGCTGCGTCACGAAt                          >  3:190755/1‑148 (MQ=255)
                                                                                                                      cATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATCACACCGCCGCCCAGATaaa                                                                        <  7:68927/101‑1 (MQ=255)
                                                                                                                      cATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATCACACCGCCGCCCAGATaaa                                                                        >  8:68927/1‑101 (MQ=255)
                                                                                                                      cATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATCACACCGCCGCCCAGATAAAGCACCGGGCGTTTGGCAGCGTTAATCATCGCCGCTGCGTCACGAATg                         >  4:131229/1‑148 (MQ=255)
                                                                                                                      cATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATCACACCGCCGCCCAGATAAAGCACCGGGCGTTTGGCAGCGTTAATCATCGCCGCTGCGTCACGAATg                         >  7:181776/1‑148 (MQ=255)
                                                                                                                       aTGCCCAGCGCCATTAAAGTCACGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATCACACCGCCGCCCAGATAAAGCACCGGGCGTTTGGCAGCGTTAATCATCGCCGCTGCGTCACGAATGc                        <  2:155128/148‑1 (MQ=255)
                                                                                                                          cccAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATCACACCGCCGCCCAGATAAAGCACCGGGCGTTTGGCAGCGTTAATCATCGCCGCTGCGTCACGAATGCtt                      <  7:388770/147‑1 (MQ=255)
                                                                                                                          cccAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATCACACCGCCGCCCAGATAAAGCACCGGGCGTTTGGCAGCGTTAATCATCGCCGCTGCGTCACGAATGCtt                      <  3:131229/147‑1 (MQ=255)
                                                                                                                            cAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATCACACCGCCGCCCAGATAAAGCACCGGGCGTTTGGCAGCGTTAATCATCGCCGCTGCGTCACGAATGCTttct                   >  5:368572/1‑148 (MQ=255)
                                                                                                                             aGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATCACACCGCCGCCCAGATAAAGCACCGGGCGTTTGGCAGCGTTAATCATCGCCGCTGCGTCACGAATGCTttctt                  >  4:375717/1‑148 (MQ=255)
                                                                                                                               cgcCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATCACACCGCCGCCCAGATAAAGCACCGGGCGTTTGGCAGCGGTAATCATCGCCGCTGCGTCACGAATGCTTTCTTCg                <  3:199874/148‑1 (MQ=255)
                                                                                                                                   aTTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCAGCCAGTTCACGCACCCGTGCGGGCGCATTGATCACACCGCCGCCCAGATAAAGCACCGGGCGTTTGGCAGCGTTAATCATCGCCGCTCCGTCACGAATGCTTTCTTCGCTa             >  4:3553/1‑147 (MQ=255)
                                                                                                                                     tAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATCACACCGCCGCCCAGATAAAGCACCGGGCGTTTGGCAGCGTTAATCATCGCCGCTGCGTCACGAATGc                        <  4:355086/134‑1 (MQ=255)
                                                                                                                                     tAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATCACACCGCCGCCCAGATAAAGCACCGGGCGTTTGGCAGCGTTAATCATCGCCGCTGCGTCACGAATGc                        >  3:355086/1‑134 (MQ=255)
                                                                                                                                     tAAAGTCACGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCa                                                                                                 >  7:252292/1‑61 (MQ=255)
                                                                                                                                     tAAAGTCACGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCa                                                                                                 <  8:252292/61‑1 (MQ=255)
                                                                                                                                      aaaGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATCACACCGCCGCCCAGATaa                                                                         >  4:297379/1‑84 (MQ=255)
                                                                                                                                      aaaGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATCACACCGCCGCCCAGATaa                                                                         <  3:297379/84‑1 (MQ=255)
                                                                                                                                            aTGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATCACACCGCCGCCCAGATAAAGCACCGGGCGTTTGGCAGCGTTAATCATCGCCGCTGCGTCACGAATGCTTTCTTCGCTAAAGGCgggg    <  3:335802/147‑1 (MQ=255)
                                                                                                                                             tggtggTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATCACACCGCCGCCCAGATAAAGCACCGGGCGTTTGGCAGCGTTAATCATCGCCGCTGCGTCACGAATGCTTTCTTCGCTAAAGGCGGGggc  >  6:321683/1‑148 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
GGTCATCAAAACGCGCACCGAGCACTATCAACAAATCCGCCTCCTGCAAAATATAGTTGGTGCTGCGCACGCCGTGCATCCCCAGCATACCCAGCGACAACGGATGCGCTTTTGGCAACATGCCCAGCGCCATTAAAGTCATGGTGGTAGGCAGTTGCGCTTTCTCCGCCAGTTCACGCACCCGTGCGGGCGCATTGATCACACCGCCGCCCAGATAAAGCACCGGGCGTTTGGCAGCGTTAATCATCGCCGCTGCGTCACGAATGCTTTCTTCGCTAAAGGCGGGGGC  >  NC_000913/3851920‑3852208

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