Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,811,560 G→A 6.8% A209A (GCG→GCA ydjN → putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,811,5600GA 6.8% 182.4 / 2.9 59A209A (GCG→GCAydjNputative transporter
Reads supporting (aligned to +/- strand):  ref base G (22/33);  new base A (2/2);  total (24/35)
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.82e-01

TCTGAGCGTTCCGCAGCTGGTCTTGTCCTTTATCCCGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTAATGACCAAAGTGGTTGCAGGTTCTAACTTGCAAGACATCATC  >  NC_000913/1811413‑1811707
                                                                                                                                                   |                                                                                                                                                   
tctgaGCGTTCCGCAGCTGGTCTTGTCCTTTATCCCGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATgcgc                                                                                                                                                    >  5:117035/1‑149 (MQ=255)
    aGCGTTCCGCAGCTGGTCTTGTCCTTTATCCCGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCTTAGCTGCATTGAAACTGCTGAAGGATGATGCGcc                                                                                                                                                   <  1:187862/146‑1 (MQ=255)
                   gTCTTGTCCTTTATCCCGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCt                                                                                                                                  >  7:124791/1‑148 (MQ=255)
                   gTCTTGTCCTTTATCCCGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCt                                                                                                                                  <  3:49716/148‑1 (MQ=255)
                    tCTTGTCCTTTATCCCGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTa                                                                                                                                <  6:117035/149‑1 (MQ=255)
                      ttGTCCTTTATCCCGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCACCGAAAGGTGAACGCGTCTTAGc                                                                                                                              >  6:133955/1‑149 (MQ=255)
                      ttGTCCTTTATCCCGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCACCGAAAGGTGAACGCGTCTTAGc                                                                                                                              >  6:137213/1‑149 (MQ=255)
                        gTCCTTTATCCCGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAg                                                                                                                               <  8:16245/146‑1 (MQ=255)
                        gTCCTTTATCCCGAAAAACCCGTTTGACGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAACCTGATGAAGCATGATGCGCCGAAAGGTGAACGCCTATTAGCCg                                                                                                                            >  5:116011/1‑149 (MQ=255)
                         tCCTTTATCCCGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGc                                                                                                                           >  8:169457/1‑149 (MQ=255)
                           cTTTATCCCGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTa                                                                                                                         <  8:118535/149‑1 (MQ=255)
                           cTTTATCCCGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTa                                                                                                                         >  4:58038/1‑149 (MQ=255)
                            tttATCCCGAAAACCCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTAGCCGCATTCCTCGGCGTATCTGCATTGAAACTGCTGAAGTATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTAt                                                                                                                        <  7:130511/149‑1 (MQ=255)
                                tCCCGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGa                                                                                                                     >  3:52324/1‑148 (MQ=255)
                                  ccGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATa                                                                                                                   <  4:194170/148‑1 (MQ=255)
                                  ccGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATAc                                                                                                                  >  5:201826/1‑149 (MQ=255)
                                  ccGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATAc                                                                                                                  >  5:112319/1‑149 (MQ=255)
                                  ccGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATAc                                                                                                                  <  4:36976/149‑1 (MQ=255)
                                  ccGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATAc                                                                                                                  >  3:36976/1‑149 (MQ=255)
                                   cGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATAcc                                                                                                                 >  5:64897/1‑149 (MQ=255)
                                              tttGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTg                                                                                                      <  4:213768/149‑1 (MQ=255)
                                                ttcagATCTTACCTGAGCCTATGCGACGTCAATTATCAGCGTTGTAATTTTTGCCGCATTCCTCGGCGTAGCTGGATTGAAACTTCTGAAGGATGATGCGCCGAAAGTTGAACGCGTCTTAGCCGCTATCTATACCCTAGAAAGCTgg                                                                                                     <  4:22289/144‑1 (MQ=255)
                                                  ccGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGt                                                                                                   <  7:169457/148‑1 (MQ=255)
                                                          aCCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGtgat                                                                                                <  1:100267/143‑1 (MQ=255)
                                                          aCCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGtga                                                                                                 <  2:37267/142‑1 (MQ=255)
                                                          aCCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGaaa                                                                                            <  8:134907/147‑1 (MQ=255)
                                                            cGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTgg                                                                                        >  4:126748/1‑149 (MQ=255)
                                                               aGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTc                                                                                     <  3:58038/149‑1 (MQ=255)
                                                                   aaTCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCt                                                                                 >  7:37394/1‑149 (MQ=255)
                                                                     tCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTgg                                                                               >  6:179148/1‑149 (MQ=255)
                                                                       cGACGTCAAATATCAGCGGGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAATGGTGAACGCGTCTTAGCCGCTATAGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGGCTGGTc                                                                             <  6:195726/149‑1 (MQ=255)
                                                                        gACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCa                                                                            <  2:35592/149‑1 (MQ=255)
                                                                              aaTTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGt                                                                      >  6:176181/1‑149 (MQ=255)
                                                                               aTTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGtt                                                                     <  7:213846/149‑1 (MQ=255)
                                                                                  aTCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGa                                                                   <  5:234141/148‑1 (MQ=255)
                                                                                            tAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGATCCCTTACGGc                                                        >  7:228833/1‑149 (MQ=255)
                                                                                            tAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGc                                                        >  1:51551/1‑149 (MQ=255)
                                                                                               tttttGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGtt                                                     <  6:64897/149‑1 (MQ=255)
                                                                                                     ccGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGct                                               <  5:176181/149‑1 (MQ=255)
                                                                                                      cGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGctc                                              <  5:165600/149‑1 (MQ=255)
                                                                                                         aTTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTaa                                           <  6:112319/149‑1 (MQ=255)
                                                                                                          ttCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTAAt                                          <  7:232243/149‑1 (MQ=255)
                                                                                                             cTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTAATGa                                        <  8:124791/148‑1 (MQ=255)
                                                                                                             cTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGGTAACGCGTCTTAGCCGCGATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGTTCTAATGAc                                       <  5:35768/149‑1 (MQ=255)
                                                                                                                          cATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTAATGACCAAAGTGGTTg                            <  5:179148/147‑1 (MQ=255)
                                                                                                                          cATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTAATGACCAAAGTGGTTGc                           <  4:52324/148‑1 (MQ=255)
                                                                                                                           aTTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTAATGACCAAAGTGGTTGCAg                         >  2:33213/1‑149 (MQ=255)
                                                                                                                           aTTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCAGGCTCTAATGACCAAAGTCGTTGCAg                         >  2:91703/1‑149 (MQ=255)
                                                                                                                             tGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTaa                                           >  4:194648/1‑129 (MQ=255)
                                                                                                                             tGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTaa                                           <  3:194648/129‑1 (MQ=255)
                                                                                                                                      tgAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTAATGACCAAAGTGGTTGCAGGTTCTAACTTg              >  6:177838/1‑149 (MQ=255)
                                                                                                                                        aaGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTAATGACCAAAGTGGtt                             <  7:43441/132‑1 (MQ=255)
                                                                                                                                        aaGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTAATGACCAAAGTGGtt                             >  8:43441/1‑132 (MQ=255)
                                                                                                                                         aGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGATCCCTTACGGCGTTCTGGCTCTAATGACCAAAGTGGTTGCAGGTTCTAACTTGCaa           <  8:228833/149‑1 (MQ=255)
                                                                                                                                         aGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTAATGACCAAAGTGGTTGCAGGTTCTAACTTGCaa           <  6:201826/149‑1 (MQ=255)
                                                                                                                                         aGGATGATGCACCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTAATGACCAAAGTGGTTGCAGGTTCTAACTTGCaa           <  5:133955/149‑1 (MQ=255)
                                                                                                                                         aGGATGATGCACCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTAATGACCAAAGTGGTTGCAGGTTCTAACTTGCaa           <  5:137213/149‑1 (MQ=255)
                                                                                                                                           gatgatGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTAATGACCAAAGTGGTTGCAGGTTCTAACTTGCAAGa         >  7:4687/1‑149 (MQ=255)
                                                                                                                                                tgCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTAATGACCAAAGTGGTTg                            <  6:2924/125‑1 (MQ=255)
                                                                                                                                                tgCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTAATGACCAAAGTGGTTg                            >  5:2924/1‑125 (MQ=255)
                                                                                                                                                   gcCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTAATGACCAAAGTGGTTGCAGGTTCTAACTTGCAAGACatcatc  >  4:102116/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TCTGAGCGTTCCGCAGCTGGTCTTGTCCTTTATCCCGAAAAACCCGTTTGCCGATCTTACCGGAGCCAATCCGACGTCAATTATCAGCGTGGTAATTTTTGCCGCATTCCTCGGCGTAGCTGCATTGAAACTGCTGAAGGATGATGCGCCGAAAGGTGAACGCGTCTTAGCCGCTATCGATACCCTACAAAGCTGGGTGATGAAACTGGTTCGCCTGGTCATGCAGTTGACCCCTTACGGCGTTCTGGCTCTAATGACCAAAGTGGTTGCAGGTTCTAACTTGCAAGACATCATC  >  NC_000913/1811413‑1811707

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: