Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,606,114 T→A 6.3% W122R (TGG→AGG)  yhhN → TMEM86 family putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,606,1140TA 6.3% 199.0 / 3.7 63W122R (TGG→AGG) yhhNTMEM86 family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (32/27);  new base A (2/2);  total (34/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01

CAACGTCTGATGTACGCCATCGGCGCGTTTTTCCTCTCGCACCTGCTGTACACCATCTATTTCGCCAGTCAGATGACGCTCTCTTTCTTCTGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGTTTATCGGCATGACGCTGGTGATGGTGTGGCTGGCAGGTGAACTGTGGTTCTTCCGTCCGACCGCTCCGGCGCTCTCTGCGTTTGTTGGCGCTTCGTTGCTGTTT  >  NC_000913/3605970‑3606257
                                                                                                                                                |                                                                                                                                               
cAACGTCTGATGTACGCCATCGGCGCGTTTTTCCTCTCGCACCTGCTGTACACCATCTATTTCGCCAGTCAGATGACGCTCTCTTTCTTCTGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGAc                                                                                                                                             >  3:226367/1‑149 (MQ=255)
 aaCGTCTGATGTACGCCATCGGCGCGTGTTTCCTCTCGCAACTGCTGTACACCATCTAGTTCGCCAGTCATATGACGCTCTCGTGCTTCTGGCCTCGACCACTGGTGCTGCTGGGTATGGGTGCGCTGTTACTGGCGATTATCTGGAcg                                                                                                                                            >  4:265462/1‑149 (MQ=255)
        gATGTACGCCATCGGCGCGTTTTTCCTCTCGCACCTGCTGTACACCATCTATTTCGCCAGTCAGATGACGCTCTCTTTCTTCTGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTg                                                                                                                                      >  7:249367/1‑148 (MQ=255)
        gATGTACGCCATCGGCGCGTTTTTCCTCTCGCACCTGCTGTACACCATCTATTTCGCCAGTCAGATGACGCTCTCTTTCTTCTGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTg                                                                                                                                      >  3:125537/1‑148 (MQ=255)
         aTGTACGCCATCGGCGCGTTTTTCCTCTCGCACCTGCTGTACACCATCTATTTCGCCAGTCAGATGACGCTCTCTTTCTTCTGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTg                                                                                                                                      <  4:154490/147‑1 (MQ=255)
         aTGTACGCCATCGGCGCGTTTTTCCTCTCGCACCTGCTGTACACCATCTATTTCGCCAGTCAGATGACGCTCTCTTTCTTCTCGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTgg                                                                                                                                     <  2:203322/148‑1 (MQ=255)
          tGTACGCCATCGGCGCGTTTTTCCTCTCGCACCTGCTGTACACCATCTATTTCGCCAGTCAGATGCCGCTCTCTTTCTTCTGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTTCTGGCGATTATCTGGACGCGCCTGGat                                                                                                                                   >  4:128675/1‑148 (MQ=255)
          tGTACGCCATCGGCGCGTTTTTCCTCTCGCACCTGCTGTACACCATCTATTTCGCCAGTCAGATGACGCTCTCTTTCTTCTGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGaa                                                                                                                                   <  5:363562/149‑1 (MQ=255)
            tACGCCATCGGCGCGTTTTTCCTCTCGCACCTGCTGAACACCATCTATTTCGCCAGTCAGATGACGCTCTCTTTCTTCTGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGaa                                                                                                                                   >  7:172299/1‑147 (MQ=255)
             aCGCCATCGGCGCGTTTTTCCTCTCGCACCTGCTGTACACCATCTATTTCGCCAGTCAGATGACGCTCTCTTTCTTCTGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGa                                                                                                                                    <  4:325888/145‑1 (MQ=255)
             aCGCCATCGGCGCGTTTTTCCTCTCGCACCTGCTGTACACCATCTATTTCGCCAGTCAGATGACGCTCTCTTTCTTCTGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGa                                                                                                                                    >  3:325888/1‑145 (MQ=255)
                 cATCGGCGCGTTTTTCCTCTCGCACCTGCTGTACACCATCTATTTCGCCAGTCAGATGACGCTCTCTTTCTTCTGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTAcc                                                                                                                            >  3:339569/1‑149 (MQ=255)
                     ggCGCGTTTTTCCTCTCGCACCTGCTGTACACCATCTATTTCGCCAGTCAGATGACGCTCTCTTTCTTCTGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGtt                                                                                                                         <  8:346543/148‑1 (MQ=255)
                     ggCGCGTTTTTCCTCTCGCACCTGCTGTACACCATCTATTTCGCCAGTCAGATGACGCTCTCTTTCTTCTGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTg                                                                                                                        <  4:339569/149‑1 (MQ=255)
                     ggCGCGTTTTTCCTCTCGCACCTGCTGTACACCATCTATTTCGCCAGTCAGATGACGCTCTCTTTCTTCTGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTg                                                                                                                        <  8:76946/149‑1 (MQ=255)
                       cgcgTTTTTCCTCTCGCACCTGCTGTACACCATCTATTTCGCCAGTCAGATGACGCTCTCTTTCTTCTGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTGGc                                                                                                                      <  7:196237/149‑1 (MQ=255)
                       cgcgTTTTTCCTCTCGCACCTGCTGTACACCATCTATTTCGCCAGTCAGATGACGCTCTCTTTCTTCTGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTGGc                                                                                                                      <  7:191204/149‑1 (MQ=255)
                           tttttCCTCTCGCACCTGCTGTACACCATCTATTTCGCCAGTCAGATGACGCTCTCTTTCTTCTGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTGGCCTAt                                                                                                                  <  5:3812/149‑1 (MQ=255)
                           gttttCCTCTCGCACCTGCTGTACACCATCTATTTCGCCAGTCAGATGACGCTCTCTTTCTTCTGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTGGCCTAt                                                                                                                  >  7:256092/2‑149 (MQ=255)
                             tttCCTCTCGCACCTGCTGTACACCATCTATTTCGCCAGTCAGATGACGCTCTCTTTCTTCTGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTGGCCTATCt                                                                                                                >  8:129252/1‑149 (MQ=255)
                                 ctctCGCACCTGCTGTACACCATCTATTTCGCCAGTCAGATGACGCTCTCTTTCTTCTGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCAc                                                                                                            <  8:256092/149‑1 (MQ=255)
                                      gCACCTGCTGTACACCATCTATTTCGCCAGTCAGATGACGCTCTCTTTCTTCTGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTTTCTGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGTTTa                                                                                                       >  3:218743/1‑149 (MQ=255)
                                      gCACCTGCTGTACACCATCTATTTCGCCAGTCAGATGACGCTCTCTTTCTTCTGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGttt                                                                                                        <  6:177554/148‑1 (MQ=255)
                                      gCACCTGCTGTACACCATCTATTTCGCCAGTCAGATGACGCTCTCTTTCTTCTGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGTTTa                                                                                                       >  3:93523/1‑149 (MQ=255)
                                      gCACCTGCTGTACACCATCTATTTCGCCAGTCAGATGACGCTCTCTTTCTTCTGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGTTTa                                                                                                       <  4:93523/149‑1 (MQ=255)
                                             ctgTACACCATCTATTTCGCCAGTCAGATGACGCTCTATTTCTTCTGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGTTTATCGGc                                                                                                  <  4:226367/147‑1 (MQ=255)
                                                    ccATCTATTTCGCCAGTCAGATGACGCTCTCTTTCTTCTGGCCTCTACCACTGGTGCTGCTGGTTCTGGGGGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGTTTATCGGCATGACGCTg                                                                                         <  2:335856/149‑1 (MQ=255)
                                                     cATCTATTTCGCCAGTCAGATGACGCTCTCTTTCTTCTGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGTTTATCGGCATGACGCTgg                                                                                        <  7:129252/149‑1 (MQ=255)
                                                      aTCTATTTCGCCAGTCAGATGACGCTCTCTTTCTTCTGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGTTTATCGGCATGACGCTGGt                                                                                       >  1:105440/1‑149 (MQ=255)
                                                       tCTATTTCGCCAGTCAGATGACGCTCTCTTTCTTCTGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGTTTATCGGCATGACGCTGGt                                                                                       >  3:209358/1‑148 (MQ=255)
                                                         tATTTCGCCAGTCAGATGACGCTCTCTTTCTTCTGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTTCTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGTTTATCGGCATGACGCTGGTGAt                                                                                    >  2:308853/1‑149 (MQ=255)
                                                                ccAGTCAGATGACGCTCTCTTTCTTCTGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGTTTATCGGCATGACGCTGGTGATGGTGtgg                                                                             >  3:231399/1‑149 (MQ=255)
                                                                      aGATGACGCTCTCTTTCTTCTGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGTTTATCGGCATGACGCTGGt                                                                                       >  5:208919/1‑133 (MQ=255)
                                                                      aGATGACGCTCTCTTTCTTCTGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGTTTATCGGCATGACGCTGGt                                                                                       <  6:208919/133‑1 (MQ=255)
                                                                             gCTCTCTTTCTTCTGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCAGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGTTTATCGGCATGACGCTGGTGATGGTGTGGCTGGCAGGTGAAc                                                                >  7:215184/1‑149 (MQ=255)
                                                                             gCTCTCTTTCTTCTGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCAGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGTTTATCGGCATGACGCTGGTGATGGTGTGGCTGGCAGGTGAAc                                                                >  5:7717/1‑149 (MQ=255)
                                                                               tctctTTCTTCTGGCGTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGTTTATCGGCATGACGCTGGTGATGGTGTGGCTGGCAGGTGAACtg                                                              >  4:39558/1‑149 (MQ=255)
                                                                                 tctttcttCTGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGTTTATCGGCATGACGCTGGTGATGGTGTGGCTGGCAGGTGAACtgtg                                                            >  7:193848/1‑149 (MQ=255)
                                                                                   tttcttCTGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGTTTATCGGCATGACGCTGGTGATGGTGTGGCTGGCAGGTGAACTGTGGt                                                          >  7:322749/1‑149 (MQ=255)
                                                                                         ctGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGTTTATCGGCATGACGCTGGTGATGgt                                                                                 >  1:39949/1‑120 (MQ=255)
                                                                                         ctGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTAACGTTGGCCTATCTGCACGTTTATCGGCATGACGCTGGTGATGgt                                                                                 <  2:39949/120‑1 (MQ=255)
                                                                                          tGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACgcgc                                                                                                                                         <  8:220615/63‑1 (MQ=255)
                                                                                          tGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACgcgc                                                                                                                                         >  7:220615/1‑63 (MQ=255)
                                                                                               tctACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGTTTATCGGCATGACGCTGGTGATGGTGTGGCTGGCAGGTGAACTGTGGTTCTTCCGTCCGa                                              >  3:211832/1‑149 (MQ=255)
                                                                                                 tACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGTTTATCGGCATGACGCTGGTGATGGTGTGGCTGGCAGGTGAACTGTGGTTCTTCCGTCCGAcc                                            >  7:221847/1‑149 (MQ=255)
                                                                                                    cACTGGTGCTGCTGGTTCTTGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGTAAGAGTACCGTTGGCCTATCTGCACGTTTATCGGCATGACGCTGTTGATGGTGTGGCTGGCAGGTGAACTGTGGTTCTTCCGTCCGACGGc                                          <  8:172299/148‑1 (MQ=255)
                                                                                                      cTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGTTTATCGGCATGACGCTGGTGATGGTGTGGCTGGCAGGTGAACTGTGGTTCTTCCGTCCGACCGCTcc                                       <  8:221847/149‑1 (MQ=255)
                                                                                                            ctgctgGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGGTTATCGGCATGACGCTGGTGATGGTGTGGCTGGCAGGTGATCTGTGGTTCTTCCGTCCGACCGCTCCGGCGct                                 >  2:324214/1‑149 (MQ=255)
                                                                                                                  gTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGTTTATCGGCATGACGCTGGTGATGGTGTGGCTGGCAGGTGAACTGTGGTTCTTCCGTCCGACCGCTCCGGCGCTCTCTGc                           >  4:9110/1‑149 (MQ=255)
                                                                                                                    tCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGTTTATCGGCATGACGCTGGTGATGGTGTGGCTGGCAGGTGAACTGTGGTTCTTCCGTCCGACCGCTCCGGCGCTCTCTGCGt                         >  3:214056/1‑149 (MQ=255)
                                                                                                                      tGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGTTTATCGGCATGACGCTGGTGATGGTGTGGCTGGCAGGTGAACtg                                                              <  1:18409/110‑1 (MQ=255)
                                                                                                                      tGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGTTTATCGGCATGACGCTGGTGATGGTGTGGCTGGCAGGTGAACtg                                                              >  2:18409/1‑110 (MQ=255)
                                                                                                                          tGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGTTTATCGGCATGACGCTGGTGATGGTGTGGCTGGCAGGTGAACTGTGGTTCTTCCGTCCGACCGCTCCGGCGCTCTCTGCGTttgttg                   <  8:249367/149‑1 (MQ=255)
                                                                                                                          tGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGTTTATCGGCATGACGCTGGTGATGGTGTGGCTGGCAGGTGAACTGTGGTTCTTCCGTCCGACCGCTCCGGCGCTCTCTGCGTttgttg                   <  4:125537/149‑1 (MQ=255)
                                                                                                                          tGCGCTGTTACTGGCGATTATCAGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGTTTATCGGCATGACGCTGGTGATGGTGTGGCTGGCAGGTGAACTGTGGTTCTTCCGTCCGACCGCTCCGGCGCTCTCTGCGTttgttg                   <  6:7717/149‑1 (MQ=255)
                                                                                                                          tGCGCTGTTACTGGCGATTATCAGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGTTTATCGGCATGACGCTGGTGATGGTGTGGCTGGCAGGTGAACTGTGGTTCTTCCGTCCGACCGCTCCGGCGCTCTCTGCGTttgttg                   <  8:215184/149‑1 (MQ=255)
                                                                                                                           gcgcTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGTTTATCGGCATGACGCTGGTGATGGTGTGGCTGGCAGGTGAACTGTGGTTCTTCCGTCCGACCGCTCCGGCGCTCTCTGCGTTTGTTgg                  <  4:209358/149‑1 (MQ=255)
                                                                                                                              cTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGTTTATCg                                                                                                    <  6:124575/64‑1 (MQ=255)
                                                                                                                              cTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGTTTATCg                                                                                                    >  5:124575/1‑64 (MQ=255)
                                                                                                                                    cTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGTTTATCGGCATGACGCTGGTGATGGTGTGGCTGGCAGGTGAACTGTGGTTCTTCCGTCCGACCGCTCCGGCGCTCTCTGCGTTTGTTGGCGCTTCGtt         >  5:299507/1‑149 (MQ=255)
                                                                                                                                           ttATCTGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGTTTATCGGCATGACGCTGGTGATGGTGTGGCTGGCAGGTGAACTGTGGTTCTTCCGTCCGACCGCTCCGGCGCTCTCTGCGTTTGTTGGCGCTTCGTTGCTGtt   >  4:227436/1‑148 (MQ=255)
                                                                                                                                           ttATCTGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGTTTATCGGCATGACGCTGGTGATGGTGTGGCTGGCAGGTGAACTGTGGGTCTTCCGTCCGACCGCTCCGGCGCTCTCTGCGTTTGTTGGCGCTTCGTTGCTGttt  <  8:193848/149‑1 (MQ=255)
                                                                                                                                           ttATCTGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGTTTATCGGCATCACGCTGGTGATCGTGTGGCTCGCACGTGAACTGTGGTTCTTCCGTCCGACCGCTCCGGCGCTCTCTGCGTTTGTTGGCGCTTCGTTGCTGtt   >  6:164499/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                               
CAACGTCTGATGTACGCCATCGGCGCGTTTTTCCTCTCGCACCTGCTGTACACCATCTATTTCGCCAGTCAGATGACGCTCTCTTTCTTCTGGCCTCTACCACTGGTGCTGCTGGTTCTGGGTGCGCTGTTACTGGCGATTATCTGGACGCGCCTGGAAGAGTACCGTTGGCCTATCTGCACGTTTATCGGCATGACGCTGGTGATGGTGTGGCTGGCAGGTGAACTGTGGTTCTTCCGTCCGACCGCTCCGGCGCTCTCTGCGTTTGTTGGCGCTTCGTTGCTGTTT  >  NC_000913/3605970‑3606257

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