Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,850,154 C→A 6.6% Q94H (CAG→CAT)  hycA ← regulator of the transcriptional regulator FhlA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,850,1540CA 6.6% 206.8 / 3.7 60Q94H (CAG→CAT) hycAregulator of the transcriptional regulator FhlA
Reads supporting (aligned to +/- strand):  ref base C (24/32);  new base A (2/2);  total (26/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 = 6.27e-01

TTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTTCAAAAAGGACGAAACAGAGTTCTTTGTCCGGCGCGCAGCTCAT  >  NC_000913/2850009‑2850300
                                                                                                                                                 |                                                                                                                                                  
ttAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTgc                                                                                                                                                 >  2:111768/1‑149 (MQ=255)
     ttCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAAt                                                                                                                                            >  2:275159/1‑149 (MQ=255)
        aCCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCAGGCGCaactga                                                                                                                                         <  1:255105/149‑5 (MQ=255)
          cGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCt                                                                                                                                       <  2:331409/149‑1 (MQ=255)
             cacaTTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCAGGCGCaactg                                                                                                                                          <  4:312441/143‑4 (MQ=255)
             cacaTTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCAGGCGCaactg                                                                                                                                          >  3:312441/1‑140 (MQ=255)
              acaTTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTGCCAGCACTTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGttt                                                                                                                                   <  2:291362/149‑1 (MQ=255)
                aTTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGt                                                                                                                                 >  1:52842/1‑149 (MQ=255)
                 ttGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTc                                                                                                                                >  2:143285/1‑149 (MQ=255)
                  tGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTc                                                                                                                                >  4:205052/1‑148 (MQ=255)
                    ttttCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTc                                                                                                                             <  1:321034/149‑1 (MQ=255)
                             gCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCttt                                                                                                                    <  2:203435/149‑1 (MQ=255)
                                   tataAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTc                                                                                                              >  1:208432/1‑149 (MQ=255)
                                     ggatcagTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTCTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGa                                                                                                            <  2:359683/143‑1 (MQ=255)
                                       aGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACg                                                                                                          >  1:121449/1‑149 (MQ=255)
                                       aGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACg                                                                                                          >  2:58748/1‑149 (MQ=255)
                                         cTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTa                                                                                                        <  2:17227/149‑1 (MQ=255)
                                         cTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTa                                                                                                        >  2:190074/1‑149 (MQ=255)
                                          tGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTa                                                                                                        <  4:83344/148‑1 (MQ=255)
                                          tGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTa                                                                                                        <  3:205052/148‑1 (MQ=255)
                                          tGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTa                                                                                                        <  3:291121/148‑1 (MQ=255)
                                            tCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATa                                                                                                     >  2:213978/1‑149 (MQ=255)
                                               tAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAGCCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGTCTGCGCaa                                                                                                                                             <  4:139702/106‑1 (MQ=255)
                                               tAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCaa                                                                                                                                             >  3:139702/1‑106 (MQ=255)
                                               tAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTc                                                                                                  >  2:197282/1‑149 (MQ=255)
                                                aaaCGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTc                                                                                                  >  4:273285/1‑148 (MQ=255)
                                                      cAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTAACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGgtg                                                                                           >  2:118909/1‑149 (MQ=255)
                                                       aGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGgtgt                                                                                          <  2:136371/149‑1 (MQ=255)
                                                        gCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGtgtg                                                                                         <  1:116248/149‑1 (MQ=255)
                                                         cGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTgg                                                                                        <  2:252026/149‑1 (MQ=255)
                                                          gATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGc                                                                                       >  2:65053/1‑149 (MQ=255)
                                                          gATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGc                                                                                       >  2:239290/1‑149 (MQ=255)
                                                           aTTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGc                                                                                       <  3:386532/148‑1 (MQ=255)
                                                                 tCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTaaa                                                                                <  2:208432/149‑1 (MQ=255)
                                                                 tCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTaa                                                                                 <  3:273285/148‑1 (MQ=255)
                                                                  cGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAg                                                                               <  1:213978/149‑1 (MQ=255)
                                                                  cGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAg                                                                               <  1:58748/149‑1 (MQ=255)
                                                                          tgtTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCg                                                                        <  4:354860/148‑1 (MQ=255)
                                                                            ttCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTTTAACAGCCacaccatc                                                                     <  1:243858/149‑8 (MQ=255)
                                                                            ttCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCa                                                                     >  2:156196/1‑149 (MQ=255)
                                                                             tCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAg                                                                    <  1:126050/149‑1 (MQ=255)
                                                                                 gAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGACATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGt                                                                 <  3:240375/148‑1 (MQ=255)
                                                                                  aaCCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGa                                                               >  1:35998/1‑149 (MQ=255)
                                                                                    ccTGTTCGCGATCGTGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGAc                                                              >  3:199890/1‑148 (MQ=255)
                                                                                             gATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGa                                                    <  1:239290/149‑1 (MQ=255)
                                                                                                      tGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTTCa                                           <  1:65053/149‑1 (MQ=255)
                                                                                                           aCCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTTCaaaaa                                       >  3:15809/1‑148 (MQ=255)
                                                                                                                  cATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAAt                                                >  2:279427/1‑132 (MQ=255)
                                                                                                                  cATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAAt                                                <  1:279427/132‑1 (MQ=255)
                                                                                                                  cATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTTCAAAAAGGACGaa                                <  3:198901/148‑1 (MQ=255)
                                                                                                                   atcaatcaTGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTTCAAAAAGGACGaaa                               >  4:294897/1‑148 (MQ=255)
                                                                                                                   atcaatcaTGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTTCAAAAAGGACGAAAc                              >  1:301607/1‑149 (MQ=255)
                                                                                                                    tcaatcaTGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTaaa                                                                                <  1:178517/98‑1 (MQ=255)
                                                                                                                    tcaatcaTGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTaaa                                                                                >  2:178517/1‑98 (MQ=255)
                                                                                                                      aatcgTGCAATCAATGCTCAGTGGCGAATAAGGAATCCTATGTTTATCTTAGCGATCTTTTGTTTCGACTTATTAATCGATGGTGTGGCTGTTAAAGGCGACCGCCCGTGTGACCTCAATGCGAAAAGGTGCTAAAAGAACGAAACgga                           >  2:44996/1‑149 (MQ=255)
                                                                                                                       atcaTGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTTCAAAAAGGACGAAACaga                           <  3:253048/148‑1 (MQ=255)
                                                                                                                        tcaTGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAGAATGTTCAAAAAGGACGAAACAGAGt                         <  1:156196/149‑1 (MQ=255)
                                                                                                                                tCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTTCAAAAAGGACGAAACAGAGTTCTTTGTc                 <  1:190074/149‑1 (MQ=255)
                                                                                                                                    tGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGTCGTAAATGCGAAAATGTTCAAAAAGGACGAAACAGAGTTCTTTGTCCGgc             <  2:121449/149‑1 (MQ=255)
                                                                                                                                            tGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTCCAAAAAGGACGAAACAGAGTTTTTTGTCCGGCGCGCAGCt     <  1:118909/149‑1 (MQ=255)
                                                                                                                                                ccTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTTCAAAAAGGACGAAACAGAGTTCTTTGTCCGGCGCGCAGCTCAt  <  4:48521/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCATCTTTTGTTTCGACGTAATACTCGATGGTGTGGCTGTTAAAGCCATCCGCCAGGGTGACGTAAATGCGAAAATGTTCAAAAAGGACGAAACAGAGTTCTTTGTCCGGCGCGCAGCTCAT  >  NC_000913/2850009‑2850300

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