Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,340,505 C→T 6.6% A88V (GCG→GTG)  yciS → DUF1049 family inner membrane protein, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,340,5050CT 6.6% 185.6 / 2.6 60A88V (GCG→GTG) yciSDUF1049 family inner membrane protein, function unknown
Reads supporting (aligned to +/- strand):  ref base C (32/24);  new base T (2/2);  total (34/26)
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.88e-01

CCGTATTTCCACATTGCTGGCGGTATTGTTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTGTTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGCAACAAAACAAGCAAGATG  >  NC_000913/1340359‑1340648
                                                                                                                                                  |                                                                                                                                               
ccGTATTTCCACAGTGCTGGCGGTATTGTTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTGTTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGCCTGGCACACCAGCTTTCACCCGCGa                                                                                                                                               >  3:93925/1‑149 (MQ=255)
       tCCACATGCCGGGCGGTCTTGTTTGCTGCGGGGTGTGCGATAGGTTGGGTGATTTGGGGCCGGTTCTGGCGGCGCGTTCGTCGTTCCCTGGCGCGCGCTGAACGTACAATAAAGCGCCGGGAAAACCCGCTTTCACCCGCGACTGCCGt                                                                                                                                        <  4:150684/149‑1 (MQ=255)
        ccACATTGCTCGCGGAATTCTATGCTCCGGGGTTTGCTATCTGTTCGTTCATTTGTGTCCACTTCTGGCTTCGAGTTCGTGTTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGCAATACCAGCTTTCACCCGCGACTGACGt                                                                                                                                        <  4:155785/148‑1 (MQ=255)
          acaTTGCTGGCGGTATTGTTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTCGTGGCCTGTTCTGGCTGCGAGTTCGTGTTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGc                                                                                                                                     >  3:142967/1‑149 (MQ=255)
           caTTGCTGGCGGTATTGTTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTGTTGCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCt                                                                                                                                    >  8:104710/1‑149 (MQ=255)
           caTTGCTGGCGGTATTGTTTGCTGCGGGGATTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTGTTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCt                                                                                                                                    >  4:36395/1‑149 (MQ=255)
                cTGGCGGTATTGTTTGCTGCGGGGTTTGCTATCGGGTGGGTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTGTTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAg                                                                                                                                <  2:115367/148‑1 (MQ=255)
                     ggTATTGTTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGAGAGTTCGTGTTTCCCTGGCGCGCGCTGAACGTAATATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGc                                                                                                                          >  1:63527/1‑149 (MQ=255)
                         ttGTTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTGTTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTc                                                                                                                      >  1:224474/1‑149 (MQ=255)
                            tttGCTGCGGGGTTTGCTCTCGGTTGGTTGTTTTGTGGCGTGTTCTGGATGCGAGTTCGTG‑TTACCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCa                                                                                                                  <  7:330833/149‑1 (MQ=255)
                                   cGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTGTTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAgcggcg                                                                                                            <  6:50230/149‑1 (MQ=255)
                                           gCTATCGGTTGGTTGATTTGTGGACTGTTCTGGCTGCGAGTTAGTGTTTCCCTGGCGAGCGCTGAACCTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAAt                                                                                                     >  1:127434/1‑148 (MQ=255)
                                            cTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTGTTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATaa                                                                                                   >  5:233143/1‑149 (MQ=255)
                                            cTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTGTTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATaa                                                                                                   <  3:233091/149‑1 (MQ=255)
                                              aTCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTGTTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGTGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACt                                                                                                 >  3:144132/1‑149 (MQ=255)
                                              aTCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTGTTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACACGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACt                                                                                                 >  5:137798/1‑149 (MQ=255)
                                                 ggttggttGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTGTTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACttt                                                                                               <  3:319893/148‑1 (MQ=255)
                                                      gttgCTTGGTGGCCTGTTCTGGCTGCGAGTTCGTGGTTCCCTGGCGCGCGCTGAACGTAAAATACAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGGGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGc                                                                                         <  4:93925/149‑1 (MQ=255)
                                                                 gCCTGTTCTGGCTGCGAGTTCGTGTTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAgttgttg                                                                              >  6:152254/1‑149 (MQ=255)
                                                                 gCCTGTTCTGGCTGCGAGTTCGTGTTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAgttgttg                                                                              <  4:142967/149‑1 (MQ=255)
                                                                 gCCTGTTCTGGCTGCGAGTTCGTGTTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAgttgtt                                                                               <  4:85995/148‑1 (MQ=255)
                                                                 gCCTGTTCTGGCTGCGAGTTCGTGTTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAgt                                                                                   >  7:334638/1‑144 (MQ=255)
                                                                 gCCTGTTCTGGCTGCGAGTTCGTGTTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAgt                                                                                   <  8:334638/144‑1 (MQ=255)
                                                                     gTTCTGGCTGCGAGTTCGTGTTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAg                                                                                                         >  7:236164/1‑118 (MQ=255)
                                                                     gTTCTGGCTGCGAGTTCGTGTTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAg                                                                                                         <  8:236164/118‑1 (MQ=255)
                                                                        cTGGCTGCGAGTTCGTGTTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGc                                                                       >  6:129492/1‑149 (MQ=255)
                                                                         tGGCTGCGAGTTCGTGTTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCt                                                                      >  1:263492/1‑149 (MQ=255)
                                                                           gCTGCGAGTTCGTGTTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCtt                                                                     <  5:129492/148‑1 (MQ=255)
                                                                            cTGCGAGTTCGTGTTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCtttt                                                                   >  1:143064/1‑149 (MQ=255)
                                                                               cGAGTTCGTGTTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGc                                                                 <  3:36395/148‑1 (MQ=255)
                                                                                gAGTTCGTGTTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCt                                                               >  7:21370/1‑149 (MQ=255)
                                                                                  gTTCGTGTTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTTCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTg                                                              >  3:290063/1‑148 (MQ=255)
                                                                                     cGTGTTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGAGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAga                                                          >  8:242240/1‑148 (MQ=255)
                                                                                       tgtTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGCCGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCg                                                        >  1:14898/1‑149 (MQ=255)
                                                                                         tttCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCt                                                      <  7:227427/149‑1 (MQ=255)
                                                                                            cccTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGTGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGc                                                    <  4:144132/148‑1 (MQ=255)
                                                                                              cTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTa                                                 >  8:264355/1‑149 (MQ=255)
                                                                                                     gcgcTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAATGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTgg                                          <  7:93269/149‑1 (MQ=255)
                                                                                                      cgcTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTTCTTTTGCCTGTAGCCGCTGCCTATGGCTGGt                                         <  6:137798/149‑1 (MQ=255)
                                                                                                      cgcTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGACTGGt                                         <  4:323363/149‑1 (MQ=255)
                                                                                                         tGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGtata                                      >  1:87092/1‑149 (MQ=255)
                                                                                                         tGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGtata                                      >  7:119416/1‑149 (MQ=255)
                                                                                                          gAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGTGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCtttt                                                                   >  7:185290/1‑119 (MQ=255)
                                                                                                          gAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGTGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCtttt                                                                   <  8:185290/119‑1 (MQ=255)
                                                                                                          gAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTatat                                     >  7:344655/1‑149 (MQ=255)
                                                                                                           aaCGTTAAATAAAGCGACTGGAATACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGAGAAGGAATATCTTTCTATGCTGGAGTTTTTGTTTCTGCTTTTGCCTGTATCCGCTGCCTATGGCTGGTATTTg                                    <  2:127434/149‑1 (MQ=255)
                                                                                                              gTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGc                                 >  5:64429/1‑149 (MQ=255)
                                                                                                              gTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGc                                 >  3:1330/1‑149 (MQ=255)
                                                                                                                           aCTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTgcg                     >  1:335556/1‑148 (MQ=255)
                                                                                                                                    ccAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGGTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGCaacaaaac            >  5:174830/1‑148 (MQ=255)
                                                                                                                                     cAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGATGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGCaacaaaacaa          >  6:93819/1‑149 (MQ=255)
                                                                                                                                       gCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGCAACAAAAcaag         >  7:195991/1‑148 (MQ=255)
                                                                                                                                       gCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGCAACAAAAcaag         >  5:72018/1‑148 (MQ=255)
                                                                                                                                           tCACCCGCGACTGACGTGGCTGTAGTTCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGCAACAAAAcaagcaa      <  4:290063/147‑1 (MQ=255)
                                                                                                                                           tCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGCAACAAAAcaagcaag     <  8:21370/148‑1 (MQ=255)
                                                                                                                                           tCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGCAACAAAAcaagcaa      <  6:72018/147‑1 (MQ=255)
                                                                                                                                           tCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGCAACAAAAcaagcaa      <  2:87092/147‑1 (MQ=255)
                                                                                                                                             aCCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGCAACAAAACAAGCAAGAt   <  1:268653/148‑1 (MQ=255)
                                                                                                                                              cccGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGCAACAAAACAAGCAAGATg  >  1:345946/1‑148 (MQ=255)
                                                                                                                                               ccGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGCAACAAAACAAGCGAGATg  >  7:282724/1‑147 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
CCGTATTTCCACATTGCTGGCGGTATTGTTTGCTGCGGGGTTTGCTATCGGTTGGTTGATTTGTGGCCTGTTCTGGCTGCGAGTTCGTGTTTCCCTGGCGCGCGCTGAACGTAAAATAAAGCGACTGGAAAACCAGCTTTCACCCGCGACTGACGTGGCTGTAGTGCCGCACTCGTCAGCGGCGAAGGAATAACTTTCTATGCTGGAGTTGTTGTTTCTGCTTTTGCCTGTAGCCGCTGCCTATGGCTGGTATATGGGCCGCAGAAGTGCGCAACAAAACAAGCAAGATG  >  NC_000913/1340359‑1340648

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