Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 545,486 A→C 76.3% I293S (ATT→AGT)  allD ← ureidoglycolate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913545,4860AC76.3% 272.1 / 87.6 142I293S (ATT→AGT) allDureidoglycolate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (19/15);  new base C (47/61);  total (66/76)
Fisher's exact test for biased strand distribution p-value = 2.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.66e-01

Only 100 of 143 total aligned reads displayed.

CATCAACAATTTCGATGCCTTCGACGGCGGCTTTGCGTTGTTTAATATCCTGATCCTGTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTAATGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTGGAGGAGAAAAAGTTCGGATTAATAACTATATGTAATTGCCCCA  >  NC_000913/545388‑545582
                                                                                                  |                                                                                                
catcaACAATTTCGATGCCTTCGACGGCGGCTTTGCGTTGTTTAATATCCTGATCCTGTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTAAt                                                                                                 <  1:2273682/100‑1 (MQ=255)
 atcaACAATTTCGATGCCTTCGACGGCGGCTTTGCGTTGTTTAATATCCTGATCCTGTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTACTg                                                                                                >  2:668719/1‑100 (MQ=255)
 atcaACAATTTCGATGCCTTCGACGGCGGCTTTGCGTTGTTTAATATCCTGATCCTGTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTACTg                                                                                                <  1:1785729/100‑1 (MQ=255)
  tcaACAATTTCGATGCCTTCGACGGCGGCTTTGCGTTGTTTAATATCCTGATCCTGTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTACTgg                                                                                               <  2:3091483/100‑1 (MQ=255)
   caacaaTTTCGATGCCTTCGACGGCGGCTTTGCGTTGTTTAATATCCTGATCCTGTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTAATGGc                                                                                              >  2:3851717/1‑100 (MQ=255)
    aacaaTTTCGATGCCTTCGACGGCGGCTTTGCGTTGTTTAATATCCTGATCCTGTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTACTGGCa                                                                                             >  1:2627497/1‑100 (MQ=255)
    aacaaTTTCGATGCCTTCGACGGCGGCTTTGCGTTGTTTAATATCCTGATCCTGTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTACTGGCa                                                                                             <  2:3856056/100‑1 (MQ=255)
    aacaaTTTCGATGCCTTCGACGGCGGCTTTGCGTTGTTTAAGATCCTGATCCTGTCCGGGACAATAAACCTGATTAAAACCGGGCGCGGGGGTACTGGCa                                                                                             >  2:2272085/1‑100 (MQ=255)
      caaTTTCGATGCCTTCGACGGGGGCTTTGCGTTGTTTAATTTCCTGATCCTGGCCGGGATAAGAAACCTGATTAAAACCGGGCGCGGGGGTACTGGCAtt                                                                                           <  2:3783429/100‑1 (MQ=255)
      caaTTTCGATGCCTTCGACGGCGGCTTTGCGTTGTTTAATATCCTGATCCTGTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTAATGGCAtt                                                                                           >  1:91354/1‑100 (MQ=255)
         tttCGATGCCTTCGACGGCGGCTTTGCGTTGTTTAATATCCTGATCCTGTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTaa                                                                                        <  2:1912809/100‑1 (MQ=255)
             gATGCCTTCGACGGCGGCTTTGCGTTGTTTAATATCCTGATCCTGTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTcg                                                                                    >  1:1870038/1‑100 (MQ=255)
               tGCCTTCGACGGCGGCTTTGCGTTGTTTAATATCCTGATCCTGTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTcgcg                                                                                  <  1:1419192/100‑1 (MQ=255)
               tGCCTTCGACGGCGGCTTTGCGTTGTTTAATATCCTGATCCTGTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTAATGGCATTTAATTcgcg                                                                                  <  1:2501106/100‑1 (MQ=255)
                  cTTCGACGGCGGCTTTGCGTTGTTTAATATCCTGATCCTATCCGGGATATTAAACCTGATTAAAACCGCGCGTGGTGGTAATGGCATATAATTCGCGGAt                                                                               >  2:675231/1‑100 (MQ=255)
                   ttCGACGGCGGCTTTGCGTTGTTTGATATCCTGATCCTGTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATg                                                                              <  1:1022836/100‑1 (MQ=255)
                   ttCAACGGCGGCTTTGCGTTGTTTAATATCCTGATCCTGTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTAATGGCATTTACTTCGCGCa                                                                                <  2:3884869/98‑1 (MQ=255)
                    tCGACGGCGGCTTTGCGTTGTTTAATATCCTGATCCTGTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTAATGGCATTTAATTCGCGCATgg                                                                             <  1:854199/100‑1 (MQ=255)
                      gACGGCGGCTTTGCGTTGTTTAATATCCTGATCCTGTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTAATGGCATTTAATTCGCGCATGGTc                                                                           <  1:1142316/100‑1 (MQ=255)
                          gcggcTTTGCGTTGTTTAATATCCTGATCCTGTCCGGGATAATAACCCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGc                                                                       <  1:828484/100‑1 (MQ=255)
                          gcggcTTTGCGTTGTTTAATATCCTGATCCTGTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGc                                                                       <  2:3032272/100‑1 (MQ=255)
                          gcggcTTTGCGTTGTTTAATATCCTGATCCTGTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGc                                                                       >  1:2437345/1‑100 (MQ=255)
                            ggcTTTGCGTTGTTTAATATCCTGATCCTGTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTa                                                                     <  2:322389/100‑1 (MQ=255)
                            ggcTTTGCGTTGTTTAATATCCTGATCCTGTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTa                                                                     <  1:1356547/100‑1 (MQ=255)
                             gcTTTGCGTTGTTTAATATCCTGATCCTGTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTa                                                                     <  2:1952631/99‑1 (MQ=255)
                                 tGCGTTGTTTAATATCCTGATCCTGTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATg                                                                >  1:2567083/1‑100 (MQ=255)
                                   cGTTGTTTAATATCCTGATCCTGTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATACGCGAATGGTCTGGCTAAGATGtt                                                              >  2:3587439/1‑100 (MQ=255)
                                    gttgttTAATATCCTGATCCTGTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTg                                                             <  1:992170/100‑1 (MQ=255)
                                    gttgttTAATATCCTGATCCTGTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTAATGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTg                                                             >  1:755658/1‑100 (MQ=255)
                                      tgttTAACATCCTGATCCTGTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGAc                                                           >  2:2178960/1‑100 (MQ=255)
                                       gttTAATATCCTGATCCTGTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGa                                                            >  1:1229954/1‑98 (MQ=255)
                                       gttTAATATCCTGATCCTGTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACg                                                          >  1:278162/1‑100 (MQ=255)
                                       gttTAATATCCTGATCCTGTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACg                                                          >  2:1327999/1‑100 (MQ=255)
                                        tttAATATCCTGATCCTGTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAGTTCGCGCATGGTCTGGCTAAGATGTTGACgg                                                         <  1:1430208/100‑1 (MQ=255)
                                        tttAATATCCTGATCCTGTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACgg                                                         >  2:4103647/1‑100 (MQ=255)
                                         ttAATATCCTGATCCTGTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGa                                                        >  2:666645/1‑100 (MQ=255)
                                         ttAATAACCTGATGCTGTCCGGGATAAGAAACCTGATTAAAACCGGGGGCGGGGGTAATGGCATTTAATTCGCGCATGGGCTGGCTAAGATGTTGACgg                                                         <  2:110003/99‑1 (MQ=255)
                                           aaTATCCTGATCCTGTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGaat                                                      >  2:3352099/1‑100 (MQ=255)
                                            atatCCTGATCCTGTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGaata                                                     >  2:2355970/1‑100 (MQ=255)
                                             tatCCTGATCCTGTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTAATGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGaataa                                                    >  1:2092339/1‑100 (MQ=255)
                                               tcctgatcctgTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAAtt                                                  <  1:3374349/100‑1 (MQ=255)
                                                cctgatcctgTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTc                                                 <  2:2171170/100‑1 (MQ=255)
                                                cctgagcctgTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAAtt                                                  <  2:1854904/99‑1 (MQ=255)
                                                 ctgatcctgTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCg                                                <  2:2022715/100‑1 (MQ=255)
                                                 ctgatcctgTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCg                                                <  1:3538211/100‑1 (MQ=255)
                                                  tgatcctgTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGc                                               <  1:1253388/100‑1 (MQ=255)
                                                  tgatcctgTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGc                                               >  1:148656/1‑100 (MQ=255)
                                                  tgatcctgTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGc                                               <  2:2822650/100‑1 (MQ=255)
                                                    atcctgTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTg                                             >  2:1575602/1‑100 (MQ=255)
                                                      cctgTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTAATGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTgga                                           <  1:2737664/100‑1 (MQ=255)
                                                      cctgTCAGTGATAATAAAGCTGATTCAATATGGTCGGGGCTGTGCTGGCATTAAATTCACGCATGGTCTCGCTAAGATGTTGTCGGAATAATTCGCTgga                                           <  2:929226/100‑1 (MQ=255)
                                                       ctgTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTggag                                          <  1:2303439/100‑1 (MQ=255)
                                                       ctgTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTAATGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTggag                                          <  1:338523/100‑1 (MQ=255)
                                                       ctgTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTAATGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTggag                                          >  1:247675/1‑100 (MQ=255)
                                                        tgTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTggagg                                         <  1:1331590/100‑1 (MQ=255)
                                                        tgTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGGCTGGCTAAGATGTTGACGGAATAATTCGCTggagg                                         <  2:1586355/100‑1 (MQ=255)
                                                           ccGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTGGAGgaga                                      <  2:3887261/100‑1 (MQ=255)
                                                                 taataaACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTGGAGGAGAAAAAGt                                >  1:1010382/1‑100 (MQ=255)
                                                                   ataaACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTGGAGGAGAAAAAGTTc                              >  2:1894713/1‑100 (MQ=255)
                                                                   ataaACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTGGAGGAGAAAAAGTTc                              <  1:2265332/100‑1 (MQ=255)
                                                                    taaACGTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTGGAGGAGAAAAAGTTCg                             <  1:545442/100‑1 (MQ=255)
                                                                    taaACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTGGAGGAGAAAAAGTTCg                             <  2:496233/100‑1 (MQ=255)
                                                                    taaACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTGGAGGAGAAAAAGTTCg                             <  1:3448494/100‑1 (MQ=255)
                                                                    taaACCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTGGAGGAGAAAAAGTTCg                             <  1:1951183/100‑1 (MQ=255)
                                                                     aaaCCTGCTTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGGTGTTGACGGAATAATTCGCTGGAGGAGAAAAAGTTCgg                            >  2:521189/1‑100 (MQ=255)
                                                                     aaaCCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTGGAGGAGAAAAAGTTCgg                            <  1:3625350/100‑1 (MQ=255)
                                                                      aaCCTGATTAAAGCCGGGCGCGGGGGTACTGGCATTTAATTCTCGCATGGACGGGCTAAGATGTTGACGGAATAATTCGCTGGAGGAGAAAAAGTTCGGa                           <  2:352132/100‑1 (MQ=255)
                                                                      aaCCTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTGGAGGAGAAAAAGTTCGGa                           <  1:3972074/100‑1 (MQ=255)
                                                                      aaCCTGATTAAAACCGGGCGCGGGGGTAATGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTGGAGGAGAAAAAGTTCGGa                           >  1:495747/1‑100 (MQ=255)
                                                                      aaCCTGATTAAAACCGGGCGCGGGGGTAATGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTGGAGGAGAAAAAGTTCGGa                           >  1:3730506/1‑100 (MQ=255)
                                                                       aCCTGATTAAAACCGGGCGCGGGGGTAATGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTGGAGGAGAAAAAGTTCGGAt                          <  1:3100250/100‑1 (MQ=255)
                                                                       aCCTGATTAAAACCGGGCGCGGGGGTAATGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTGGAGGAGAAAAAGTTCGGAt                          <  1:2765551/100‑1 (MQ=255)
                                                                        ccTGATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTGGAGGAGAAAAAGTTCGGAtt                         >  1:431681/1‑100 (MQ=255)
                                                                        ccTGATTAAAACCGGGCGCGGGGGTAATGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTGGAGGAGAAAAAGTTCGGAtt                         >  1:2634859/1‑100 (MQ=255)
                                                                           gATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTGGAGGAGAAAAAGTTCGGATtaat                      <  1:2723666/100‑1 (MQ=255)
                                                                           gATTAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTGGAGGAGAAAAAGTTCGGATtaat                      >  1:770461/1‑100 (MQ=255)
                                                                            aTTAACACCGGGCGCGGGGTTACTGTCATGGAATTCGCGCAGGGTCTGGAGAAGATGTTGACGGAATAATTCGCTGGAGGAGAAAAAGTTCGGAtt                         <  2:4153894/96‑1 (MQ=255)
                                                                            aTTAAAACCGGGCGCGGGGGTACTGGCATTTTATTCGCGCATGGGCTGGCTAAGATGTTGACGGAATAATTCGCTGGAGGAGAAAAAGTTCGGATtaata                     <  2:2811802/100‑1 (MQ=255)
                                                                             ttAAAACCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTGGAGGAGAAAAAGTTCGGATtaataa                    >  2:2938326/1‑100 (MQ=255)
                                                                               aaaaCCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGCTGTTGACGGAATAATTCGCTGGAGGAGAAAAAGTTCGGATTAATAACt                  >  2:3280067/1‑100 (MQ=255)
                                                                               aaaaCCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTGGAGGAGAAAAAGTTCGGATTAATAACt                  <  1:533121/100‑1 (MQ=255)
                                                                                aaaCCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTGGAGGAGAAAAAGTTCGGATTAATAACta                 >  2:1046755/1‑100 (MQ=255)
                                                                                aaaCCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTGGAGGAGAAAAAGTTCGGATTAATAACta                 >  2:690149/1‑100 (MQ=255)
                                                                                 aaCCGGGCGCGGGGGTAATGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTGGAGGAGAAAAAGTTCGGATTAATAACt                  <  2:1518967/98‑1 (MQ=255)
                                                                                  aCCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTGGAGGAGAAAAAGTTCGGATTAATAACtata               <  2:148654/100‑1 (MQ=255)
                                                                                  aCCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTGGAGGAGAAAAAGTTCGGATTAATAACtata               >  1:3890201/1‑100 (MQ=255)
                                                                                  aCCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTGGAGGAGAAAAAGTTCGGATTAATAACtata               >  2:788991/1‑100 (MQ=255)
                                                                                  aCCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTGGAGGAGAAAAAGTTCGGATTAATAACtata               >  1:2934325/1‑100 (MQ=255)
                                                                                  aCCGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTGGAGGAGAAAAAGTTCGGATTAATAACtata               <  1:263184/100‑1 (MQ=255)
                                                                                   ccGGGCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTGGAGGAGAAAAAGTTCGGATTAATAACTatat              >  1:3203875/1‑100 (MQ=255)
                                                                                    cGGGCGCGGGGATACTGGCATTTAATTCGCTCATTGTATGGCTAAGATGTTGACGGAATAATTCTCTGGAGGAGAAAAAGTTCCAATTACTAACTATATg             <  2:3538731/100‑1 (MQ=255)
                                                                                     gggCGGAGGGGTACTGGCATTTAATTCGAGAATGGTAAGGCTAAGATGTTGTCGTAATAATTAAAAAAAAAAAAAAAAGTTCGGATTAATAACTatat              <  2:430088/98‑1 (MQ=255)
                                                                                     gggCGCGGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTGGAGGAGAAAAAGTTCGGATTAATAACTATATGt            <  1:3758436/100‑1 (MQ=255)
                                                                                        cgcgGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTGGAGGAGAAAAAGTTCGGATTAATAACTATATGTAAt         <  2:350431/100‑1 (MQ=255)
                                                                                         gcgGGGGTTATGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAGTTCGCTGGAGGAGAAAAAGTTCGGGTTAATAACGATATGTAAt         <  2:2913868/99‑1 (MQ=255)
                                                                                         gcgGGGGTACTGGCATTTAATTCGCGCATGGTAGGGATAAGATGTTGGCTGAATAATTAAAAAAAAAAAAAAAAGTTCGGATTAATAACTATATGTAAt         <  2:2649348/99‑1 (MQ=255)
                                                                                          cgGGGGTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTGGAGGAGAAAAAGTTCGGATTAATAACTATATGTAATTg       <  1:566147/100‑1 (MQ=255)
                                                                                             gggTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTGGAGGAGAAAAAGTTCGGATTAATAACTATATGTAATTGccc    <  2:1490257/100‑1 (MQ=255)
                                                                                               gTACTGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTGGAGGAGAAAAAGTTCGGATTAATAACTATATGTAATTGCCCCa  <  2:1264250/100‑1 (MQ=255)
                                                                                                 aaTGGAATATAATTAGCGAAAGGACAAGCTAAGATGTTGACGTAATAATTAAAAAAAAAAAAAAAAGTTCGGATTAATAACTATATGTAATTGccc    <  2:1103021/96‑1 (MQ=255)
                                                                                                  |                                                                                                
CATCAACAATTTCGATGCCTTCGACGGCGGCTTTGCGTTGTTTAATATCCTGATCCTGTCCGGGATAATAAACCTGATTAAAACCGGGCGCGGGGGTAATGGCATTTAATTCGCGCATGGTCTGGCTAAGATGTTGACGGAATAATTCGCTGGAGGAGAAAAAGTTCGGATTAATAACTATATGTAATTGCCCCA  >  NC_000913/545388‑545582

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