Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,768,131 A→T 100% pseudogene (911/945 nt) rhsJ → Rhs‑family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,768,1310AT100.0% 714.5 / NA 199pseudogene (911/945 nt)rhsJRhs‑family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (105/94);  total (105/94)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 199 total aligned reads displayed.

ACTAATGGTAAGTTAAAAATGGATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCAAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTC  >  NC_000913/3767989‑3768274
                                                                                                                                              |                                                                                                                                               
aCTAATGGTAAGTTAAAAATGGATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTAc                                                                                                                                           >  3:92185/1‑149 (MQ=255)
    aTGGTAAGTTAAAAATGGATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGc                                                                                                                                         >  6:85294/1‑147 (MQ=255)
        tAAGTTAAAAATGGATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGaa                                                                                                                                    >  8:186332/1‑148 (MQ=255)
          aGTTAAAAATGGATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTCGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGATAGGTACGCCAGaaaaa                                                                                                                                 >  7:91352/1‑149 (MQ=255)
              aaaaaTGGATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGTATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAAccc                                                                                                                              <  5:168472/148‑1 (MQ=255)
               aaaaTGGATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCa                                                                                                                             <  3:335611/148‑1 (MQ=255)
                aaaTGGATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGt                                                                                                                           >  4:136699/1‑149 (MQ=255)
                aaaTGGATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGt                                                                                                                           <  6:147662/149‑1 (MQ=255)
                  aTGGATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAg                                                                                                                                               >  5:212574/1‑127 (MQ=255)
                   tGGATACAGTGCTAACAGCACCAACGAACCTTAACTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCa                                                                                                                             <  5:145052/144‑1 (MQ=255)
                    ggATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGata                                                                                                                       >  7:193098/1‑149 (MQ=255)
                      aTACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGatat                                                                                                                      >  4:52078/1‑148 (MQ=255)
                       tACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATaa                                                                                                                    <  8:103036/149‑1 (MQ=255)
                       tACAGTGCTAACAGCACCAAAGAACCTTAATTCATATATGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGTTAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATaa                                                                                                                    <  6:173890/149‑1 (MQ=255)
                           gTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGCGCACCCCATTTATAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGa                                                                                                                  >  2:195218/1‑147 (MQ=255)
                                 ccAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTg                                                                                                           >  7:272171/2‑148 (MQ=255)
                                 aCAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTg                                                                                                           >  6:300257/1‑148 (MQ=255)
                                 aCAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTg                                                                                                           >  2:271375/1‑148 (MQ=255)
                                  cAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTg                                                                                                           >  8:50730/1‑147 (MQ=255)
                                    gCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGa                                                                                                       <  4:345606/149‑1 (MQ=255)
                                       ccAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATg                                                                                                     >  2:342198/1‑148 (MQ=255)
                                         aaCGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTAt                                                                                                  <  5:341955/149‑1 (MQ=255)
                                         aaCGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGAACTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTAt                                                                                                  <  1:160937/149‑1 (MQ=255)
                                            gAACCTTAATTCAAATAAGGATCTACTGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGata                                                                                                                       >  8:137961/1‑125 (MQ=255)
                                            gAACCTTAATTCAAATAAGGATCTACTGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGata                                                                                                                       <  7:137961/125‑1 (MQ=255)
                                             aaCCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGaaa                                                                                              <  5:300257/149‑1 (MQ=255)
                                              aCCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGaaaa                                                                                                                                  >  2:61767/1‑112 (MQ=255)
                                              aCCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGaaaa                                                                                                                                  <  1:61767/112‑1 (MQ=255)
                                              aCCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCa                                                                                                                             <  6:21000/117‑1 (MQ=255)
                                                cTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCa                                                                                           >  1:130666/1‑149 (MQ=255)
                                                 ttAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAACCCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCAt                                                                                          >  5:231983/1‑149 (MQ=255)
                                                 ttAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATaa                                                                                                                    >  5:162137/1‑123 (MQ=255)
                                                 ttAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCAt                                                                                          >  8:131249/1‑149 (MQ=255)
                                                   aaTTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCAt                                                                                          <  3:136699/147‑1 (MQ=255)
                                                    aTTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCt                                                                                       >  6:35674/1‑149 (MQ=255)
                                                     ttCAAATAAGGATCTAATGGCTATAACAACAAAAGCCCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCtt                                                                                      >  2:172371/1‑149 (MQ=255)
                                                     ttCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAg                                                                                                        >  8:207794/1‑131 (MQ=255)
                                                        aaaTAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGtt                                                                                   <  4:261483/149‑1 (MQ=255)
                                                              ggATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTAtt                                                                             >  3:174337/1‑149 (MQ=255)
                                                               gATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTAttt                                                                            <  1:271375/149‑1 (MQ=255)
                                                                aTCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTAtttt                                                                           >  1:265430/1‑149 (MQ=255)
                                                                aTCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTCATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTAtttt                                                                           >  6:54845/1‑149 (MQ=255)
                                                                    aaTGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTa                                                                       >  6:76243/1‑149 (MQ=255)
                                                                     aTGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTAttttt                                                                          <  5:117274/145‑1 (MQ=255)
                                                                     aTGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTa                                                                       <  2:130666/148‑1 (MQ=255)
                                                                     aTGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTa                                                                       <  7:185710/148‑1 (MQ=255)
                                                                     aTGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTAt                                                                      <  7:131249/149‑1 (MQ=255)
                                                                     aTGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTAt                                                                      <  5:85294/149‑1 (MQ=255)
                                                                     aTGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCAAGCTTGTTACTATTTTTTTAt                                                                      >  7:154355/1‑149 (MQ=255)
                                                                      tGGCTATAACAACAAAAGACCTTGGTAATTGCATAACATATGGGCCATCTGGGCACCCCATTTCTAATAAGCTAGGAACGCCAGAAAAACCAAGTGAAATAAGATAGGTTGTAGATGTATCAAAGAATGCTAGTTACTATTTTTTTAtt                                                                     >  6:199796/1‑149 (MQ=255)
                                                                        gCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGAt                                                                                                      <  1:383293/114‑1 (MQ=255)
                                                                        gCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTAttttt                                                                          <  8:84648/142‑1 (MQ=255)
                                                                         cTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTtct                                                                  >  2:144873/1‑149 (MQ=255)
                                                                           ataACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTa                                                                                                         >  3:206186/1‑108 (MQ=255)
                                                                           ataACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTa                                                                                                         <  4:206186/108‑1 (MQ=255)
                                                                           ataACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTa                                                                       >  1:39438/1‑142 (MQ=255)
                                                                               caacaaAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAgg                                                                                                              >  7:343563/1‑99 (MQ=255)
                                                                               caacaaAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAgg                                                                                                              <  8:343563/99‑1 (MQ=255)
                                                                                aacaaAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGt                                                                                                                           <  2:387623/85‑1 (MQ=255)
                                                                                aacaaAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATTCTTGTTACTATTTTTTTATTTCTCTTTAAt                                                           >  6:253597/1‑149 (MQ=255)
                                                                                aacaaAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTAt                                                           >  6:369018/1‑149 (MQ=255)
                                                                                aacaaAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTAt                                                           <  7:234924/149‑1 (MQ=255)
                                                                                  caaAAGACCTTGGTAATTGGTTAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTAttt                                                         <  8:251767/149‑1 (MQ=255)
                                                                                  caaAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTAttt                                                         <  6:107431/149‑1 (MQ=255)
                                                                                   aaaaGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGaaa                                                                                              >  5:66138/1‑111 (MQ=255)
                                                                                   aaaaGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTg                                                        >  7:58055/1‑149 (MQ=255)
                                                                                     aaGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGt                                                       >  7:293718/1‑148 (MQ=255)
                                                                                         ccTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAAtt                                                  <  8:234606/149‑1 (MQ=255)
                                                                                          cTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTg                                                 >  3:128074/1‑149 (MQ=255)
                                                                                              gTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATg                                             <  5:363512/149‑1 (MQ=255)
                                                                                                 aTTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCtt                                                               <  6:315201/128‑1 (MQ=255)
                                                                                                 aTTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCtt                                                               >  5:315201/1‑128 (MQ=255)
                                                                                                  ttGGATAACATATGGGCCATCTGGGCCCCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTTTGAAAGCATGCTTGTTACTATTTTTTTATTTCTGTTTTATTTGTCAATTGCATGGTGc                                         >  6:119337/1‑149 (MQ=255)
                                                                                                     gATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCCGTGAGATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTTCTATTTTTTTATTTCTCTGTTATTTGTCAATTGCCTGGTGCTGa                                      >  2:315763/1‑149 (MQ=255)
                                                                                                     gATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGt                                                                                                                           <  3:217882/64‑1 (MQ=255)
                                                                                                     gATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGt                                                                                                                           >  4:217882/1‑64 (MQ=255)
                                                                                                     gATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGa                                      >  1:187198/1‑149 (MQ=255)
                                                                                                     gATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCTGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGAt                                     >  1:101536/1‑149 (MQ=255)
                                                                                                       tAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGTTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGAtg                                    >  6:181567/1‑149 (MQ=255)
                                                                                                          cATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTtc                                                                   >  1:291288/1‑115 (MQ=255)
                                                                                                          cATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCa                                                     <  3:251346/129‑1 (MQ=255)
                                                                                                          cATATGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCa                                                     >  4:251346/1‑129 (MQ=255)
                                                                                                            tatGGGCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGaa                                >  7:100636/1‑148 (MQ=255)
                                                                                                                ggCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATaa                           <  4:128074/149‑1 (MQ=255)
                                                                                                                 gCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCtttt                                                             >  4:95775/1‑114 (MQ=255)
                                                                                                                 gCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATaaa                          <  3:241403/149‑1 (MQ=255)
                                                                                                                 gCCATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATaaa                          <  8:154355/149‑1 (MQ=255)
                                                                                                                  ccATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAATCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATaaaa                         <  5:41056/149‑1 (MQ=255)
                                                                                                                  ccATCTGGGCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTAttttttt                                                                        <  2:377147/102‑1 (MQ=255)
                                                                                                                         ggCACCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATaaaacaaaac                   <  1:144873/148‑1 (MQ=255)
                                                                                                                            aCCCCATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGaaa               <  5:205351/149‑1 (MQ=255)
                                                                                                                             ccccATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAg              >  2:225341/1‑149 (MQ=255)
                                                                                                                              cccATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATTCTTGTTACTATTTTTTTATTTCTATTTTATTTGTCAATTGCATGGTGATGATGTGAAAATAAAACAAAACGAAAGt             >  8:274225/1‑149 (MQ=255)
                                                                                                                              cccATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGAtgt                                   >  8:177304/1‑127 (MQ=255)
                                                                                                                              cccATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATaaaac                        <  7:211466/138‑1 (MQ=255)
                                                                                                                              cccATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAg              >  4:180518/1‑148 (MQ=255)
                                                                                                                              cccATTTCTAATATGCTAGGTACGCCAGAAAAACACAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAg              >  4:103616/1‑148 (MQ=255)
                                                                                                                                cATTTCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTat           <  8:293718/149‑1 (MQ=255)
                                                                                                                                    tCTAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTatgatg       >  2:454/1‑149 (MQ=255)
                                                                                                                                      tAATATGCTAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATggg     >  2:333498/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                               
ACTAATGGTAAGTTAAAAATGGATACAGTGCTAACAGCACCAACGAACCTTAATTCAAATAAGGATCTAATGGCTATAACAACAAAAGACCTTGGTAATTGGATAACATATGGGCCATCTGGGCACCCCATTTCTAATATGCAAGGTACGCCAGAAAAACCCAGTGATATAAGATAGGTTGTAGATGTATGAAAGCATGCTTGTTACTATTTTTTTATTTCTCTTTTATTTGTCAATTGCATGGTGCTGATGTGAAAATAAAACAAAACGAAAGTATGATGGGGTC  >  NC_000913/3767989‑3768274

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: