Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,490,088 A→G 54.2% G31G (GGT→GGC argD ← bifunctional acetylornithine aminotransferase and succinyldiaminopimelate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,490,0880AG54.2% 55.3 / 200.3 169G31G (GGT→GGCargDbifunctional acetylornithine aminotransferase and succinyldiaminopimelate aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (35/42);  new base G (50/41);  total (85/84)
Fisher's exact test for biased strand distribution p-value = 2.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

Only 100 of 172 total aligned reads displayed.

CACCAGCGCAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGA  >  NC_000913/3489989‑3490187
                                                                                                   |                                                                                                   
cACCAGCGCAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGa                                                                                                     >  2:2737571/1‑100 (MQ=255)
  ccAGCGCAGGCTGGGAATGGCCCAGCTCCGTAACTGCAATGCCACCCGCGAAATCGACAGAGTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGcc                                                                                                   <  1:3738950/100‑1 (MQ=37)
  ccAGCGCAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGcc                                                                                                   >  2:3944442/1‑100 (MQ=255)
   cAGCGCAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCt                                                                                                  >  1:2938590/1‑100 (MQ=255)
     gcgcAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCGGATTCGGCTGCCCTGGCCttt                                                                                                <  2:4070228/100‑1 (MQ=255)
        cAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTAc                                                                                             >  2:2418771/1‑100 (MQ=255)
         aGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTAcc                                                                                            <  2:1864465/100‑1 (MQ=255)
          ggacgacAATGGCCCAACGCCATAACTCCAATGCCCCCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCACTGGCCTTTTACCg                                                                                           <  1:2183310/94‑1 (MQ=255)
           gATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCgg                                                                                          <  2:1681390/100‑1 (MQ=255)
           gATGGCAATGGCCCAACGCCGAAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCgg                                                                                          <  1:2666455/100‑1 (MQ=255)
              ggCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAAt                                                                                       <  2:2661045/100‑1 (MQ=255)
              ggCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAAt                                                                                       >  2:1401620/1‑100 (MQ=255)
              ggCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAAt                                                                                       <  1:2691694/100‑1 (MQ=255)
               gCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATa                                                                                      <  2:4035538/100‑1 (MQ=255)
                  aTGGCCCAACGCCGTCACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAAc                                                                                   <  2:4213319/100‑1 (MQ=255)
                  aTGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAAc                                                                                   <  1:3237903/100‑1 (MQ=255)
                   ttgCCCAACGCCGTAACTGCCATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATCCGGCGGCCCTGACCTTTTACCGAAAAAAAc                                                                                   <  1:4264250/97‑1 (MQ=255)
                     gCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAActct                                                                                >  2:788152/1‑100 (MQ=255)
                     gCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAActct                                                                                >  2:577155/1‑100 (MQ=255)
                      cccAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCATATTCGCCTGCCCTGGCCTTTTACCGGAATAAACTCTg                                                                               >  1:2873781/1‑100 (MQ=255)
                      cccAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTg                                                                               >  2:2378472/1‑100 (MQ=255)
                       ccAACGCCGTAACTGCGATGCCACGCGCGAAATCGACAGATTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGc                                                                              <  1:3692569/100‑1 (MQ=255)
                       ccAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGc                                                                              <  1:591430/100‑1 (MQ=255)
                         aaCGCCGTAACTGCAATGCCACCCGCGAAATCGCCATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCg                                                                            <  2:1065556/100‑1 (MQ=255)
                            gCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAg                                                                         <  1:2401918/100‑1 (MQ=255)
                            gCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAg                                                                         >  2:2231840/1‑100 (MQ=255)
                            gCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAg                                                                         >  1:1123592/1‑100 (MQ=255)
                             ccGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGc                                                                        <  1:1368692/100‑1 (MQ=255)
                             ccGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGc                                                                        >  1:2240664/1‑100 (MQ=255)
                              cGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCa                                                                       <  2:3072510/100‑1 (MQ=255)
                                tAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATa                                                                     >  2:2038293/1‑100 (MQ=255)
                                     gCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATc                                                                 >  2:2002196/1‑99 (MQ=255)
                                      cAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCgg                                                               <  2:2785646/100‑1 (MQ=255)
                                      cAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCgg                                                               <  1:1624222/100‑1 (MQ=255)
                                      cAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCgg                                                               <  1:1623277/100‑1 (MQ=255)
                                        aTGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCa                                                             >  2:3703809/1‑100 (MQ=255)
                                          gCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAgg                                                           <  2:3115677/100‑1 (MQ=255)
                                           ccACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAGATCggaagag                                                          >  1:2521023/1‑95 (MQ=255)
                                            cACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGAt                                                         <  2:1432556/100‑1 (MQ=255)
                                            cACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGAt                                                         >  1:3410520/1‑100 (MQ=255)
                                              cccGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCa                                                       >  1:3904832/1‑100 (MQ=255)
                                              cccGCGAAATCGACATACTCCTTGCCTTGCTCATCCCACATTCGGCTGCCCTGACCTTTTCCCGGAATAAACTCTGCCGGCGCATAAATCGGCAGGATc                                                        <  1:2994407/99‑1 (MQ=255)
                                                cgcgAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACt                                                     >  2:359730/1‑100 (MQ=255)
                                                  cgAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTc                                                   <  2:2596493/100‑1 (MQ=255)
                                                  cgAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTc                                                   <  2:3793165/100‑1 (MQ=255)
                                                    aaaTCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCAt                                                 >  2:3608974/1‑100 (MQ=255)
                                                    aaaTCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCATGGCCTTTTACCGGAATAAACTCTGCCGGAGCATACATCGGCAGGAt                                                         >  1:4087018/1‑92 (MQ=255)
                                                      aTCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCg                                               <  1:529101/100‑1 (MQ=255)
                                                        cGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGaa                                             >  2:3367001/1‑100 (MQ=255)
                                                         gACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGaaa                                            >  1:1172471/1‑100 (MQ=255)
                                                         gACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGTCGGAGCATAAATCGGCAGGATCACTTCATCGaaa                                            <  2:3349879/100‑1 (MQ=255)
                                                           cATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGt                                          >  2:2654973/1‑100 (MQ=255)
                                                           cATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTGACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAg                                           <  2:657531/99‑1 (MQ=255)
                                                               cTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGCAAGTcg                                        >  1:2928497/1‑98 (MQ=255)
                                                               cTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTcgcg                                      <  2:3173745/100‑1 (MQ=255)
                                                               cTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTcgcg                                      >  2:1316162/1‑100 (MQ=255)
                                                                tCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTcgc                                       >  1:973369/1‑98 (MQ=255)
                                                                 ccttgccttgCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTcgcgcg                                    >  2:61790/1‑100 (MQ=255)
                                                                 ccttgccttgCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTcgcgcg                                    >  1:2769508/1‑100 (MQ=255)
                                                                 ccttgccttgCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTcgcgcg                                    <  2:2082875/100‑1 (MQ=255)
                                                                 ccttgccttgCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTcgcgcg                                    <  1:3898469/100‑1 (MQ=255)
                                                                    tgccttgcTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCgtgt                                 <  2:2674877/100‑1 (MQ=255)
                                                                     gccttgcTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTa                                >  2:215607/1‑100 (MQ=255)
                                                                       cttgcTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAAt                              <  2:631101/100‑1 (MQ=255)
                                                                       cttgcTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAAt                              <  2:1499562/100‑1 (MQ=255)
                                                                       cttgcTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAAt                              <  1:194764/100‑1 (MQ=255)
                                                                        ttgcTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAAtt                             >  1:3449073/1‑100 (MQ=255)
                                                                          gcTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGc                           >  1:4197138/1‑100 (MQ=255)
                                                                          gcTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGc                           <  1:3809521/100‑1 (MQ=255)
                                                                           cTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCt                          >  2:1456716/1‑100 (MQ=255)
                                                                            tGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCtg                         <  1:1460360/100‑1 (MQ=255)
                                                                            tGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCtg                         >  1:1729903/1‑100 (MQ=255)
                                                                            tGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCtg                         <  2:2979487/100‑1 (MQ=255)
                                                                            tGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCtg                         <  2:3717121/100‑1 (MQ=255)
                                                                            tGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCtg                         <  1:4220673/100‑1 (MQ=255)
                                                                            tGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCtg                         >  1:4225324/1‑100 (MQ=255)
                                                                              aTCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTCCATCGAAAGTCGCGCGTGTAATTGCtgtt                       >  1:2565232/1‑100 (MQ=255)
                                                                              aTCCCAGATTCGGCTGCCCTGGCCGTTTCCCGGAATAAACGCCGCCGGAGCATAAATCGGCAGGATCACCTCATCGAAAGTCGCGCGTGTAATTGCtgtt                       <  1:693631/100‑1 (MQ=255)
                                                                                nccAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCtgtttg                     >  1:3526721/2‑100 (MQ=255)
                                                                                cccAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCtgtttg                     <  2:3585563/100‑1 (MQ=255)
                                                                                 ccAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCtgtttgt                    >  2:3263064/1‑100 (MQ=255)
                                                                                   aGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGGGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTc                  <  2:2725216/100‑1 (MQ=255)
                                                                                   aGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTc                  <  2:1271047/100‑1 (MQ=255)
                                                                                    gATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCa                 <  1:313852/100‑1 (MQ=255)
                                                                                     aTTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCaa                <  2:3813638/100‑1 (MQ=255)
                                                                                        cGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTg             >  2:3053243/1‑100 (MQ=255)
                                                                                        cGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTg             <  1:3189933/100‑1 (MQ=255)
                                                                                        cGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTg             <  2:1169601/100‑1 (MQ=255)
                                                                                         ggCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGc            >  2:2936372/1‑100 (MQ=255)
                                                                                           cTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCa          <  2:2763772/100‑1 (MQ=255)
                                                                                           cTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCa          <  2:3335972/100‑1 (MQ=255)
                                                                                           cTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCa          >  2:186015/1‑100 (MQ=255)
                                                                                             gCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCAtc        >  1:1566563/1‑100 (MQ=255)
                                                                                             gCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCAtc        <  2:1124240/100‑1 (MQ=255)
                                                                                              cccTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCAtct       >  1:357065/1‑100 (MQ=255)
                                                                                              cccTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCAtct       >  1:1242353/1‑100 (MQ=255)
                                                                                              cccTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCAtct       <  1:1133423/100‑1 (MQ=255)
                                                                                              cccTGACCGTTTACCGGACTAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTACTTGCAGTTTGTTTATTTGACAtat       >  1:3188692/1‑98 (MQ=255)
                                                                                               ccTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCAtctc      >  2:138133/1‑100 (MQ=255)
                                                                                                   aCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGa  <  1:1961686/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CACCAGCGCAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCATGA  >  NC_000913/3489989‑3490187

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: