Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,801,909 G→A 6.2% T221M (ACG→ATG)  thrS ← threonyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,801,9090GA 6.2% 205.2 / 2.5 64T221M (ACG→ATG) thrSthreonyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (35/25);  new base A (2/2);  total (37/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

GTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGAAATGATGGCAGAAACGCATGTTCGGTACGTGCGGACCGCGGCACATATCGACATATTCTTCATGGAAGTACAG  >  NC_000913/1801764‑1802057
                                                                                                                                                 |                                                                                                                                                    
gtcgtTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGccgtac                                                                                                                                                   <  2:276972/149‑1 (MQ=255)
   gtTGTGCCAGAATACCCTACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGccgtaccgt                                                                                                                                                <  1:115629/149‑1 (MQ=255)
   gtTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGccgtaccgt                                                                                                                                                <  5:11416/149‑1 (MQ=255)
            gAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGtt                                                                                                                                       <  1:71706/149‑1 (MQ=255)
                   ataccCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGTGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACAt                                                                                                                                >  3:128846/1‑149 (MQ=255)
                   ataccCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACAt                                                                                                                                >  2:74630/1‑149 (MQ=255)
                     accCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACAttt                                                                                                                              >  8:206399/1‑149 (MQ=255)
                       ccGGCGCTTCTTCCTGCATATGGTACAGGTCGCGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCAAGCAGTACCGTAAATACGTTTCAACATTttg                                                                                                                            >  7:262034/1‑149 (MQ=255)
                        cGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACAttttat                                                                                                                           <  3:316443/149‑3 (MQ=255)
                              ttcttcCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGc                                                                                                                     >  6:267318/1‑149 (MQ=255)
                                cttcCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCt                                                                                                                    <  5:254489/148‑1 (MQ=255)
                                    cTGCATATGGTACAGGTCGAACTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTcgc                                                                                                               >  4:196411/1‑149 (MQ=255)
                                      gCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTcgcca                                                                                                             >  2:45580/1‑149 (MQ=255)
                                       cATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTcgcca                                                                                                             <  3:218979/148‑1 (MQ=255)
                                           tGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTcgccacgcca                                                                                                        <  7:222321/149‑1 (MQ=255)
                                            ggTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCACCACGCCAg                                                                                                       >  3:36545/1‑149 (MQ=255)
                                             gTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAg                                                                                                       >  4:265334/1‑148 (MQ=255)
                                                   gTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGccc                                                                                                <  5:267318/149‑1 (MQ=255)
                                                       aGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACAtt                                                                                                                               <  1:137696/114‑1 (MQ=255)
                                                       aGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACAtt                                                                                                                               >  2:137696/1‑114 (MQ=255)
                                                        gCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGcc                                                                                           >  2:83528/1‑149 (MQ=255)
                                                        gCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGcc                                                                                           >  4:152218/1‑149 (MQ=255)
                                                           gTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGtt                                                                                        >  2:57179/1‑149 (MQ=255)
                                                           gTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGtt                                                                                        >  1:167828/1‑149 (MQ=255)
                                                              tACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGtttt                                                                                      <  3:152218/148‑1 (MQ=255)
                                                                ccGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTTAGTGCTGTTTTGTCTGCCCACGCCATACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCAt                                                                                   >  6:225304/1‑149 (MQ=255)
                                                                ccGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCa                                                                                    <  7:328621/148‑1 (MQ=255)
                                                                 cGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCAtt                                                                                  <  1:205902/149‑1 (MQ=255)
                                                                   aTTTTACGGTGGACGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGACTGCCCACGCCGTACCGTAAATACGATGCAACATTTTCTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTATCATAAg                                                                                >  7:217221/1‑149 (MQ=255)
                                                                         cGGTGGTCGCGTATCCCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTTAGTGCTTTTTTGTCTGCCCACGCCGTACCGAAAAAACGATGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGa                                                                           >  6:220387/1‑148 (MQ=255)
                                                                               tCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGAAatgatg                                                                    <  3:265334/149‑1 (MQ=255)
                                                                                cgcgTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGAAATGATgg                                                                   <  1:263621/149‑1 (MQ=255)
                                                                                 gcgTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCACCTGCCGTTTTCATTAGTTTGAAATGATGGc                                                                  <  6:321331/149‑1 (MQ=255)
                                                                                         cgGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGaa                                                                          >  1:207339/1‑133 (MQ=255)
                                                                                         cgGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGaa                                                                          <  2:207339/133‑1 (MQ=255)
                                                                                            cttcttCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGAAATGATgg                                                                   >  7:211588/1‑137 (MQ=255)
                                                                                            cttcttCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGAAATGATgg                                                                   <  8:211588/137‑1 (MQ=255)
                                                                                                  ccAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGAAATGATGGCAGAAACGCATGTTCg                                                   <  7:206399/147‑1 (MQ=255)
                                                                                                       cccTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGAAATGATGGCAGAAACGCATGTTCGGTACGTg                                            >  4:69427/3‑149 (MQ=255)
                                                                                                        gcTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGAAATGATGGCAGAAACGCATGTTCGGTACGTGc                                           >  5:312110/1‑149 (MQ=255)
                                                                                                            cAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGAAATGATGGCAGAAACGCATGTTCGGTACGTGCGGAc                                       >  1:337693/1‑149 (MQ=255)
                                                                                                              ggTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCACCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGAAATGATGGCAGAAACGCATGTTCGGTACGTGCGGACcg                                     <  4:36545/149‑1 (MQ=255)
                                                                                                               gTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGAAATGATGGCAGAAACGCATGTTCGGTACGTGCGGACcgc                                    >  4:179320/1‑149 (MQ=255)
                                                                                                               gTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGAAATGATGGCAGAAACGCATGTTCGGTACGTGCGGACcgc                                    >  1:296786/1‑149 (MQ=255)
                                                                                                                tAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGAAATGATGGCAGAAACGCATGTTCGGTACGTGCGGACcgcg                                   <  2:167828/149‑1 (MQ=255)
                                                                                                                     gTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGAAATGATGGCAGAAACGCATGTTCGGTACGTGCGGACCGCGGcaca                              >  6:303702/1‑149 (MQ=255)
                                                                                                                     gTTAAGTGCTTTTTTGGCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGAAATGATGGCAGAAACGCATGTTCGGTCCGCACGTACcgaacatgc                              <  1:35516/149‑8 (MQ=255)
                                                                                                                          gTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGAAATGATGGCAGAAACGCATGTTCGGTACGTGCGGACCGCGGCACata                            >  6:189389/1‑146 (MQ=255)
                                                                                                                          gTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGAAATGATGGCAGAAACGCATGTTCGGTACGTGCGGACCGCGGCACATATc                          >  7:277693/1‑148 (MQ=255)
                                                                                                                          gTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCATTCCTTTTTAATTAGTTTGAAATGATGTCAGAAACGCATGTTCTTTACTTGCGGACCGATGCACATATc                          >  2:252808/1‑148 (MQ=255)
                                                                                                                           tGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGAAatgatg                                                                    <  6:28530/105‑1 (MQ=255)
                                                                                                                           tGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGAAatgatg                                                                    >  5:28530/1‑105 (MQ=255)
                                                                                                                           tGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGAAATGATGGCAGAAACGCATGTTCGGTACGTGCGGACCGCGGCACATATCGa                        >  2:317297/1‑149 (MQ=255)
                                                                                                                              ttttttGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGAAATGATGGCAGAAACGCATGTTCGGTACGTGCGGACCGCGGCACATATCGACa                      >  8:141217/1‑148 (MQ=255)
                                                                                                                                ttttGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGAAATGATGGCAGAAACGCATGTTCGGTACGTGCGGACCGCGGCACATATCGACata                    <  1:57179/148‑1 (MQ=255)
                                                                                                                                   tGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGAAATGATGGCAGAAACGCATGTTCGGTACGTGCGGACCGCGGCACATATCGACATAttc                 <  1:83528/148‑1 (MQ=255)
                                                                                                                                   tGTCTGCCCACGCCATACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGAAATGATGGCAGAAACGCATGTTCGGTACGTGCGGACCGCGGCACATATCGACATAttc                 <  5:225304/148‑1 (MQ=255)
                                                                                                                                       tGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGAAATGATGGCAGAAACGCATGTTCGGTACGTGCGGACCGCGGCACATATCGACATATTCTTCAt            >  5:333694/1‑149 (MQ=255)
                                                                                                                                         cccACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGAAATGATGGCAGAAACGCATGTTCGGTACGTGCGGACCGCGGCACATATCGACATATTCTTCa             <  1:45580/146‑1 (MQ=255)
                                                                                                                                         cccACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGAAATGATGGCAGAAACGCATGTTCGGTACGTGCGGACCGCGGCACATATCGACATATTCTTCAt            <  6:312110/147‑1 (MQ=255)
                                                                                                                                         cccACGCCATACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGAAATGATGGCAGAAACGCATGTTCGGTACGTGCGGACCGCGGCACATATCGACATATTCTTCATg           >  6:118923/1‑148 (MQ=255)
                                                                                                                                          ccACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGAAATGATGGCAGAAACGCATGTTCGGTACGTGCGGACCGCGGCACATATCGCCATATTCTTCATGGa         >  1:75649/1‑149 (MQ=255)
                                                                                                                                          ccACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGAAATGATGGCAGAAACGCATGTTAGGTACGTGCGGACCGCGGCACATATCGACATATTCTTCATgg          >  5:87422/1‑148 (MQ=255)
                                                                                                                                           cccGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGAAATGATGGCAGAAACGCATGTTCGGTACGTGCGGACCGCGGCACATATCGACATATTCTTCATgg          >  8:184598/3‑147 (MQ=255)
                                                                                                                                           cACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGAAATGATGGCAGAAACGCATGTTCGGTACGTGCGGACCGCGGCACATATCGACATATTCTTCATgg          >  3:94588/1‑147 (MQ=255)
                                                                                                                                           cACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGAAATGATGGCAGAAACGCATGTTCGGTACGTGCGGACCGCGGCACATATCGACATATTCTTCATGGaa        <  3:196411/149‑1 (MQ=255)
                                                                                                                                           cACGCCATACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGAAATGATGGCAGAAACGCATGTTCGGTACGTGCGGACCGCGGCACATATCGACATATTCTTCATGGaa        <  5:118923/149‑1 (MQ=255)
                                                                                                                                                 gtaccgtaAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGAAATGATGGCAGAAACGCATGTTCGGTACGTGCGGACCGCGGCACATATCGACATATTCTTCATGGAAGTACAg  >  8:156799/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GTCGTTGTGCCAGAATACCATACCCGGCGCTTCTTCCTGCATATGGTACAGGTCGAGCTGTTTACCGATTTTACGGTGGTCGCGTTTCGCGGCTTCTTCCAGGCGCTGCAGGTAAGCGTTAAGTGCTTTTTTGTCTGCCCACGCCGTACCGTAAATACGTTGCAACATTTTGTTGTTGCTGTCGCCACGCCAGTAAGCCCCTGCCGTTTTCATTAGTTTGAAATGATGGCAGAAACGCATGTTCGGTACGTGCGGACCGCGGCACATATCGACATATTCTTCATGGAAGTACAG  >  NC_000913/1801764‑1802057

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: ‑