Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 988,789 A→C L66R (CTG→CGG)  asnS ← asparaginyl tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913988,7890AC99.5% 773.0 / ‑4.5 217L66R (CTG→CGG) asnSasparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (97/119);  total (98/119)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

Only 100 of 217 total aligned reads displayed.

TTCAACCTTGCTGGCCTGAATTTCAAATTGTTGCCCCTGGCCCGGCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCAGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAGCAGGAACCGTCATAAACGGCGAGGAAGGAGATGCCAG  >  NC_000913/988692‑988885
                                                                                                 |                                                                                                
    aCCTTGCTGGCCTGAATTTCAAATTGTTGCCCCTGGCCCGGCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGt                                                                                            >  2:2780685/1‑100 (MQ=255)
     ccTTGCTGGCCTGAATTTCAAATTGTTGCCCCTGGCCCGGCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTc                                                                                           >  2:1168337/1‑100 (MQ=255)
        tGCTGGCCTGAATTTCAAATTGTTGCCCCTGGCCCGGCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTc                                                                                        >  2:1785009/1‑100 (MQ=255)
            ggCCTGAATTTCAAATTGTTGCCCCTGGCCCGGCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTg                                                                                    >  2:2851847/1‑100 (MQ=255)
            ggCCTGAAGTTCAAATTGTTGCCCCTGGCCCGGCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTg                                                                                    <  1:729393/100‑1 (MQ=255)
              ccTGAATTTCAAATAGTTGCCCCTGGACCGGCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGGCGTCTTCGTTGTa                                                                                  >  2:5129625/1‑100 (MQ=255)
               cTGAATTTCAAATTGTTGCCCCTGGCCCGGCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTaa                                                                                 >  2:3993298/1‑100 (MQ=255)
                tGAATTTCAAATTGTTGCCCCTGGCCCGGCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAAt                                                                                <  1:3122876/100‑1 (MQ=255)
                tGAATTTCAAATTGTTGCCCCTGGCCCGGCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAAt                                                                                >  1:2050671/1‑100 (MQ=255)
                  aaTTTCAAATTGTTGCCCCTGGCCCGGCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTg                                                                              <  2:880396/100‑1 (MQ=255)
                  aaTTTCAAATTGTTGCCCCTGGCCCGGCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTg                                                                              <  2:4784894/100‑1 (MQ=255)
                    tttCAAATTGTTGCGCCTGGCCCGGCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTggg                                                                            <  2:4556930/100‑1 (MQ=255)
                    tttCAAATTGTTGCCCCTGGCCCGGCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTggg                                                                            >  1:4285415/1‑100 (MQ=255)
                     ttCAAATTGTTGCCCCTGGCCCGGCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGc                                                                           <  2:1017428/100‑1 (MQ=255)
                     ttCAAATTGTTGCCCCTGGCCCGGCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGc                                                                           <  1:5518974/100‑1 (MQ=255)
                     ttCAAATTGTTGCCCCTGGCCCGGCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGc                                                                           <  1:493725/100‑1 (MQ=255)
                     ttCAAATTGTTGCCCCTGGCCCGGCGACGCCACGACGTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGc                                                                           <  1:2711096/100‑1 (MQ=255)
                      tCAAATTGTTGCCCCTGGCCCGGCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGTCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCa                                                                          <  1:4862763/100‑1 (MQ=255)
                        aaaTTGTTGCCCCTGGCCCGGCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCaga                                                                        <  1:3113396/100‑1 (MQ=255)
                          aTTGTTGCCCCTGGCCCGGCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAgaga                                                                      <  1:4988644/100‑1 (MQ=255)
                               tgCCCCTGGCCCGGCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAattat                                                                 >  1:1666626/1‑100 (MQ=255)
                                gCCCCTGGCCCGGCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAattatt                                                                >  2:2786765/1‑100 (MQ=255)
                                 ccccTGGCCCGGCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATtg                                                               <  2:3133060/100‑1 (MQ=255)
                                 ccccTGGCCCGGCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATtg                                                               <  2:2718256/100‑1 (MQ=255)
                                 ccccTGGCCCGGCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATtg                                                               >  1:2158747/1‑100 (MQ=255)
                                  cccTGGCCCGGCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATtga                                                              >  1:3209418/1‑100 (MQ=255)
                                   ccTGGCCCGGCTACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCAGTGTAATTGGGCAGAGAATTATtgat                                                             <  2:4380954/100‑1 (MQ=255)
                                   ccTGGCCCGGCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATtgat                                                             >  1:2265556/1‑100 (MQ=255)
                                    ntGGCCCGGCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATtgatg                                                            >  1:4429379/2‑100 (MQ=255)
                                    cgggCCCGGCGACGCGACGACTTGACCCGTCACAATGAACGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGAAGAGAATTATtgatg                                                            <  2:1949780/98‑1 (MQ=255)
                                    cTGGCCCGGCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATtgatg                                                            <  1:4723595/100‑1 (MQ=255)
                                    cTGGCCCGGCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATtgatg                                                            >  1:635240/1‑100 (MQ=255)
                                    cTGGCCCGGCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATtgatg                                                            <  2:2499/100‑1 (MQ=255)
                                    cTGGCCCGGCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATtgatg                                                            <  1:4478749/100‑1 (MQ=255)
                                      ggCCCGGCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGAc                                                          <  2:3633216/100‑1 (MQ=255)
                                      ggCCCGGCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGAc                                                          <  2:3816479/100‑1 (MQ=255)
                                      ggCCCGGCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGAc                                                          <  2:853413/100‑1 (MQ=255)
                                       gCCCGGCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACa                                                         <  1:342158/100‑1 (MQ=255)
                                        cccGGCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAg                                                        >  1:1202621/1‑100 (MQ=255)
                                        cccGGCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAg                                                        <  1:939741/100‑1 (MQ=255)
                                            gCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTg                                                    <  2:1564517/100‑1 (MQ=255)
                                             cGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGt                                                   <  1:4810119/100‑1 (MQ=255)
                                             cGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGt                                                   <  1:299075/100‑1 (MQ=255)
                                              gACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGGATTATTGATGACAGCCTGTa                                                  >  1:500907/1‑100 (MQ=255)
                                                  ccACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAgg                                              >  2:2289694/1‑100 (MQ=255)
                                                  ccACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAgg                                              >  2:524240/1‑100 (MQ=255)
                                                  ccACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAgg                                              >  2:3512140/1‑100 (MQ=255)
                                                  ccACGACTTTACCCGTCACAATGACCGAGCAGCCGGAGGTCAGACGCCGGACCTCTTCATTGTAATTGGACAGAGAATTATTGATGACAGCCAGTACAgg                                              >  2:3069974/1‑100 (MQ=255)
                                                   cACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGa                                             >  2:5142260/1‑100 (MQ=255)
                                                    aCGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGAt                                            <  2:2408232/100‑1 (MQ=255)
                                                    aCGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGAt                                            >  2:3049120/1‑100 (MQ=255)
                                                    aCGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGACAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGAt                                            >  2:3048664/1‑100 (MQ=255)
                                                    aCGACTTTACCCGTCACAATGACCAAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGan                                            <  1:955423/100‑2 (MQ=255)
                                                     cGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATc                                           >  1:1001502/1‑100 (MQ=255)
                                                      gACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCa                                          <  2:1461233/100‑1 (MQ=255)
                                                        cTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCaaa                                        >  1:2847774/1‑100 (MQ=255)
                                                        cTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCaaa                                        <  1:2493087/100‑1 (MQ=255)
                                                         tttACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAg                                       <  1:1163184/100‑1 (MQ=255)
                                                         tctACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAg                                       <  2:2616965/98‑1 (MQ=255)
                                                           tACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAGCa                                     >  1:3592811/1‑100 (MQ=255)
                                                           tACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAGCa                                     <  2:3842898/100‑1 (MQ=255)
                                                            aCCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAGCAg                                    <  1:56226/100‑1 (MQ=255)
                                                            aCCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATGGTAATTGGGCAGATAATTATTGATGACAGCCTGTACAGGATCAAAGCAg                                    >  2:5415092/1‑100 (MQ=255)
                                                             cccGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAGCAgg                                   >  1:3395820/1‑100 (MQ=255)
                                                             cccGTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAGCAgg                                   <  1:3157727/100‑1 (MQ=255)
                                                                gTCACACTGACCGAGCAGACGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAGCAGGAAc                                >  1:4646882/1‑100 (MQ=255)
                                                                gTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAGCAGGAAc                                <  2:2242318/100‑1 (MQ=255)
                                                                gTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAGCAGGAAc                                <  2:4740991/100‑1 (MQ=255)
                                                                gTCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAGCAGGAAc                                >  1:5190261/1‑100 (MQ=255)
                                                                 tCACAATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAGCAGGAAcc                               >  1:3962505/1‑100 (MQ=255)
                                                                    caATGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAGCAGGAACCGTc                            >  2:3241229/1‑100 (MQ=255)
                                                                      aTGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAGCAGGAACCGTCAt                          <  2:3368008/100‑1 (MQ=255)
                                                                      aTGACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAGCAGGAACCGTCAt                          >  1:847228/1‑100 (MQ=255)
                                                                        gACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAGCAGGAACCGTCATaa                        >  1:2204310/1‑100 (MQ=255)
                                                                        gACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAGCAGGAACCGTCATaa                        >  2:580752/1‑100 (MQ=255)
                                                                        gACCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCATGTACAGGATCAAAGCAGGAACCGTCATaa                        >  2:528393/1‑100 (MQ=255)
                                                                         aCCGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAGCAGGAACCGTCATaaa                       >  1:3494668/1‑100 (MQ=255)
                                                                          ccGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAGCAGGAACCGTCATAAAc                      <  1:1413644/100‑1 (MQ=255)
                                                                           cGAGCAGCCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAGCAGGAACCGTCATAAACg                     <  1:1805475/100‑1 (MQ=255)
                                                                             agcagcCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAGCAGGAACCGTCATAAACGGc                   <  2:219888/100‑1 (MQ=255)
                                                                                agcCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAGCAGGAACCGTCATAAACGGCGag                <  2:627642/100‑1 (MQ=255)
                                                                                agcCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAGCAGGAACCGTCATAAACGGCGag                <  1:3077752/100‑1 (MQ=255)
                                                                                agcCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAGCAGGAACCGTCATAAACGGCGag                <  1:3077293/100‑1 (MQ=255)
                                                                                agcCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTCGTGACAGCCTGTACAGGATGAAAGCAGGAACCGTCATAAACGGCGag                <  2:3561470/100‑1 (MQ=255)
                                                                                 gcCGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAGCAGGAACCGTCATAAACGGCGagg               >  1:2252883/1‑100 (MQ=255)
                                                                                  ccGGTGGTCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAGCAGGAACCGTCATAAACGGCGagga              <  2:1096814/100‑1 (MQ=255)
                                                                                    ggtggtCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAGCAGGAACCGTCATAAACGGCGaggaag            >  1:2895764/1‑100 (MQ=255)
                                                                                     gtggtCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGGATTATTGATGACAGCCTGTACAGGATCAAAGCAGGAACCGTCATAAACGGCGaggaagg           <  2:1033792/100‑1 (MQ=255)
                                                                                     gctgtCATACGCCGGACGTCTTCATTGGAATGGGGCATAGAATTATTGATGACAGCCTGTACAGGATCAAAGCAGGAACCGTCATAAACGGCGaggaagg           <  2:4147622/97‑1 (MQ=255)
                                                                                        gtCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAGCAGGAACCGTCATAAACGGCGAGGAAGgaga        <  2:1348400/100‑1 (MQ=255)
                                                                                        gtCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAGCAGGAACCGTCATAAACGGCGAGGAAGgaga        >  1:4185367/1‑100 (MQ=255)
                                                                                         tCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAGCAGGAACCGTCATAAACGGCGAGGAAGGAGAt       <  1:2636216/100‑1 (MQ=255)
                                                                                         tCAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAGCAGGAACCGTCATAAACGGCGAGGAAGGAGAt       >  1:3275425/1‑100 (MQ=255)
                                                                                          cAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAGCAGGAACCGTCATAAACGGCGAGGAAGGAGATg      <  2:4453003/100‑1 (MQ=255)
                                                                                          cAGACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAGCAGGAACCGTCATAAACGGCGAGGAAGGAGATg      <  2:2845986/100‑1 (MQ=255)
                                                                                            gACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAGCAGGAACCGTCATAAACGGCGAGGAAGGAGATGcc    >  1:2869867/1‑100 (MQ=255)
                                                                                            gACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAGCAGGAACCGTCATAAACGGCGAGGAAGGAGATGcc    >  2:21383/1‑100 (MQ=255)
                                                                                            gACGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAGCAGGAACCGTCATAAACGGCGAGGAAGGAGATGcc    <  2:3466948/100‑1 (MQ=255)
                                                                                             aCGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAGCAGGAACCGTCATAAACGGCGAGGAAGGAGATGCCa   <  2:3179764/100‑1 (MQ=255)
                                                                                              cGCCGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAGCAGGAACCGTCATAAACGACGAGGAAGGAGATGCCAg  >  2:1879404/1‑100 (MQ=255)
                                                                                                 |                                                                                                
TTCAACCTTGCTGGCCTGAATTTCAAATTGTTGCCCCTGGCCCGGCGACGCCACGACTTTACCCGTCACAATGACCGAGCAGCCGGTGGTCAGACGCAGGACGTCTTCATTGTAATTGGGCAGAGAATTATTGATGACAGCCTGTACAGGATCAAAGCAGGAACCGTCATAAACGGCGAGGAAGGAGATGCCAG  >  NC_000913/988692‑988885

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