Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 214,739 T→C 5.7% M150T (ATG→ACG)  yaeQ → PDDEXK superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913214,7390TC 5.7% 212.9 / 2.2 70M150T (ATG→ACG) yaeQPDDEXK superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (33/33);  new base C (2/2);  total (35/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 = 9.19e-01

AAGTGGCGCTGTTTACCTATAATAGTCGGGCGGCGCAAATCTGGTGGCAGCAAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAACCGCCTGGCAACAACCTTCATGATTGTGATTTCCCGACATGTTGCTATTCCCGATGGTGAGCTTGAGATCACC  >  NC_000913/214592‑214886
                                                                                                                                                   |                                                                                                                                                   
aaGTGGCGCTGTTTACCTATAATAGTCGGGCGGCGCAAATCTGGTGGCAGCAAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCAt                                                                                                                                                     >  2:318914/1‑148 (MQ=255)
                cTATAATAGTCGGGCGGCGCAAATCTGGTGGCAGCAAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGCTATCTGGACGATGAACAACTGGGCTAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGa                                                                                                                                    >  6:123414/1‑149 (MQ=255)
                  ataataGTCGGGCGGCGCAAATCTGGTGGCAGCAAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGAtt                                                                                                                                  >  7:28684/1‑149 (MQ=255)
                   taataGTCGGGCGGCGCAAATCTGGTGGCAGCAAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGAAGCTGCAGGCAACGAtt                                                                                                                                  >  5:308988/1‑148 (MQ=255)
                     ataGTCGGGCGGCGCAAATCTGGTGGCAGCAAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCa                                                                                                                                <  4:175089/148‑1 (MQ=255)
                       aGTCGGGCGGCGCAACTCCGGTGGCAGCAAACTCCGCGCAAATGTGTGCAGTTTGCCACGCTTTCCGGCTGGTATCTGGACGAGGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGa                                                                                                                             <  2:71541/149‑1 (MQ=255)
                        gTCGGGCGGCGCAAATCTGGTGGCAGCAAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGa                                                                                                                             <  3:171441/148‑1 (MQ=255)
                        gTCGGGCGGCGCAAATCTGGTGGCAGCAAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGa                                                                                                                             <  5:164188/148‑1 (MQ=255)
                         tCGGGCGGCGCAAATCTGGTGGCAGCAAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATg                                                                                                                           >  1:129368/1‑149 (MQ=255)
                               ggcgCAAATCTGGTGGCAGCAAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGa                                                                                                                     <  8:1541/149‑1 (MQ=255)
                                gcgcAAATCTGGTGGCAGCAAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGAt                                                                                                                            <  8:196363/141‑1 (MQ=255)
                                gcgcAAATCTGGTGGCAGCAAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGAt                                                                                                                            >  7:196363/1‑141 (MQ=255)
                                   cAAATCTGGTGGCAGCAAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTg                                                                                                                 <  3:43039/149‑1 (MQ=255)
                                    aaaTCTGGTGGCAGCAAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGc                                                                                                                         <  2:240999/140‑1 (MQ=255)
                                    aaaTCTGGTGGCAGCAAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGc                                                                                                                         >  1:240999/1‑140 (MQ=255)
                                         tggtggCAGCAAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATc                                                                                                           >  7:324080/1‑149 (MQ=255)
                                               cagcaAAATCAGAGCAAATGTGTGCAGTTTGCGAATCTTTCCGTCTGGTATCGGGACGATGACCAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCACGACGCTGCAGGCAACGCGTCAGGATGGCGTGATCTGGTGATCGgatgc                                                                                                     >  6:223317/1‑148 (MQ=255)
                                               cagcaAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTTGTTATCGgatga                                                                                                     <  7:92751/149‑1 (MQ=255)
                                               cagcaAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACGTGCGAAAGTAAGCGCCTTTGCCGATCGTACCACGACGCAGCCGGCAACGATTCAGGATGGCGTGATCTGGTTATCGgatg                                                                                                      >  6:223318/1‑148 (MQ=255)
                                                         agagCAAATGTGTGCAGGTTGCCAATCTTTCCGTCTGGTCTCTTGACGATGAACAACTGTCGAAAGTAAGCGACTTGGCCGATCGTACCACGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGaataat                                                                                           <  5:223317/149‑1 (MQ=255)
                                                         agagCAAATGTGGGCAGTTTGCCAAGCTTTCCGTCTGGTATCTTGACGAAGAACAACTTGCGAAAGTAAGCGCCTTTGACGATCGTACCACGACGCTGCAGGCAACTCTTCAGGATGGCGTGATCTGGTTATCGGATGATAAGaataag                                                                                           <  5:223318/149‑2 (MQ=255)
                                                          gagCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTTGTTATCGGATGATAAGAATAATc                                                                                          >  8:96500/1‑149 (MQ=255)
                                                             cAAATGTGTGCAGTTTGCCAATCTTTCCGTCTCGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATAGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTgg                                                                                       <  5:123414/149‑1 (MQ=255)
                                                             aaaaTGTGTGCAGTTTGCCAATCTTTCCGTCTTGTATCTGGACGATTAACAACTGTCGAAATTAAGCGCATTTGCCGATCGTACCATGACGCTGCATGCAACGATTCAGGATGTCGTTATCTGGTTATCGGATGATAAGAATAATCTgg                                                                                       <  6:308988/148‑1 (MQ=255)
                                                                 tgtgtgCAGTTTGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCGTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGt                                                                                   <  6:295018/149‑1 (MQ=255)
                                                                   tgtgCAGTTTGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGa                                                                                 <  8:324080/149‑1 (MQ=255)
                                                                    gtgCAGTTTGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGaa                                                                                >  1:333298/1‑149 (MQ=255)
                                                                     tgCAGTTTGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGaa                                                                                <  7:26841/148‑1 (MQ=255)
                                                                     tgCAGTTTGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCATTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAAc                                                                               <  6:248614/149‑1 (MQ=255)
                                                                       cAGTTTGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTg                                                                                                                                              >  4:256434/1‑84 (MQ=255)
                                                                        aGTTTGCCAATCTTTCCTTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCaa                                                                                                                                       <  6:153407/90‑1 (MQ=255)
                                                                        aGTTTGCCAATCTTTCCATCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCaa                                                                                                                                       >  5:153407/1‑90 (MQ=255)
                                                                          tttGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAAc                                                                          <  2:62361/149‑1 (MQ=255)
                                                                             gCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAACCg                                                                        <  1:182189/148‑1 (MQ=255)
                                                                              ccAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAACCGcc                                                                      >  1:225836/1‑149 (MQ=255)
                                                                              ccAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAACCGc                                                                       <  2:333298/148‑1 (MQ=255)
                                                                               cAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGTATGATAAGAATAATCTGGAAGTGAACTTAACCGCCt                                                                     >  8:185926/1‑149 (MQ=255)
                                                                               cAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAACCGCCt                                                                     >  3:230077/1‑149 (MQ=255)
                                                                                    tttCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAACCGCCTGGcaa                                                                >  2:301521/1‑149 (MQ=255)
                                                                                        cGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAACCGCCTGGCaacaac                                                            >  7:296516/1‑149 (MQ=255)
                                                                                          tCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAACCGCCTGGCAACAACCt                                                          <  8:296516/149‑1 (MQ=255)
                                                                                              gTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAACCGCCTGGCAACAACCTTc                                                        <  8:269740/147‑1 (MQ=255)
                                                                                              gTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAACCGCCTGGCAACAACCTTc                                                        <  8:28684/147‑1 (MQ=255)
                                                                                              gTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAACCGCCTGGCAACAACCTTCAt                                                      >  7:156328/1‑149 (MQ=255)
                                                                                                 tCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAACCGCCTGGCAACAACCTTCATGAt                                                   >  1:82301/1‑149 (MQ=255)
                                                                                                     gACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAACCGCCTGGCAACAACCTTCATGATtgt                                                <  6:131097/148‑1 (MQ=255)
                                                                                                        gATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGCTCGTACCATGACGCTGCAGGCAACTATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAACCGCCTGGCAACAACCTTCATGATTGTGAtt                                            >  6:302776/1‑149 (MQ=255)
                                                                                                            aacaacTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAACCGCCt                                                                     >  4:283179/1‑120 (MQ=255)
                                                                                                            aacaacTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAACCGCCt                                                                     <  3:283179/120‑1 (MQ=255)
                                                                                                             acaacTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAACCGCCTGGCAACAACCTTCATGATTGTGATTTccc                                        >  1:134563/1‑148 (MQ=255)
                                                                                                               aacTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGTAGGCAACGATTCa                                                                                                                                >  2:276110/1‑58 (MQ=255)
                                                                                                               aacTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGTAGGCAACGATTCa                                                                                                                                <  1:276110/58‑1 (MQ=255)
                                                                                                               aacTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAACCGCCTGGCAACAACCTTc                                                        <  2:297348/130‑1 (MQ=255)
                                                                                                               aacTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAACCGCCTGGCAACAACCTTc                                                        >  1:297348/1‑130 (MQ=255)
                                                                                                               aacTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAACCGCCTGGCAACAACCTTCATGATTGTGATTTCCCGa                                      >  8:18924/1‑148 (MQ=255)
                                                                                                               aacTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAACCGCCTGGCAACAACCTTCATGATTGTGATTTCCCGa                                      >  6:264296/1‑148 (MQ=255)
                                                                                                                 cTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCCGCAGGCAACGATTCAGCAT‑GCGTGATCTGGTTATAGGATGATAAGAATAATCTGGAAGTGAACATAACCGCCTGGCAACAACCTTCATGATTGTGATTTCCCGACATg                                  >  8:286174/1‑149 (MQ=255)
                                                                                                                    gCGAAAGTAAGCGCCTTTGCCGATCGTACCATGCCGCTGCAGGCAACGATTCAGGATTGCGTGACCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAACCGCCTGGCAACAACCTTCATGATTGTGATTTCCCGACATgat                                <  4:225562/149‑3 (MQ=255)
                                                                                                                    gCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAACCGCCTGGCAACAACCTTCATGATTGTGATTTCCCGACATGtt                                >  3:32781/1‑149 (MQ=255)
                                                                                                                                     tGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAACCGCCTGGCAACAACCTTCATGATTGTGATTTCCCGACATGTTGCTATTCCCGATGGt                 <  1:301521/147‑1 (MQ=255)
                                                                                                                                       ccGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAACCGCCTGGCAACAACCTTCATGATTGTGATTTCCCGACATGTTGCTATTCCCGATGGTGAGc             >  2:69717/1‑149 (MQ=255)
                                                                                                                                       ccGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGGGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAACCGCCTGGCAACAACCTTCATGATTGTGATTTCCCGACATGTTGCTATTCCCGATGGTGAGc             <  2:225836/149‑1 (MQ=255)
                                                                                                                                        cGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAACCGCCTGGCAACAACCTTCATGATTGTGATTTCCCGACATGTTGCTATTCCCGATGGTGAGCt            >  4:233880/1‑149 (MQ=255)
                                                                                                                                        cGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAACCGCCTGGCAACAACCTTCATGATTGTGATTTCCCGACATGTTGCTATTCCCGATGGTGAGCt            >  8:31687/1‑149 (MQ=255)
                                                                                                                                           tCGTACCCTGACGCTGCAGGCAACGATTCAGGATGGCGTGAGCTGGGTATCGGATGATAAGAATAATCTGGAAGTGAACTTAACCGCCTGGCAACAACCTTCATGATTGTGAATTCCCGACATGTTGCTATTCCCGATGGTGAGCTTga         <  4:66620/149‑1 (MQ=255)
                                                                                                                                           tCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAACCGCCTGGCAACAACCTTCATGATTGTGATTTCCCGACATGTTGCTATTCCCGATGGTGAGCTTga         <  7:18924/149‑1 (MQ=255)
                                                                                                                                            cGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAACCGCCAGGCAACAACCTTCATGATTGTGATTTCCCGACATGTTGCTATTCCCGATGGTGAGCTTgag        <  8:243455/149‑1 (MQ=255)
                                                                                                                                              tACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAACCGCCTGGCAACACCCTTCATGATTGTGATTTCCCGACATGTTGCTATTCCCGATGGTGAGCtt           >  5:257305/1‑144 (MQ=255)
                                                                                                                                              tACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAACCGCCTGGCAACAACCTTCATGATTGTGATTTCCCGACATGTTGCTATTCCCGATGGTGAGCtt           >  4:35754/1‑144 (MQ=255)
                                                                                                                                                 cATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAACCGCCTGGCAACAACCTTCATGATTGTGATTTCCCGACATGTTGCTATTCCCGATGGTGAGCTTGAGATCAc   >  4:25257/1‑149 (MQ=255)
                                                                                                                                                  aTGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGGGAACTTAACCGCCTGGCAACAACCTTCATGATTGTGATTTCCCGACATGTTGCTATTCCCGATGGTGAGCTTGAGATCAcc  <  8:273585/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AAGTGGCGCTGTTTACCTATAATAGTCGGGCGGCGCAAATCTGGTGGCAGCAAAATCAGAGCAAATGTGTGCAGTTTGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTGGTTATCGGATGATAAGAATAATCTGGAAGTGAACTTAACCGCCTGGCAACAACCTTCATGATTGTGATTTCCCGACATGTTGCTATTCCCGATGGTGAGCTTGAGATCACC  >  NC_000913/214592‑214886

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