Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,665,193 T→A 5.1% Y268* (TAT→TAA yfhR → putative S9 family prolyl oligopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,665,1930TA 5.1% 255.0 / 2.9 78Y268* (TAT→TAAyfhRputative S9 family prolyl oligopeptidase
Reads supporting (aligned to +/- strand):  ref base T (38/36);  new base A (2/2);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

TACTCATTCACGGTAAAGCTGATCACGTTATCCCATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGATAATCATTCACAAAACCACCTTAAGACATGCTAATCCACTGGTCAGAACAGTTTAAGATGAGAAAAATTCTGT  >  NC_000913/2665048‑2665336
                                                                                                                                                 |                                                                                                                                               
tACTCATTCACGGTAAAGCTGATCACGTTATCCCATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATcg                                                                                                                                               <  4:265210/148‑1 (MQ=255)
 aCTCATTCACGGTAAAGCTGATCACGTTATCCCATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATGCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATcgc                                                                                                                                              <  5:198232/148‑1 (MQ=255)
     aTTCACGGTAAAGCTGATCACGTTATCCCATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGaa                                                                                                                                           <  7:131663/147‑1 (MQ=255)
     aTTCACGGTAAAGCTGATCACGTTATCCCATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGaa                                                                                                                                           >  8:70444/1‑147 (MQ=255)
     aTTCACGGTAAAGCTGATCACGTTATCCCATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCAGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGaa                                                                                                                                           >  7:55561/1‑147 (MQ=255)
      ttCACGGTAAAGCTGATCACGTTATCCCATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGaa                                                                                                                                           >  6:30490/1‑146 (MQ=255)
      ttCACGGTAAAGCTGATCACGTTATCCCATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGaa                                                                                                                                           >  7:376726/1‑146 (MQ=255)
       tCACGGTAAAGCTGATCACGTTATCCCATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGaa                                                                                                                                           >  1:49338/1‑145 (MQ=255)
         atggTATAGCTGATCACGTTATCCCATGGCAGCACATCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGTTCCTAATCCCGGATGGCGAACACATTGATTCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGAtg                                                                                                                                     <  1:254632/147‑1 (MQ=255)
          cGGTAAAGCTGATCACGTTATCCCATGGCAGCACAGCGAACAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGAtgg                                                                                                                                    <  6:304709/149‑1 (MQ=255)
          cGGTAAAGCTGATCACGTTATCCCATGGCAGCACAGCGAAAAGTGGTATATACTGGCAAAAGAGCGAAAAGGGTTGATCCTAATCCCGTATGGCGAACACATTGATGGTTTTTCCGATCGTCACGGCTATGTTTATCGCGAACAGAtgg                                                                                                                                    <  7:296724/149‑1 (MQ=255)
            gTAAAGCTGATCACGTTATCCCATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGAtggtg                                                                                                                                  <  3:279978/149‑1 (MQ=255)
                  cTGATCACGTTATCCCATGGCAGTACAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACttt                                                                                                                            <  4:42203/149‑1 (MQ=255)
                  cTGATCACGTTATCCCATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACttt                                                                                                                            <  5:30490/149‑1 (MQ=255)
                      tCACGTTATCCCATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGAtggtgg                                                                                                                                 >  1:126684/1‑140 (MQ=255)
                      tCACGTTATCCCATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGAtggtgg                                                                                                                                 <  2:126684/140‑1 (MQ=255)
                         cGTTATCCCATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCt                                                                                                                       >  8:7618/1‑147 (MQ=255)
                         cGTTATCCCATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGAGGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTa                                                                                                                     <  1:257438/149‑1 (MQ=255)
                          gTTATCCCATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAg                                                                                                                    >  6:272812/1‑149 (MQ=255)
                           ttATCCCATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAgc                                                                                                                   <  2:17237/149‑1 (MQ=255)
                                ccATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGt                                                                                                               <  2:286641/148‑1 (MQ=255)
                                  aTGGCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGAtggtg                                                                                                                                  <  1:45860/127‑1 (MQ=255)
                                  aTGGCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGAtggtg                                                                                                                                  >  2:45860/1‑127 (MQ=255)
                                  aTGGCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTg                                                                                                             >  7:355934/1‑148 (MQ=255)
                                   tGGCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGaa                                                                                                           >  1:123773/1‑149 (MQ=255)
                                   tGGCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGCCGAAGTTTATCGCGAACACATGGAGGACTTAATCCATACCTCTTTGa                                                                                                            >  4:186272/1‑148 (MQ=255)
                                       agcaCAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTGGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGGTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCg                                                                                                       <  6:116225/149‑1 (MQ=255)
                                          acaGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCg                                                                                                       >  2:124597/1‑146 (MQ=255)
                                              cGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATCACGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAAc                                                                                                 >  4:59602/1‑148 (MQ=255)
                                                         ataGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCTCGAACAGATGGTGGACTTTATCCTTAGCGCTTTGAATCCGCATAACTAACCAATGATc                                                                                     >  6:292586/1‑149 (MQ=255)
                                                             ccTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCCTGATCGCTa                                                                                 <  2:17276/149‑1 (MQ=255)
                                                             ccTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTa                                                                                 <  8:376726/149‑1 (MQ=255)
                                                              cTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAg                                                                                >  3:145264/1‑149 (MQ=255)
                                                                ggCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCa                                                                              >  7:209885/1‑149 (MQ=255)
                                                                    aaaGAGCCAAAACGGTTGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGAt                                                                          <  7:70444/149‑1 (MQ=255)
                                                                         gCCAAAACGGTTGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGATAATc                                                                      <  3:186272/148‑1 (MQ=255)
                                                                          ccAAAACGGTTGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGaa                                                                                                  >  4:227672/1‑119 (MQ=255)
                                                                          ccAAAACGGTTGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGaa                                                                                                  <  3:227672/119‑1 (MQ=255)
                                                                           cAAAACGGTTGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGATAATCAtt                                                                   <  5:272812/149‑1 (MQ=255)
                                                                            aaaaCGGTTGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGATAATCAtt                                                                   >  6:15674/1‑148 (MQ=255)
                                                                                  gTTGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATc                                                                                     >  6:244606/1‑124 (MQ=255)
                                                                                  gTTGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATc                                                                                     <  5:244606/124‑1 (MQ=255)
                                                                                  gTTGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGATAATCATTCACaaaa                                                            >  3:175200/1‑149 (MQ=255)
                                                                                   ttGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGATAATCATTCACAAAac                                                           >  1:35265/1‑149 (MQ=255)
                                                                                     gATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATCTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGATAATCATTCACAAAacca                                                         >  6:216204/1‑149 (MQ=255)
                                                                                      aTCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGATAATCATTCACAAAaccac                                                        >  8:313651/1‑149 (MQ=255)
                                                                                        ccTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTAACGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGATAATCATTCACAAAaccacc                                                       >  2:320138/1‑148 (MQ=255)
                                                                                        ccTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTAACGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGATAATCATTCACAAAACCACCt                                                      >  2:325300/1‑149 (MQ=255)
                                                                                         cTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGATAATCATTCACAAAACCACCtt                                                     <  8:355934/149‑1 (MQ=255)
                                                                                            aTCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTAACGCTAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGATAATCATTCACAAAACCACCTaaa                                                   <  1:320138/148‑1 (MQ=255)
                                                                                            aTCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTAACGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGATAATCATTCACAAAACCACCTCaa                                                   <  1:325300/148‑1 (MQ=255)
                                                                                                        gAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGATAATCATTc                                                                  >  4:297273/1‑121 (MQ=255)
                                                                                                        gAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGATAATCATTc                                                                  <  3:297273/121‑1 (MQ=255)
                                                                                                        gAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGATAATCATTCACAAAACCACCTTAAGACATGCTAATc                                       >  1:223571/1‑148 (MQ=255)
                                                                                                         aaCACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGATAATCATTCACAAAACCACCTTAAGACATGCTAATCCa                                     <  2:123773/149‑1 (MQ=255)
                                                                                                          acacATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGATAATCATTCACAAAACCACCTTAAGACATGCTAATCCAc                                    >  8:87388/1‑149 (MQ=255)
                                                                                                            acaTTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGATAATCATTCAGAAAACCACCTTAAGACATGCTAATCCACTg                                  <  8:228770/149‑1 (MQ=255)
                                                                                                                 gATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGATAATCAtt                                                                   <  2:16887/111‑1 (MQ=255)
                                                                                                                 gATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGATAATCAtt                                                                   >  1:16887/1‑111 (MQ=255)
                                                                                                                  aTGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGATAATCATTCACAAAACCACCTTAAGACATGCTAATCCACTGGTCAg                             <  1:124597/148‑1 (MQ=255)
                                                                                                                   tGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGATAATCATTCACAAAACCACCTTAAGACATGCTAATCCACTGGTCAGaa                           <  4:175200/149‑1 (MQ=255)
                                                                                                                   tGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGATAGCACGATAATCATTCACAAAACCACCTTAAGACATGCTAATCCACTGGTCAGaa                           <  5:292586/149‑1 (MQ=255)
                                                                                                                     cTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGATAATCATTCACAAAACCACCTTAAGACa                                               <  2:110928/127‑1 (MQ=255)
                                                                                                                     cTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGATAATCATTCACAAAACCACCTTAAGACa                                               >  1:110928/1‑127 (MQ=255)
                                                                                                                      tttttCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGATAATCATTCACAAAACCACCTTAAGACATGCTAATCCACTGGTCAGAACAg                        >  3:95238/1‑149 (MQ=255)
                                                                                                                      tttttCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGATAATCATTCACAAAACCACCTTAAGACATGCTAATCCACTGGTCAGAACAg                        >  7:375697/1‑149 (MQ=255)
                                                                                                                            cGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGATAATCATTCACAAAACCACCTTAAGACATGCTAATCCACTGGTCAGAACAGTTTAAg                  <  7:313651/149‑1 (MQ=255)
                                                                                                                             gATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGATAATCATTCACAAAACCACCTTAAGACATGCTAATCCACTGGTCAGAACAGttt                     >  1:356365/1‑145 (MQ=255)
                                                                                                                              aTCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGATAATCATTCACAAAACCACCTTAAGACATGCTAATCCACTGGTCAGAACAGTTTAAGAt                <  7:7618/149‑1 (MQ=255)
                                                                                                                              aTCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATACTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGATAATCATTCACAAAACCACCTTAAGACATGCTAATCCACTGGTCAGAACAGTTTAAGAt                >  1:332135/1‑149 (MQ=255)
                                                                                                                                 gTCACGTCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACTATAATCATTCACAAAACCACCTTAAGACATGCTAATCCACTGGTCAGAACAGTTTAAGATgag             <  3:59602/149‑1 (MQ=255)
                                                                                                                                 gTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGATAATCATTCACAAAACCACCTTAAGACATGCTAATCCACTGGTCAGAACAGTTTAAGATgag             <  5:15674/149‑1 (MQ=255)
                                                                                                                                  tCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGATAATCATTCACAAAACCACCTTAAGACATGCTAATCCACTGGTCAGAACAGTTTAAGATgag             >  3:273897/1‑148 (MQ=255)
                                                                                                                                    aCGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGATAATCATTCACaaaa                                                            <  2:60826/99‑1 (MQ=255)
                                                                                                                                    aCGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGATAATCATTCACaaaa                                                            >  1:60826/1‑99 (MQ=255)
                                                                                                                                      ggCGATGTTTATGGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGATAATCATTCACAAAACCACCTTAAGACATGCTAAACCACTGGTCAGAACAGTTTAAGTTTAGaaaaa        >  4:264596/1‑149 (MQ=255)
                                                                                                                                      ggCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGATAATCATTCACAAAACCACCTTAAGACATGCTAATCCACTGGTCAGAACAGTTTAAGATGAGaaaaa        >  2:337555/1‑149 (MQ=255)
                                                                                                                                            gTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGATAATCATTCACAAAACCACCTTAAGACATGCTAATCCACTGGTCAGAACAGTTTAAGATGAGAAAAATTCtgt  >  7:16226/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
TACTCATTCACGGTAAAGCTGATCACGTTATCCCATGGCAGCACAGCGAAAAGTTGTATAGCCTGGCAAAAGAGCCAAAACGGTTGATCCTAATCCCGGATGGCGAACACATTGATGCTTTTTCCGATCGTCACGGCGATGTTTATCGCGAACAGATGGTGGACTTTATCCTTAGCGCGTTGAATCCGCAGAACTAACCCATGATCGCTAGCACGATAATCATTCACAAAACCACCTTAAGACATGCTAATCCACTGGTCAGAACAGTTTAAGATGAGAAAAATTCTGT  >  NC_000913/2665048‑2665336

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: