Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 178,535 A→T 7.0% F207I (TTC→ATC)  mtn ← 5'‑methylthioadenosine/S‑adenosylhomocysteine nucleosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913178,5350AT 7.0% 217.5 / 5.0 70F207I (TTC→ATC) mtn5'‑methylthioadenosine/S‑adenosylhomocysteine nucleosidase
Reads supporting (aligned to +/- strand):  ref base A (31/34);  new base T (3/2);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGTTTCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACGGGACGTTGAAATTGTGGCAGACATGGGCGATTGCCGTCGCTTCCATCTCTACAGCAATGGCCTGTGGGAAGTTGTGGCGGATTTTCGC  >  NC_000913/178388‑178679
                                                                                                                                                   |                                                                                                                                                
gcgcggcgcgGCGTTGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGTTTCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGaa                                                                                                                                                  <  8:1882/148‑1 (MQ=255)
   cggcgcggcgTTGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGTTTCTGCACCAGTGCCTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGc                                                                                                                                                >  6:31820/1‑147 (MQ=255)
       gcggcgTTGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGTTTCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGa                                                                                                                                           <  3:230394/148‑1 (MQ=255)
              tGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGTTTCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAg                                                                                                                                            >  1:285648/1‑140 (MQ=255)
              tGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGTTTCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAg                                                                                                                                            <  2:285648/140‑1 (MQ=255)
                aGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCTATGTGCAAGTTTCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGt                                                                                                                                  <  4:187426/148‑1 (MQ=255)
                aGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGTTTCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTg                                                                                                                                   <  2:224165/147‑1 (MQ=255)
                   cacaGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGTTTCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGACCTCATCGAAGCTAAGATGAGACTGTTGa                                                                                                                               <  2:48244/148‑1 (MQ=255)
                   cacaGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGTTTCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGa                                                                                                                               >  5:32330/1‑148 (MQ=255)
                   cacaGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGTTTCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGa                                                                                                                               >  8:196084/1‑148 (MQ=255)
                     caGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGTTTCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATc                                                                                                                             >  8:273766/1‑148 (MQ=255)
                       gTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATATGCAAGTTTCTGCACCAGTGACTCAACCATCAGGATGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCgg                                                                                                                           >  1:320284/1‑148 (MQ=255)
                            gcAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGTTTCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGa                                                                                                                               >  3:136493/1‑139 (MQ=255)
                            gcAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGTTTCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGa                                                                                                                               <  4:136493/139‑1 (MQ=255)
                             cAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGTTTCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGATGCTAAGATGAGACTGTTGATCGGCCACGt                                                                                                                     >  5:120904/1‑148 (MQ=255)
                                  aaaGCCAGCGCGACCAGCGCCCTGAACAGTGCCTTAGCCATGTGCAAGTTTCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCCTCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGgag                                                                                                                <  6:248599/148‑1 (MQ=255)
                                  aaaGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGTTTCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGgag                                                                                                                <  7:199399/148‑1 (MQ=255)
                                  aaaGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGTTTCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGgag                                                                                                                <  8:98095/148‑1 (MQ=255)
                                    aGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGTTTCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGAt                                                                                                              <  3:93850/148‑1 (MQ=255)
                                           gCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGTTTCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGt                                                                                                                                  >  6:70702/1‑121 (MQ=255)
                                           gCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGTTTCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGt                                                                                                                                  <  5:70702/121‑1 (MQ=255)
                                                  gcgcCCTGAACAGTGACTTAGCCATGTGCAAGTTTCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGacaa                                                                                                >  3:199404/1‑148 (MQ=255)
                                                   cgcCCTGAACAGTGACTTAGCCATGTGCAAGTTTCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGacaa                                                                                                >  1:141779/1‑147 (MQ=255)
                                                       cTGAACAGTGACTTAGCCATGTGCAAGTTACTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACAAATCGAAGCTAAGATGAGACTGTTCATCGGCCACGTCCGAGATGGCGCGTACGACAACAAAc                                                                                           >  6:167025/1‑148 (MQ=255)
                                                          aaCAGTGACTTAGCCATGTGCAAGTTTCTGCACTAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACgg                                                                                         >  3:235684/1‑147 (MQ=255)
                                                          aaCAGTGACTTAGCCATGTGCAAGTTTCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACg                                                                                          <  2:141779/146‑1 (MQ=255)
                                                            cAGTGACTTAGCCATGTGCAAGTTTCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGt                                                                                                       >  8:2960/1‑131 (MQ=255)
                                                            cAGTGACTTAGCCATGTGCAAGTTTCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGt                                                                                                       <  7:2960/131‑1 (MQ=255)
                                                                  cTTAGCCATGTGCAAGTTTCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGATGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACGGGACGTTGa                                                                                 >  4:58486/1‑147 (MQ=255)
                                                                  cTTAGCCATGTGCAAGTTACTGCACCAGTGACTCAACCATCAGGCTCGACTGTTTAGCGGCAACAGCCAGGAACTCATCGCTGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACGGGACGTTGa                                                                                 >  4:62493/1‑147 (MQ=255)
                                                                    tAGCCATGTGCAAGTTTCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACGGGACGTTGa                                                                                 <  8:15437/145‑1 (MQ=255)
                                                                    tAGCCATGTGCAAGTTTCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACGGGACGTTGa                                                                                 >  7:15437/1‑145 (MQ=255)
                                                                      gCCATGTGCAAGTTTCTGCACTAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACGGGACGTTGAAATtg                                                                            <  4:235684/148‑1 (MQ=255)
                                                                         aTGTGCAAGTTTCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACGGGACGTTGAAATTGTgg                                                                         <  6:247284/148‑1 (MQ=255)
                                                                          tgtgCAAGTTTCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACGGGACGTTGAAATTGTGGc                                                                        <  6:32330/148‑1 (MQ=255)
                                                                                aGTTTCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACGGGACGTTGaa                                                                                <  2:161333/134‑1 (MQ=255)
                                                                                aGTTTCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACGGGACGTTGaa                                                                                >  1:161333/1‑134 (MQ=255)
                                                                                aGTTTCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGAAACAGCCAGGAACTCATCGAAGCTATGATGAGCCTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACGGGACGTTGAAATTGTGGCAGACAt                                                                  >  2:221231/1‑148 (MQ=255)
                                                                                  tttCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACggg                                                                                        <  8:21639/124‑1 (MQ=255)
                                                                                  tttCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACggg                                                                                        >  7:21639/1‑124 (MQ=255)
                                                                                   ttCTGCACCAGTGACTCAACCATCAGGCTGGACGGTTTAGCGGCAACAGCCAGGAACTCAACGAATCTAAGATGAAACTCTTTATCGGCCACGTCGGAGATCGCGCGTACTACAAAAAACCGGACCTTCAAATATTGGCAGACAAggg                                                               >  8:268131/1‑148 (MQ=255)
                                                                                        cACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACGGGACGTTGAAATTGTGGc                                                                        <  8:3634/134‑1 (MQ=255)
                                                                                        cACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACGGGACGTTGAAATTGTGGc                                                                        >  7:3634/1‑134 (MQ=255)
                                                                                        cACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACGGGACGTTGAAATTGTGGCAGACATGGGCGAtt                                                          >  8:6816/1‑148 (MQ=255)
                                                                                        cACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACGGGACGTTGAAATTGTGGCAGACATGGGCGAt                                                           <  5:31820/147‑1 (MQ=255)
                                                                                          ccAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACGGGACGTTGAAATTGTGGCAGACATGGGCGATTGc                                                        <  6:153975/148‑1 (MQ=255)
                                                                                              tGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACGGGACGTTGAAATTGTGGCAGACATGGGCGATTGCCGTc                                                    >  8:96949/1‑148 (MQ=255)
                                                                                                 cTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGAGGGCGCGTACGACAACAAACGGGACGTTGAAATTGTGGCAGACATGGGCGATTGCCGTCGCt                                                 <  8:57681/148‑1 (MQ=255)
                                                                                                   cAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGAcaacaa                                                                                             <  6:38554/102‑1 (MQ=255)
                                                                                                   cAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGAcaacaa                                                                                             >  5:38554/1‑102 (MQ=255)
                                                                                                   cAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACGGGACGTTGAAATTGTGGCAGACATGGGCGATTGCCGTCGCTTc                                               >  6:123814/1‑148 (MQ=255)
                                                                                                           aGGCTGGACTGTTTAGCGGCAACAGCAAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACGGGACGTTGAAATTGTGGCAGACATGGGCGATTGCCGTCGCTTCCATCTCTa                                       <  6:120904/148‑1 (MQ=255)
                                                                                                              cTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACGGGACGTTGAAATTGTGGCAGACATGGGCGATTGCCg                                                      >  8:25187/1‑130 (MQ=255)
                                                                                                              cTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACGGGACGTTGAAATTGTGGCAGACATGGGCGATTGCCg                                                      <  7:25187/130‑1 (MQ=255)
                                                                                                                 gACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACGGGACGTTGAAATTGTGGCAGACATGGGCGATTGCCGTCGCTTCCATCTCTACAGCaa                                 <  7:96949/148‑1 (MQ=255)
                                                                                                                  aCTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACGGGACGTTGAAATTGTGGCAGACATGGGCGATTGCCGTCGCTTCCATCTCTACAGCAAt                                >  2:96783/1‑148 (MQ=255)
                                                                                                                   cTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACGGGACGTTGAAATTGTGGCAGACATGGGCGATTGCCGTCGCTTTCATCTCTACAGCAATg                               >  8:4795/1‑148 (MQ=255)
                                                                                                                     gTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGCCCACGTCGGAGATGGCGCGTACGCCAACAAACGGGAAGTTGAAATTTTGGCAGACATGGCCGATTGCCGTCGCTACCATCTCTACAGCAATGGc                             >  6:70822/1‑148 (MQ=255)
                                                                                                                      tttAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACGGGACGTTGAAATTGTGGCAGACATGGGCGATTGCCGTCGCTTCCATCTCTACAGCAATgg                              <  7:273766/146‑1 (MQ=255)
                                                                                                                          gCGGCAACAGCCAGGAACTCATCGATGCTAAGATGAGACTGTTGATCGGCGACGGCGGAGATGGGGCGTACGACATCAAACGGGACGTTGAAATTGTGGCAGACATGGGCGATTGCCGTCGCTTCCATCTCTACAGCAATGGCCtgtg                        <  3:62493/148‑1 (MQ=255)
                                                                                                                          gCGGCAACAGCCAGGAACTCATCGATGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACGGGACGTTGAAATTGTGGCAGACATGGGCGATTGCCGTCGCTTCCATCTCTACAGCAATGGCCtgtg                        <  3:58486/148‑1 (MQ=255)
                                                                                                                                 cAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACGGGACGTTGTAATTGTGGCAGACATGGGCGATTGCCGt                                                     <  3:4508/112‑1 (MQ=255)
                                                                                                                                 cAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACGGGACGTTGAAATTGTGGCAGACATGGGCGATTGCCGt                                                     >  4:4508/1‑112 (MQ=255)
                                                                                                                                      aGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACGGGACGTTGAAATTGTGGCAGACATGGGCGATTGCCGTCGCTTCCATCTCTACAGCAATGGCCTGTGGGAAGTGGTGGc            >  6:1022/1‑148 (MQ=255)
                                                                                                                                       ggAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACGGGACGTTGAAATTGTGGCAGACATGGGCGATTGCCGTCGCTTCCATCTCTACAGCAATGGCCTGTGGGAAGTTGTGGCg           <  1:96783/148‑1 (MQ=255)
                                                                                                                                           cTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACGGGACGTTGAAATTGTGGCAGACATGGGCGATTGCCGTCGCTTc                                               >  6:183458/1‑108 (MQ=255)
                                                                                                                                           cTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACGGGACGTTGAAATTGTGGCAGACATGGGCGATTGCCGTCGCTTc                                               <  5:183458/108‑1 (MQ=255)
                                                                                                                                              aTCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACGGGACGTTGAAATTGTGGCAGACATGGGCGATTGCCGTCGCTTCCATCTCTACAGCAATGGCCTGTGGGAAGTTGTGGCGGAtttt     <  7:268131/147‑1 (MQ=255)
                                                                                                                                              aTCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACGGGACGTTGAAATTGTGGCAGACATGGGCGATTGCCGTCGCTTCCATCTCTACAGCAATGGCCTGTGGGAAGTTGTGGCGGAtttt     <  7:196084/147‑1 (MQ=255)
                                                                                                                                               tCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACGGGACGTTGAAATTGTGGCAGACATGGGCGATTGCCGTCGCTTCCATCTCTACAGCAATGGCCTGTGGGAAGTTGTGGCGGATTTTCg   <  5:1022/148‑1 (MQ=255)
                                                                                                                                                cGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACGGGACGTTGAAATTGTGGCAGACATGGGCGATTGCCGTCGCTTCCATCTCTACAGCAATGGCCTGTGGGAAGTTGTGGCGGATTTTCGc  >  3:176717/1‑148 (MQ=255)
                                                                                                                                                   aGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACg                                                                                          <  2:15202/57‑1 (MQ=255)
                                                                                                                                                   aGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACg                                                                                          >  1:15202/1‑57 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGTTTCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTGTTGATCGGCCACGTCGGAGATGGCGCGTACGACAACAAACGGGACGTTGAAATTGTGGCAGACATGGGCGATTGCCGTCGCTTCCATCTCTACAGCAATGGCCTGTGGGAAGTTGTGGCGGATTTTCGC  >  NC_000913/178388‑178679

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