Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,095,489 Δ1 bp 5.7% coding (392/705 nt) yggS → UPF0001 family protein, PLP‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,095,4890T. 5.7% 230.9 / 5.8 70coding (392/705 nt)yggSUPF0001 family protein, PLP‑binding
Reads supporting (aligned to +/- strand):  ref base T (33/32);  new base . (2/2);  total (36/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 = 1.00e+00

GCCCGTTGCAGTCTAATAAAAGCCGCCTGGTGGCAGAGCATTTCGACTGGTGTCATACCATCGACCGTTTGCGCATCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTGAACTGGACGAGCTGGCAGCTGCGGTCGCTGAACTACCGCGTTTACGTCTGCGCGGGTTGATGGCAATCCCTGCGCCTGAGTCAGAATATGTAAGGCAGTTTGAAGTTGC  >  NC_000913/3095345‑3095636
                                                                                                                                                |                                                                                                                                                   
gCCCGGTGCAGTCTAATAAAAGCCGCCTGGTGGCAGAGCATTTCGACTGGTGTCATACCATCGACCGTTTGCGCATCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAg                                                                                                                                                  <  8:186189/148‑1 (MQ=255)
  ccGTTGCAGTCTAATAAAAGCCGCCTGGTGGCAGAGCATTTCGACTGGTGTCATACCATCGACCGTTTGCGCATCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGtg                                                                                                                                                <  2:67240/148‑1 (MQ=255)
   cGTTGCAGTCTAATAAAAGCCGCCTGGTGGCAGAGCATTTCGACTGGTGTCATACCATCGACCGTTTGCGCATCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGtga                                                                                                                                               <  8:155812/148‑1 (MQ=255)
        cAGTCTAATAAAAGCCTCCTGGGGGCAGCTCATTTCGACTCGTGTCAGACCATCGACCGTTTGCGCAGAGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCACTCTTAACGTTCTGATTCAAATTAACATTAGTGATGaaa                                                                                                                                          <  2:285697/148‑1 (MQ=255)
         aGTCTAATAAAAGCCGCCTGGTGGCAGAGCATTTCGACTGGTGTCATACCATCGACCGTTTGCGCATCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGaaaa                                                                                                                                         <  5:126659/148‑1 (MQ=255)
            cTAATAAAAGCCGCCTGGTGGCAGAGCATTTCGACTGGTGTCATACCATCGACCGTTTGCGCATCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAg                                                                                                                                      >  3:8831/1‑148 (MQ=255)
              aataaAAGCCGCCTGGTGGCAGAGCATTTCGACTGGTGTCATACCATCGACCGTTTGCGCATCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGt                                                                                                                                     >  6:238105/1‑147 (MQ=255)
              aataaAAGCCGCCTGGTGGCAGAGCATTTCGACTGGTGTCATACCATCGACCGTTTGCGCATCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATAGTGATGAAAACAGt                                                                                                                                     >  3:279249/1‑146 (MQ=255)
              aataaAAGCCGCCTGGTGGCAGAGCATTTCGACTGGTGTCATACCATCGACCGTTTGCGCATCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATAGTGATGAAAACAGt                                                                                                                                     >  3:284307/1‑146 (MQ=255)
               ataaAAGCCGCCTGGTGGCAGAGCATTTCGACTGGTGTCATACCATCGACCGTTTGCGCATCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGt                                                                                                                                     <  8:357386/146‑1 (MQ=255)
               ataaAAGCCGCCTGGTGGCAGAGCATTTCGACTGGTGTCATACCATCGACCGTTTGCGCATCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTa                                                                                                                                    <  6:235238/147‑1 (MQ=255)
                  aaaGCCGCCTGGGGGCAGAGCATTTCGACTGGTGTCATACCATCGACCGTTTGCGCATCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGt                                                                                                                                 <  2:213856/147‑1 (MQ=255)
                    aGCCGCCTGGTGGCAGAGCATTTCGACTGGTGTCATACCATCGACCGTTTGCGCATCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCg                                                                                                                              <  7:5782/148‑1 (MQ=255)
                    aGCCGCCTGGTGGCAGAGCATTTCGACTGGTGTCATACCATCGACCGTTTGCGCATCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCg                                                                                                                              >  7:48464/1‑148 (MQ=255)
                      ccgccTGGTGGCAGAGCATTTCGACTGGTGTCATACCATCGACCGTTTGCGCATCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCggg                                                                                                                            <  3:315336/148‑1 (MQ=255)
                       cgccTGGTGGCAGAGCATTTCGACTGGTGTCATACCATCGACCGTTTGCGCATCGCTACCCGTCTCAACGACCAGCGCCCGGCAGCACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGa                                                                                                                           >  7:44605/1‑148 (MQ=255)
                             gtggCAGAGCATTTCGACTGGTGTCATACCATCGACCGTTTGCGCATCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCa                                                                                                                       <  6:370283/146‑1 (MQ=255)
                             gtggCAGAGCATTTCGACTGGTGTCATACCATCGACCGTTTGCGCATCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATAGTGATGAAAACAGTAAGTCCGGGATTCAAc                                                                                                                     <  4:279249/147‑1 (MQ=255)
                             gtggCAGAGCATTTCGACTGGTGTCATACCATCGACCGTTTGCGCATCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATAGTGATGAAAACAGTAAGTCCGGGATTCAAc                                                                                                                     <  4:284307/147‑1 (MQ=255)
                                    agCATTTCGACTGGTGTCATACCATCGACCGTTTGCGCATCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTg                                                                                                               >  8:50905/1‑147 (MQ=255)
                                     gCATTTCGACTGGTGTCATACCATCGACCGTTTGCGCATCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTGaa                                                                                                             >  1:179338/1‑148 (MQ=255)
                                         ttCGACTGGTGTCATACCATCGACCGTTTGCGCATCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTACCGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTGAACTg                                                                                                          >  4:226619/1‑147 (MQ=255)
                                             aCTGGTGTCATACCATCGACCGTTTGCGCATCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTGAACTGGACGa                                                                                                     <  7:96314/148‑1 (MQ=255)
                                                    tCATACCATCGACCGTTTGCGCATCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTGAACTGGACGAGCTGGc                                                                                               <  5:377851/147‑1 (MQ=255)
                                                    tCATACCATCGAACGTTTGCGCATCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTGAACTGGACGAGCTGGCa                                                                                              <  2:41550/148‑1 (MQ=255)
                                                       tACCATCGACCGTTTGCGCATCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGc                                                                                                                 <  7:211327/126‑1 (MQ=255)
                                                       tACCATCGACCGTTTGCGCATCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGc                                                                                                                 >  8:211327/1‑126 (MQ=255)
                                                         ccATCGACCGTTTGCGCATCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCg                                                                                                                              <  2:382462/111‑1 (MQ=255)
                                                         ccATCGACCGTTTGCGCATCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCg                                                                                                                              >  1:382462/1‑111 (MQ=255)
                                                           aTCGACCGTTTGCGCATCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTGAACTGGACGAGCTGGCAGCTGCgg                                                                                       <  4:8831/148‑1 (MQ=255)
                                                            tCGACCGTTTGCGCATCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTGAACTGGACGAGCTGGCAGCTGCGGt                                                                                      >  2:237041/1‑148 (MQ=255)
                                                             cGACCGTTTGCGCATCGCTACCCGTCTCAACGACCAGCGACCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTGAACTGGACGAGCTGGCAGCTGCGGTc                                                                                     >  2:83839/1‑148 (MQ=255)
                                                               aCCGTTTGCGCATCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCTCTCTTAACGTTCTGATTCAAATTAACATTAGtga                                                                                                                                               <  6:68704/88‑1 (MQ=255)
                                                               aCCGTTTGCGCATCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCTCTCTTAACGTTCTGATTCAAATTAACATTAGtga                                                                                                                                               >  5:68704/1‑88 (MQ=255)
                                                               aCCGTTTGCGCATCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTa                                                                                                                                    >  4:40891/1‑99 (MQ=255)
                                                               aCCGTTTGCGCATCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTa                                                                                                                                    <  3:40891/99‑1 (MQ=255)
                                                               aCCGTTTGCGCATCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCa                                                                                                                       <  4:294533/112‑1 (MQ=255)
                                                               aCCGTTTGCGCATCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCa                                                                                                                       >  3:294533/1‑112 (MQ=255)
                                                                         cATCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTGAACTGGACGAGCTGGCAGCTGCGGTCGCTGAACTACcg                                                                         <  7:3460/148‑1 (MQ=255)
                                                                           tCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGc                                                                                                                 >  6:357344/1‑106 (MQ=255)
                                                                           tCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGc                                                                                                                 <  5:357344/106‑1 (MQ=255)
                                                                              cTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTGAACTGGACGAGCTGGCAGCTGCGGTCGCTGAACTACCGCGttt                                                                    <  7:50905/148‑1 (MQ=255)
                                                                                      ctcAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTGAACTGGACGAGCTGGCAGCTGCGGTCGCTGAACTACCGCGTTTACGTCTgc                                                            >  4:222426/1‑148 (MQ=255)
                                                                                        cAACTACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTATGATTCAAATTAACAGTAGTGATGAAACCAGTAAGTCCGGGAGTCAACTGGCTGAACTGGCCGAGCTGGCAGCTGCGGTCGCTGAACTACCGCGTTTACGTCTgcgc                                                          >  7:289349/1‑148 (MQ=255)
                                                                                        cAACGCCCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTACGGGATTCAACTGGCTGAACTGGACGAGCTGGCTGCTGCGGTCGCTGCACTACCGCGGTTACGTCTgcgc                                                          >  7:289347/1‑148 (MQ=255)
                                                                                        cAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTGAACTGGACGAGCTGGCAGCt                                                                                           >  7:326788/1‑115 (MQ=255)
                                                                                        cAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTGAACTGGACGAGCTGGCAGCt                                                                                           <  8:326788/115‑1 (MQ=255)
                                                                                              ccAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTGAACTGGACGAGCTGGCAGCTGCGGTCGCTGAACTACCGCGTTTACGTCTGCGCGGGTTg                                                    <  6:174033/148‑1 (MQ=255)
                                                                                                  cgcCCGGCAGACCTTCCCCCTCTTACCGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCCACTGGCTGAACTGGACGAGCTGGCAGCTGCGGTCGCTGCACTACCGCGGTTACGTCTGCGCGGGTTGAAgg                                                >  6:24955/1‑148 (MQ=255)
                                                                                                  cgcCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTGAACTGGACGAGCTg                                                                                                 >  1:67508/1‑99 (MQ=255)
                                                                                                  cgcCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTGAACTGGACGAGCTg                                                                                                 <  2:67508/99‑1 (MQ=255)
                                                                                                     ccGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTGAACTGGACGAGCTGGCAGCTGCGGTCGCTGAACTACCGCGTTTACGTCTGCGCGGGTTGATGGCaa                                             >  2:171412/1‑148 (MQ=255)
                                                                                                      cGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTGAACTGGACGAGCTGGCAGCTGCGGTCGCTGAACTACCGCGTTTACGTCTGCGCGGGTTGATGGCAAt                                            >  1:143598/1‑148 (MQ=255)
                                                                                                          aGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGAGGAAAACAGTAAGTCCGGGATTCAACTGGCTGAACTGGACGAGCTGGCAGCTGCGGTCGCTGAACTACCGCGTTTACGTCTGCGCGGGTTGATGGCAATCCCt                                        <  5:24955/148‑1 (MQ=255)
                                                                                                             aCTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTGAACTGGACGAGCTGGCAGCTGCGGTCGCTGAACTACCGCGTTTACGTCTGCGCGGGTTGATGGCAATCCCTgcg                                     >  7:88119/1‑148 (MQ=255)
                                                                                                                  ccccTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTGAACTGGACGAGCTGGCAGCTGCGGTCGCTGAACTACCGCGTTTACGTCTGCGCGGGTTGATGGCAATCCCTGCGCCt                                  >  5:102157/1‑146 (MQ=255)
                                                                                                                   cccTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTGAACTGGACGAGCTGGCAGCTGCGGTCGCTGaa                                                                               <  4:142118/100‑1 (MQ=255)
                                                                                                                   cccTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTGAACTGGACGAGCTGGCAGCTGCGGTCGCTGaa                                                                               >  3:142118/1‑100 (MQ=255)
                                                                                                                     ctctTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTGAACTGGACGAGCTGGCAGCTGCGGTCGCTGAACTACCGCGTTTACGTCTGCGCGGGTTGATGGCAATCCCTGCGCCTGAGt                              >  6:149426/1‑147 (MQ=255)
                                                                                                                      tctTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTGAACTGGACGAGCTGGCAGCTGCGGTCGCTGAACTACCGCGTTTACGTCTGCGCGGGTTGATGGCAATCCCTGCGCCTGAGt                              >  1:98379/1‑146 (MQ=255)
                                                                                                                       ctTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTGAACTGGACGAGCTGGCAGCTGCGGTCGCTGAACTACCGCGTTTACGTCTGCGCGGGTTGATGGCAATCCCTGCGCCTGAGt                              >  4:10009/1‑145 (MQ=255)
                                                                                                                       ctTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTGAACTGGACGAGCTGGCAGCTGCGGTCGCTGAACTACCGCGTTTACGTCTGCGCGGGTTGATGGCAATCCCTGCGCCTGAGt                              >  5:338497/1‑145 (MQ=255)
                                                                                                                           aCGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTGAACTGGACGAGCTGGCAGCTGCGGTCGCTGAACTACCGCGTTTACGTCTGCGCGGGTTGATGGCAATCCCTGCGCCTGAGTCAGAat                        <  6:102157/147‑1 (MQ=255)
                                                                                                                             gTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTGAACTGGACGAGCTGGCAGCTGCGGTCGCTGAACTACCGCGTTTACGTCTGCGCGGGTTGATGGCAATCCCTGCGCCTGAGTCAGAATATg                     <  5:238105/148‑1 (MQ=255)
                                                                                                                               tCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTGAACTGGACGAGCTGGCAGCTGCGGTCGCTGa                                                                                <  2:165500/87‑1 (MQ=255)
                                                                                                                               tCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTGAACTGGACGAGCTGGCAGCTGCGGTCGCTGa                                                                                >  1:165500/1‑87 (MQ=255)
                                                                                                                               tCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTGAACTGGACGAGCTGGCAGCTGCGGTCGCTGAACTACCGCGTTTACGTCTGCGCGGGTTGATGGCAATCCCTGCGCCTGAGTCAGAATATGt                    >  7:75038/1‑147 (MQ=255)
                                                                                                                               tCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTGAACTGGACGAGCTGGCAGCTGCGGTCGCTGAACTACCGCGTTTACGTCTGCGCGGGTTGATGGCAATCCCTGCGCCTGAGTCAGAATATGt                    >  2:46798/1‑147 (MQ=255)
                                                                                                                                cTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTGAACTGGACGAGCTGGCAGCTGCGGTCGCTGAACTACCGCGTTTACGTCTGCGCGGGTTGATGGCAATCCCTGCGCCTGAGTCAGAATATGTa                   <  1:237041/147‑1 (MQ=255)
                                                                                                                                           tAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTGAACTGGACGAGCTGGCAGCTGCGGTCGCTGAACTACCGCGTTTACGTCTGCGCGGGTTGATGGCAATCCCTGCGCCTGAGTCAGAATATGTAAGGCAGTTTg         >  1:376265/1‑146 (MQ=255)
                                                                                                                                                ttAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTGAACTGGTAGAGCTGGCAGCTGCGGGCGCTGAACTACCGCGTTTACGTCTGCTCGGGTTGATGGCAATCCCTGCGCCTGAGTCAGAATATGTAAGGCAGTTTGAAGTTGc  <  3:226619/148‑1 (MQ=255)
                                                                                                                                                ttAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTGAACTGGACGAGCTGGCAGCTGCGGTCGCTGAACTACCGCGTTTACGTCTGcgcg                                                         <  2:228345/93‑1 (MQ=255)
                                                                                                                                                ttAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTGAACTGGACGAGCTGGCAGCTGCGGTCGCTGAACTACCGCGTTTACGTCTGcgcg                                                         >  1:228345/1‑93 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GCCCGTTGCAGTCTAATAAAAGCCGCCTGGTGGCAGAGCATTTCGACTGGTGTCATACCATCGACCGTTTGCGCATCGCTACCCGTCTCAACGACCAGCGCCCGGCAGAACTTCCCCCTCTTAACGTTCTGATTCAAATTAACATTAGTGATGAAAACAGTAAGTCCGGGATTCAACTGGCTGAACTGGACGAGCTGGCAGCTGCGGTCGCTGAACTACCGCGTTTACGTCTGCGCGGGTTGATGGCAATCCCTGCGCCTGAGTCAGAATATGTAAGGCAGTTTGAAGTTGC  >  NC_000913/3095345‑3095636

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: