Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 121,230 T→C 6.3% M108V (ATG→GTG)  aroP ← aromatic amino acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913121,2300TC 6.3% 206.4 / 4.2 63M108V (ATG→GTG) aroParomatic amino acid transporter
Reads supporting (aligned to +/- strand):  ref base T (27/32);  new base C (2/2);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.08e-01

AGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAATCACCACAAAGAATACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAA  >  NC_000913/121089‑121373
                                                                                                                                             |                                                                                                                                               
aGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAATCACCACAAAGAATACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACt                                                                                                                                          >  2:244684/1‑149 (MQ=255)
 gAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAATCACCACAAAGAATACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTa                                                                                                                                         >  1:145160/1‑149 (MQ=255)
 gAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAATCACCACAAAGAATACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTa                                                                                                                                         >  4:576921/1‑149 (MQ=255)
   aCTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAATCACCACAAAGAATACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTaaa                                                                                                                                       <  2:616680/149‑1 (MQ=255)
    cTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAATCACCACAAAGAATACGGCGGCAGAAACCCATGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTaaaa                                                                                                                                      <  3:273763/149‑1 (MQ=255)
    cTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAATCACCACAAAGAATACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTaaaa                                                                                                                                      <  4:71075/149‑1 (MQ=255)
     tCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAATCACCACAAAGAATACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAAc                                                                                                                                     >  3:29760/1‑149 (MQ=255)
     tCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAATCACCACAAAGAATACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAAc                                                                                                                                     <  2:501257/149‑1 (MQ=255)
        aTCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAATCACCACAAAGAATACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTa                                                                                                                                  <  2:617882/149‑1 (MQ=255)
        aTCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAATCACCACAAAGAATACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTa                                                                                                                                  <  4:373987/149‑1 (MQ=255)
        aTCTCGACAAACACTTTAACGTTGGTCAGGTTGATGTCGTTAATCACCACAAAGAATACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCGCGGCAGTCAGCTCAGCCATGGCAACTAAAACGTa                                                                                                                                  <  3:412686/149‑1 (MQ=255)
         tctcGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAATCACCACAAAGAATACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTAc                                                                                                                                 <  4:37105/149‑1 (MQ=255)
         tctcGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAATCACCACAAAGAATACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTAc                                                                                                                                 >  2:50435/1‑149 (MQ=255)
                  acacTTTAACGTTGGTCAGGTTGATGGCGTTAATCACCACAAAGAATACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCag                                                                                                                        <  4:14595/149‑1 (MQ=255)
                  acacTTTAACGTTGGTCAGGTTGATGGCGTTAATCACCACAAAGAATACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCag                                                                                                                        <  3:515696/149‑1 (MQ=255)
                    acTTTAACGTTGGTCAGGTTGATGGCGTTAATCACCACAAAGAATACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCagta                                                                                                                      >  2:80488/1‑149 (MQ=255)
                    acTTTAACGTTGGTCAGGTTGATGGCGTTAATCACCACAAAGAATACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCagta                                                                                                                      >  3:337661/1‑149 (MQ=255)
                    acTTTAACGTTGGTCAGGTTGATGGCGTTAATCACCACAAAGAATACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCagta                                                                                                                      >  1:145902/1‑149 (MQ=255)
                     cTTTAACGTTGGTCAGGTTGATGGCGTTAATCACCACAAAGAATACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCagtag                                                                                                                     <  1:489124/149‑1 (MQ=255)
                      tttAACGTTGGTCAGGTTGATGGCGTTAATCACCACAAAGAATACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCagtagt                                                                                                                    >  4:460379/1‑149 (MQ=255)
                      tttAACGTTGGTCAGGTTGATGGCGTTAATCACCACAAAGAATACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCagtagt                                                                                                                    <  3:260635/149‑1 (MQ=255)
                          aCGTTGGTCAGGTTGATGGCGTTAATCACCACAAAGAATACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTcca                                                                                                                <  4:197130/149‑1 (MQ=255)
                                     gTTGATGGCGTTAATCACCACAAAGAATACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCGGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGa                                                                                                     >  3:492790/1‑149 (MQ=255)
                                          tGGCGTTAATCACCACAAAGAATACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCg                                                                                                <  2:145902/149‑1 (MQ=255)
                                                     accacAAAGAATACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGcc                                                                                     <  4:139232/149‑1 (MQ=255)
                                                     accacAAAGAATACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGcc                                                                                     <  1:50435/149‑1 (MQ=255)
                                                      ccacAAAGAATACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGccc                                                                                    >  1:332240/1‑149 (MQ=255)
                                                       cacaAAGAATACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTAccc                                                                                                                          >  4:354208/1‑110 (MQ=255)
                                                       cacaAAGAATACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTAccc                                                                                                                          <  3:354208/110‑1 (MQ=255)
                                                       cacaAAGAATACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGcccc                                                                                   >  3:114861/1‑149 (MQ=255)
                                                       cacaAAGAATACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGcccc                                                                                   >  4:161931/1‑149 (MQ=255)
                                                        acaAAGAATACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCa                                                                                  >  4:492613/1‑149 (MQ=255)
                                                        acaAAGAATACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCa                                                                                  >  3:227960/1‑149 (MQ=255)
                                                         caAAGAATACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGcccc                                                                                   <  1:610659/147‑1 (MQ=255)
                                                         caAAGAATACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGcccc                                                                                   >  2:610659/1‑147 (MQ=255)
                                                          aaaGAATACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGt                                                                                >  4:528904/1‑149 (MQ=255)
                                                          aaaGAATACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGt                                                                                >  4:574299/1‑149 (MQ=255)
                                                             gAATACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTAtt                                                                             >  2:308479/1‑149 (MQ=255)
                                                                tACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCACGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTtat                                                                          >  2:149681/1‑149 (MQ=255)
                                                                tACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCACGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTtat                                                                          >  1:164421/1‑149 (MQ=255)
                                                                  cggcggCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATaa                                                                        <  3:18869/149‑1 (MQ=255)
                                                                                 aGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAg                                                         <  2:559210/149‑1 (MQ=255)
                                                                                    tGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAg                                                      <  4:227960/149‑1 (MQ=255)
                                                                                    tGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAg                                                      >  4:30719/1‑149 (MQ=255)
                                                                                             ccGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACa                                             <  2:332240/149‑1 (MQ=255)
                                                                                                      aGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCg                                    >  2:200105/1‑149 (MQ=255)
                                                                                                      aGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCACGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCg                                    <  2:164421/149‑1 (MQ=255)
                                                                                                      aGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCACGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCg                                    <  1:149681/149‑1 (MQ=255)
                                                                                                               tGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCAtt                           <  4:29760/149‑1 (MQ=255)
                                                                                                                  ctttACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAg                                                         <  3:462857/115‑1 (MQ=255)
                                                                                                                  ctttACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAg                                                         >  4:462857/2‑116 (MQ=255)
                                                                                                                      aCCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCa                    <  3:574299/149‑1 (MQ=255)
                                                                                                                        ccACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGc                  >  4:58515/1‑149 (MQ=255)
                                                                                                                         cACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCt                 >  2:547384/1‑149 (MQ=255)
                                                                                                                            ggCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGAc              >  4:291820/1‑149 (MQ=255)
                                                                                                                                 tcagctcagcCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATg         <  1:244684/149‑1 (MQ=255)
                                                                                                                                  cagctcagcCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGa        <  3:12143/149‑1 (MQ=255)
                                                                                                                                   agctcagcCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGAt       <  3:492613/149‑1 (MQ=255)
                                                                                                                                    gctcagcCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATc      <  1:80488/149‑1 (MQ=255)
                                                                                                                                    gctcagcCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATc      <  2:145160/149‑1 (MQ=255)
                                                                                                                                        agcCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGaa  <  4:337661/149‑1 (MQ=255)
                                                                                                                                            aTGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGa        <  4:420038/139‑1 (MQ=255)
                                                                                                                                            aTGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGa        >  3:420038/1‑139 (MQ=255)
                                                                                                                                             |                                                                                                                                               
AGAACTCCATCTCGCCAAACACTTTAACGTTGGTCAGGTTGATGGCGTTAATCACCACAAAGAATACGGCGGCAGAAACCCAGGTGGGGATTTCCGGATACCAGAACTGAATGTATTTACCCACGGCAGTCAGCTCAGCCATGGCAACTAAAACGTACAGTACCCAGTAGTTCCAGCCAGAGGCGAAACCGGCAAAACTGCCCCAGTATTTATAAGCAAAGTGGCTAAAGGAGCCTGCGACAGGTTCTTCGACCACCATTTCACCCAGCTGACGCATGATCAGAA  >  NC_000913/121089‑121373

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