Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 234,579 C→A 6.1% K68N (AAG→AAT) ‡ gloB ← hydroxyacylglutathione hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913234,5790CA 6.1% 346.8 / 8.5 98K68N (AAG→AAT) ‡gloBhydroxyacylglutathione hydrolase
Reads supporting (aligned to +/- strand):  ref base C (45/47);  new base A (3/3);  total (48/50)
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.95e-01

CAGATATGTCCTAAAGTGTGACCCGGCGTAGCAATTACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATGGTGGGTGAGAAATATGGCCTCCGGTTGCCAGTTATTGGCGGCAATGGCGTTTAATACTGGCTCTGCGTCTCCGGGATCGACAATCAGGCAGCGACCT  >  NC_000913/234436‑234720
                                                                                                                                               |                                                                                                                                             
cAGATATGTCCTAAAGTGTGACCCGGCGTAGCAATTACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACtttt                                                                                                                                           <  1:842420/148‑1 (MQ=255)
         ccTAAAGTGTGACCCGGCGTAGCAATTACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGATCATACACCACAATTTGTGGAAACTTTTCCACCAGtt                                                                                                                                  >  3:758389/1‑148 (MQ=255)
         ccTAAAGTGTGACCCGGCGTAGCAATTACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGtt                                                                                                                                  >  4:1138120/1‑148 (MQ=255)
               gtgtGACCCGGCGTAGCAATTACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTa                                                                                                                            <  3:1138120/148‑1 (MQ=255)
                                   tACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCgtgg                                                                                                        <  4:758389/148‑1 (MQ=255)
                                        tAAGTTCATGCCCCAAACCGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCtt                                                                                                                               <  3:8417/120‑1 (MQ=255)
                                        tAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCtt                                                                                                                               >  4:8417/1‑120 (MQ=255)
                                        tAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTa                                                                                                                            >  4:505011/1‑123 (MQ=255)
                                        tAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTa                                                                                                                            <  3:505011/123‑1 (MQ=255)
                                         aaaTTCATGCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACtttt                                                                                                                                           >  2:785180/1‑106 (MQ=255)
                                         aaaTTCATGCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACtttt                                                                                                                                           <  1:785180/106‑1 (MQ=255)
                                                   cccAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACtttt                                                                                                                                           >  4:617469/1‑97 (MQ=255)
                                                   cccAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACtttg                                                                                                                                           <  3:617469/97‑2 (MQ=255)
                                                     cAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACtt                                                                                                                                             <  4:986476/93‑1 (MQ=255)
                                                     cAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACtt                                                                                                                                             >  3:986476/1‑93 (MQ=255)
                                                     cAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCaa                                                                                                                   <  4:280240/119‑1 (MQ=255)
                                                     cAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCaa                                                                                                                   >  3:280240/1‑119 (MQ=255)
                                                     cAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGAtgg                                                                                                  <  4:830242/136‑1 (MQ=255)
                                                     cAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGAtgg                                                                                                  >  3:830242/1‑136 (MQ=255)
                                                     cAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATGGTGGTTGAGAAat                                                                                      >  1:591956/1‑148 (MQ=255)
                                                     cAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATGGTGGGTGAGAAat                                                                                      >  3:551238/1‑148 (MQ=255)
                                                     cAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAAAGTTTccac                                                                                                                                       <  2:1017673/99‑1 (MQ=255)
                                                     cAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAAAGTTTccac                                                                                                                                       >  1:1017673/1‑99 (MQ=255)
                                                      aaaaCGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACtttt                                                                                                                                           >  4:288818/1‑94 (MQ=255)
                                                      aaaaCGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACtttt                                                                                                                                           <  3:288818/94‑1 (MQ=255)
                                                      aaaaCGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGt                                                                                                                <  4:584563/121‑1 (MQ=255)
                                                      aaaaCGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGt                                                                                                                >  3:584563/1‑121 (MQ=255)
                                                      aaaaCGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATGGTGGGTGAGAAata                                                                                     <  1:43807/148‑1 (MQ=255)
                                                        aaCGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATggtgggtg                                                                                            <  4:203444/139‑1 (MQ=255)
                                                        aaCGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATggtgggtg                                                                                            >  3:203444/1‑139 (MQ=255)
                                                        aaCGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATGGTGGGTGAGAAATATg                                                                                   <  2:591956/148‑1 (MQ=255)
                                                         aCGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTAcgcc                                                                                                                        <  2:668764/110‑1 (MQ=255)
                                                         aCGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTAcgcc                                                                                                                        >  1:668764/1‑110 (MQ=255)
                                                          cGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATggtgggtg                                                                                            <  3:920526/137‑1 (MQ=255)
                                                          cGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGGGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGAGggtgggtg                                                                                            >  4:920526/1‑137 (MQ=255)
                                                          cGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGCGTCTCTTGTGGACCATACACCACAATTTGTGGAAACtttt                                                                                                                                           >  1:961061/1‑90 (MQ=255)
                                                          cGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGCGTCTCTTGTGGACCATACACCACAATTTGTGGAAACtttt                                                                                                                                           <  2:961061/90‑1 (MQ=255)
                                                                   tttCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGAt                                                                                                    >  2:366797/1‑120 (MQ=255)
                                                                   tttCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGAt                                                                                                    <  1:366797/120‑1 (MQ=255)
                                                                     tCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATGGTGGGTGAGAAATATGGCCTCCGGTTGcc                                                                      <  4:551238/148‑1 (MQ=255)
                                                                        ccATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATGGTGGGTGAGAAATATGGCCTCCGGTTGCCAGt                                                                   >  3:573883/1‑148 (MQ=255)
                                                                            cTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATggtggg                                                                                              >  4:1143538/1‑117 (MQ=255)
                                                                            cTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATggtggg                                                                                              <  3:1143538/117‑1 (MQ=255)
                                                                                    aCCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATGGTGGGTGAGAAata                                                                                     >  1:29310/1‑118 (MQ=255)
                                                                                    aCCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATGGTGGGTGAGAAata                                                                                     <  2:29310/118‑1 (MQ=255)
                                                                                      cTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGAtg                                                                                                   <  3:181233/102‑1 (MQ=255)
                                                                                      cTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGAtg                                                                                                   >  4:181233/1‑102 (MQ=255)
                                                                                      cTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATGGTGGGTGAGAAata                                                                                     <  1:198851/116‑1 (MQ=255)
                                                                                      cTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATGGTGGGTGAGAAata                                                                                     >  2:198851/1‑116 (MQ=255)
                                                                                      cTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAAAGTTTcc                                                                                                                                         >  4:1123393/1‑64 (MQ=255)
                                                                                      cTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAAAGTTTcc                                                                                                                                         <  3:1123393/64‑1 (MQ=255)
                                                                                       tgtgTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGAtg                                                                                                   >  4:754413/1‑101 (MQ=255)
                                                                                       tgtgTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGAtg                                                                                                   <  3:754413/101‑1 (MQ=255)
                                                                                        gtgtTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCtttt                                                                                                                             <  4:52669/74‑1 (MQ=255)
                                                                                        gtgtTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCtttt                                                                                                                             >  3:52669/1‑74 (MQ=255)
                                                                                        gtgtTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATGGTGGGTGAGAAATATGGCCTCCGGTTGCCAGTTATTggcggc                                                         >  2:1082713/1‑142 (MQ=255)
                                                                                        gtgtTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATGGTGGGTGAGAAATATGGCCTCCGGTTGCCAGTTATTggcggc                                                         <  1:1082713/142‑1 (MQ=255)
                                                                                              ttCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCtt                                                                                                                               >  2:353837/1‑66 (MQ=255)
                                                                                              ttCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCtt                                                                                                                               <  1:353837/66‑1 (MQ=255)
                                                                                               tCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGAtggtg                                                                                                <  3:623839/96‑1 (MQ=255)
                                                                                               tCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGAtggtg                                                                                                >  4:623839/1‑96 (MQ=255)
                                                                                               tCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATGGTGGGTGAGAAATATGGCCTCCGGTTGCCAGTTATTGGCGGCAATGGCGTTTAAt                                            <  4:573883/148‑1 (MQ=255)
                                                                                                  cTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCtt                                                                                                                               <  3:1068946/62‑1 (MQ=255)
                                                                                                  cTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCtt                                                                                                                               >  4:1068946/1‑62 (MQ=255)
                                                                                                    tATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATGGTGGGTGAGAAATATGGCCTCCGGTTGCCAGTTATTGGCGGCAATg                                                     >  4:666211/1‑134 (MQ=255)
                                                                                                    tATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATGGTGGGTGAGAAATATGGCCTCCGGTTGCCAGTTATTGGCGGCAATg                                                     <  3:666211/134‑1 (MQ=255)
                                                                                                    tATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATGGTGGGTGAGAAATATGGCCTCCGGTTGCCAGTTATTGGCGGCAATGGCGTTTAATACTgg                                       >  3:66134/1‑148 (MQ=255)
                                                                                                            gtCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCtt                                                                                                                               >  3:792880/1‑52 (MQ=255)
                                                                                                            gtCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCtt                                                                                                                               <  4:792880/52‑1 (MQ=255)
                                                                                                             tctcTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCa                                                                                                                    >  1:118232/1‑62 (MQ=255)
                                                                                                             tctcTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCa                                                                                                                    <  2:118232/62‑1 (MQ=255)
                                                                                                              ctctTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATGGTGGGTGAGAAATATGGCCTCCGGt                                                                          >  1:555782/1‑103 (MQ=255)
                                                                                                              ctctTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATGGTGGGTGAGAAATATGGCCTCCGGt                                                                          <  2:555782/103‑1 (MQ=255)
                                                                                                              ctctTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATGGTGGGTGAGAAATATGGCCTCCGGTTGCCAGTTATTGGCGGCAATGGCGTTTAATACTGGCTCTGCGtct                             >  1:780215/1‑148 (MQ=255)
                                                                                                                ctTGTGGACCATACTCCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGAtggt                                                                                                 <  3:14120/78‑1 (MQ=255)
                                                                                                                ctTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCgtgg                                                                                                        >  4:642847/1‑71 (MQ=255)
                                                                                                                ctTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCgtgg                                                                                                        <  3:642847/71‑1 (MQ=255)
                                                                                                                ctTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGAtggt                                                                                                 >  4:14120/1‑78 (MQ=255)
                                                                                                                       aCCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGAtggtg                                                                                                <  3:566996/72‑1 (MQ=255)
                                                                                                                       aCCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGAtggtg                                                                                                >  4:566996/1‑72 (MQ=255)
                                                                                                                             caccacAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATGGTGGGTGAGAAATATGGCCTCCGGTTGCCAGTTATTGGCGGCAATGGCGTTTAATACTGGCTCTGCGTCTCCGGGATCGACAATc              <  4:66134/148‑1 (MQ=255)
                                                                                                                               ccacAATTTGTGGAAACTTTTCCACCAGTTCTTTTAcg                                                                                                                          >  4:481891/1‑38 (MQ=255)
                                                                                                                               ccacAATTTGTGGAAACTTTTCCACCAGTTCTTTTAcg                                                                                                                          <  3:481891/38‑1 (MQ=255)
                                                                                                                               ccacAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATggtgggtg                                                                                            <  4:306842/68‑1 (MQ=255)
                                                                                                                               ccacAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATggtgggtg                                                                                            >  3:306842/1‑68 (MQ=255)
                                                                                                                                cacaATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATGGTGGGTGAGAAATATGGCCTCCGGTTGCCAGTTATTGGCGGCAATGGCGTTTAATACTGGCTCTGCGTCTCCGGGGTCGACAATCAgg           >  1:614907/1‑148 (MQ=255)
                                                                                                                                 acaATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGAtggtg                                                                                                >  4:766781/1‑62 (MQ=255)
                                                                                                                                 acaATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGAtggtg                                                                                                <  3:766781/62‑1 (MQ=255)
                                                                                                                                 acaATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATGGTGGGTGAGAAATATGGCCTCCGGTTGCCAGTTATTGGCGGCAATGGCGTTTAATACTGGCTCTGCGTCTCCGGGATCGACAATCAGGc          <  2:780215/148‑1 (MQ=255)
                                                                                                                                   aaTTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATGGTGGGTGAGAAATATGGCCTCCGGTTGCCAGTTATTGGCGGCAATGGCGTTTAATACTGGCTCTGCGTCTCCGGGGTCGACAATAAGGCAg        <  2:614907/148‑1 (MQ=255)
                                                                                                                                         tgGAAAAGTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATGGTGGGTGAGAAATATGGCCTCCGGTTGCCAGTTATTgg                                                             >  3:974134/1‑89 (MQ=255)
                                                                                                                                         tgGAAAAGTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATGGTGGGTGAGAAACATGGCCTCCGGTTGCCAGTTGTTgg                                                             <  4:974134/89‑1 (MQ=255)
                                                                                                                                            aaaCTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATGGTGGGTGAGAAATATGGCCTCCGGTTGCCAGTTATTGGCGGCAATGGCGt                                                 >  4:1057249/1‑98 (MQ=255)
                                                                                                                                            aaaCTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATGGTGGGTGAGAAATATGGCCTCCGGTTGCCAGTTATTGGCGGCAATGGCGt                                                 <  3:1057249/98‑1 (MQ=255)
                                                                                                                                            aaaCTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATGGTGGGTGAGAAATATGGCCTCCGGTTGCCAGTTATTGGCGGCAATGGCGTTTAATACTGGCTCTGCGTCTCCGGGATCGACAATCAGGCAGCGACCt  >  3:931287/1‑145 (MQ=255)
                                                                                                                                            aaaCTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATGGTGGGTGAGAAATATGGCCTCCGGTTGCCAGTTATTGGCGGCAATGGCGTTTAATACTGGCTCTGCGTCTCCGGGATCGACAATCAGGCAGCGACCt  <  4:931287/145‑1 (MQ=255)
                                                                                                                                              aCTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATGGTGGGTGAGAAATATGGCCTCCGGTTGCCAGTTATTGGCGGCAATGGCGTTTAATACTGGCTCTg                                  <  3:209507/111‑1 (MQ=255)
                                                                                                                                              aCTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATGGTGGGTGAGAAATATGGCCTCCGGTTGCCAGTTATTGGCGGCAATGGCGTTTAATACTGGCTCTg                                  >  4:209507/1‑111 (MQ=255)
                                                                                                                                               |                                                                                                                                             
CAGATATGTCCTAAAGTGTGACCCGGCGTAGCAATTACACTAAATTCATGCCCCAAAACGAAGGCAGTTTCGCCATCTTTGACTACCTGTGTTGTTCCCTTATCTTGTGTCTCTTGTGGACCATACACCACAATTTGTGGAAACTTTTCCACCAGTTCTTTTACGCCGCCAACGTGATCGTGGTGATGGTGGGTGAGAAATATGGCCTCCGGTTGCCAGTTATTGGCGGCAATGGCGTTTAATACTGGCTCTGCGTCTCCGGGATCGACAATCAGGCAGCGACCT  >  NC_000913/234436‑234720

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