Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,198,900 G→A 6.0% G37G (GGG→GGA nfi → endonuclease V; deoxyinosine 3' endonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,198,9000GA 6.0% 214.5 / 2.7 67G37G (GGG→GGAnfiendonuclease V; deoxyinosine 3' endonuclease
Reads supporting (aligned to +/- strand):  ref base G (36/27);  new base A (2/2);  total (38/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 = 9.11e-01

ATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGC  >  NC_000913/4198752‑4199048
                                                                                                                                                     |                                                                                                                                                    
aTCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTAAGCGCTCAACAAATCGAACTGGCTTCTTATGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCggg                                                                                                                                                      >  3:5708/1‑149 (MQ=255)
   gACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGCTCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGttt                                                                                                                                                   <  1:173048/149‑1 (MQ=255)
    aCTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTg                                                                                                                                                  <  2:356773/149‑1 (MQ=255)
     cTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGttt                                                                                                                                                   >  7:323259/1‑147 (MQ=255)
     cTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGa                                                                                                                                                 <  5:307827/149‑1 (MQ=255)
     cTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGa                                                                                                                                                 <  8:323259/149‑1 (MQ=255)
        tCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAg                                                                                                                                                >  1:130173/1‑147 (MQ=255)
        tCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGATTTGAg                                                                                                                                                >  7:275981/1‑147 (MQ=255)
        tCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGATTTGAg                                                                                                                                                >  7:281163/1‑147 (MQ=255)
          tGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAgg                                                                                                                                            >  7:376499/1‑149 (MQ=255)
          tGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAgg                                                                                                                                            >  2:47739/1‑149 (MQ=255)
                     aCCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGt                                                                                                                                  >  5:329586/1‑148 (MQ=255)
                     aCCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGt                                                                                                                                  >  8:306058/1‑148 (MQ=255)
                     aCCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTg                                                                                                                                 <  4:117874/149‑1 (MQ=255)
                            aGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGATTTGAGCAGGGCGGAGAAGTGAcgcg                                                                                                                            <  8:275981/147‑1 (MQ=255)
                            aGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGATTTGAGCAGGGCGGAGAAGTGAcgcg                                                                                                                            <  8:281163/147‑1 (MQ=255)
                             ggAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAgc                                                                                                                         >  1:289550/1‑149 (MQ=255)
                                   aTTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGAt                                                                                                                   <  8:150091/149‑1 (MQ=255)
                                    ttATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGaa                                                                                                                                    <  8:109146/131‑1 (MQ=255)
                                    ttATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGaa                                                                                                                                    >  7:109146/1‑131 (MQ=255)
                                    ttATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATg                                                                                                                  <  2:289550/149‑1 (MQ=255)
                                       tGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGttt                                                                                                                                                   <  2:311306/113‑1 (MQ=255)
                                       tGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGttt                                                                                                                                                   >  1:311306/1‑113 (MQ=255)
                                       tGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGaa                                                                                                                                    >  2:180543/1‑128 (MQ=255)
                                       tGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGaa                                                                                                                                    <  1:180543/128‑1 (MQ=255)
                                        ggATCTCGCGTCATTACGCGCTCACAAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGtgc                                                                                                              >  5:28516/1‑149 (MQ=255)
                                        ggATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGa                                                                                                                                                 <  4:303690/114‑1 (MQ=255)
                                        ggATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGa                                                                                                                                                 >  3:303690/1‑114 (MQ=255)
                                          aTCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGtgctg                                                                                                            >  7:275073/1‑149 (MQ=255)
                                             tcGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGa                                                                                                                                                 >  8:103714/1‑109 (MQ=255)
                                             tcGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGa                                                                                                                                                 <  7:103714/109‑1 (MQ=255)
                                             tcGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTgctgct                                                                                                          >  6:24257/1‑148 (MQ=255)
                                                cgTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATg                                                                                                                  >  5:263057/1‑137 (MQ=255)
                                                cgTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATg                                                                                                                  <  6:263057/137‑1 (MQ=255)
                                                       aCGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATccc                                                                                                <  7:34008/148‑1 (MQ=255)
                                                        cgcgCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCt                                                                                               <  7:306058/148‑1 (MQ=255)
                                                               aacaaATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAg                                                                                       >  1:235392/1‑149 (MQ=255)
                                                               aacaaATCGAACTGGCTTCTTCTGTGATCAGCGAGGATCGACTCGATAAAGATACACCGGATCTGATCTCCGGAGCCGATGTCGGGTTTGAGCAGGGCGCAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGc                                                                                            >  5:89990/1‑144 (MQ=255)
                                                                acaaATTCGAACTGGCTTGTTCTGTGCTGCGCGAGGATCGACTCGATCTAGTTCGACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGGTGAAATATCCCTCGc                                                                                            <  6:89990/144‑1 (MQ=255)
                                                                 caaATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCt                                                                                     >  7:51517/1‑149 (MQ=255)
                                                                        gAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGa                                                                               <  6:28516/149‑1 (MQ=255)
                                                                          aCTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGt                                                                             <  6:329586/149‑1 (MQ=255)
                                                                           cTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCAACCTGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGt                                                                             <  4:5708/148‑1 (MQ=255)
                                                                            tGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGta                                                                            >  4:59716/1‑148 (MQ=255)
                                                                              gCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGCAATATCCCTCGCTTGAGCTGGTc                                                                                 >  7:180151/1‑141 (MQ=255)
                                                                              gCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGCAATATCCCTCGCTTGAGCTGGTc                                                                                 <  8:180151/141‑1 (MQ=255)
                                                                              gCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGtata                                                                          <  5:24257/148‑1 (MQ=255)
                                                                                     ctGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGtt                                                                     <  8:275073/146‑1 (MQ=255)
                                                                                     ctGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGcc                                                                  >  5:92687/1‑149 (MQ=255)
                                                                                        tgATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATaaa                                                                        >  4:328362/1‑140 (MQ=255)
                                                                                        tgATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATaaa                                                                        <  3:328362/140‑1 (MQ=255)
                                                                                             cgcgAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGc                                                          >  5:136102/1‑149 (MQ=255)
                                                                                              gcgAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGcc                                                         >  2:130790/1‑149 (MQ=255)
                                                                                                  ggATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGccacca                                                     >  2:321565/1‑149 (MQ=255)
                                                                                                                aaGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACAtt                                       >  7:122097/1‑149 (MQ=255)
                                                                                                                 aGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTc                                      >  8:124310/1‑149 (MQ=255)
                                                                                                                  gATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTc                                      >  6:209267/1‑148 (MQ=255)
                                                                                                                   aTCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCa                                    <  3:59716/149‑1 (MQ=255)
                                                                                                                   aTCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCa                                    <  8:122097/149‑1 (MQ=255)
                                                                                                                    tCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCt                                                                                               >  8:68831/1‑89 (MQ=255)
                                                                                                                     ccACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAgg                                  >  2:127981/1‑149 (MQ=255)
                                                                                                                      cACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGt                                                                             <  7:18634/105‑1 (MQ=255)
                                                                                                                      cACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGt                                                                             >  8:18634/1‑105 (MQ=255)
                                                                                                                      cACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGt                                 <  7:124310/149‑1 (MQ=255)
                                                                                                                       aCCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGtt                                >  1:349554/1‑149 (MQ=255)
                                                                                                                                gATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTttcct                       >  4:339845/1‑149 (MQ=255)
                                                                                                                                             gATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATcc          >  2:314709/1‑149 (MQ=255)
                                                                                                                                             gATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATcc          >  2:263757/1‑149 (MQ=255)
                                                                                                                                                   gggTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCgct    >  2:75616/1‑149 (MQ=255)
                                                                                                                                                     gTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCgctgc  <  5:209267/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAGTATAAAGTTGCCCGCATCGCCACCACCATGCCTTACATTCCAGGTTTTCTTTCCTTCCGCGAATATCCTGCGCTGC  >  NC_000913/4198752‑4199048

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: