Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,490,088 A→G 15.7% 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,0880AG15.7% 379.4 / 56.1 174G31G (GGT→GGCargDbifunctional acetylornithine aminotransferase and succinyldiaminopimelate aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (71/73);  new base G (13/14);  total (86/88)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.31e-01

Only 100 of 179 total aligned reads displayed.

CACCAGCGCAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCAT  >  NC_000913/3489989‑3490185
                                                                                                   |                                                                                                 
cACCAGCGCAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGa                                                                                                   <  2:1329406/100‑1 (MQ=255)
cACCAGCGCAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGa                                                                                                   >  1:2791886/1‑100 (MQ=255)
cACCAGCGCAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGa                                                                                                   >  1:2410076/1‑100 (MQ=255)
 aCCAGCGCAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGAc                                                                                                  <  1:695355/100‑1 (MQ=255)
  ccAGCGCAGGATGGCAATGGCCCACCGTCGTAACTGCAATGCCACCCGCGAACTCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGAc                                                                                                  <  2:3446494/99‑1 (MQ=255)
  ccAGCGCAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGAcc                                                                                                 >  2:1043834/1‑100 (MQ=255)
  ccAGCGCAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGAcc                                                                                                 <  2:3317116/100‑1 (MQ=255)
    aGCGCAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCtt                                                                                               >  2:2964146/1‑100 (MQ=255)
    aGCGCAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCtt                                                                                               >  1:1271932/1‑100 (MQ=255)
    aGCGCAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCtt                                                                                               >  2:3008725/1‑100 (MQ=255)
        cAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGTCCTGACCtttttc                                                                                           <  1:2810927/100‑3 (MQ=255)
        cAGGATGACCATGGCCCAACGCCGTAACCGCACAGCAAATCGCCAACTCGACCTACTCCCTGCCTTGCTGCTCCCCGATTCGGCTGCCCTCACCATTTAc                                                                                           <  2:2921862/100‑1 (MQ=255)
         aGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTAcc                                                                                          <  1:2701765/100‑1 (MQ=255)
          ggATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCg                                                                                         >  1:443163/1‑100 (MQ=255)
            ntGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATCCTCCTTGCCATGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGa                                                                                       >  2:1784353/2‑100 (MQ=255)
            aTGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGa                                                                                       <  1:398717/100‑1 (MQ=255)
             tGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGaa                                                                                      <  1:2644497/100‑1 (MQ=255)
               gCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATa                                                                                    <  1:2145164/100‑1 (MQ=255)
               gCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATa                                                                                    <  1:3481835/100‑1 (MQ=255)
               gCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATa                                                                                    <  2:1871411/100‑1 (MQ=255)
               gCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATa                                                                                    <  1:133492/100‑1 (MQ=255)
                cAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATaa                                                                                   <  1:3658949/100‑1 (MQ=255)
                    ggCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAActc                                                                               >  1:3838280/1‑100 (MQ=255)
                    ggCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAActc                                                                               >  2:1242547/1‑100 (MQ=255)
                    ggCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAActc                                                                               <  2:1682096/100‑1 (MQ=255)
                     gCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAActct                                                                              >  1:304149/1‑100 (MQ=255)
                      cccTACGACGTAACTGCAATGTCACCCGGGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGTTGCCCTGTCCTTTTGCCGGAGCAAACTTTg                                                                             >  2:2190336/1‑100 (MQ=255)
                      cccAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTg                                                                             <  2:2791870/100‑1 (MQ=255)
                      cccAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTg                                                                             >  1:1253209/1‑100 (MQ=255)
                       ccAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGc                                                                            >  1:3497016/1‑100 (MQ=255)
                        cAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGcc                                                                           <  1:788857/100‑1 (MQ=255)
                        cAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGcc                                                                           >  2:1430333/1‑100 (MQ=255)
                           cGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGa                                                                        <  2:1840021/100‑1 (MQ=255)
                            gCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAg                                                                       <  2:1951607/100‑1 (MQ=255)
                             ccGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGc                                                                      <  1:2915311/100‑1 (MQ=255)
                             ccGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGc                                                                      <  1:84975/100‑1 (MQ=255)
                             ccGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGc                                                                      <  2:3572477/100‑1 (MQ=255)
                              ngTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCa                                                                     >  2:229609/2‑100 (MQ=255)
                               gTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCAt                                                                    <  2:1271925/100‑1 (MQ=255)
                                   cTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATaaa                                                                 <  2:2462395/99‑1 (MQ=255)
                                   cTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAAt                                                                >  2:3131566/1‑100 (MQ=255)
                                   cTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCCCCTCCCCTGACCTTTTCCCGGAATCAACTCTGCCGGAGCCTaa                                                                  >  2:22414/1‑98 (MQ=37)
                                    tGCAATGCCACCCGCGAAATCGACCTACTCCTTGCCTTGCCGACCCCCGATTCGGCTGCCGTGACCTTTTACCGGAATAAACTCTGACGGCGCATAAATc                                                               >  2:1313817/1‑100 (MQ=255)
                                    tGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATc                                                               <  1:3871391/100‑1 (MQ=255)
                                    tGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATc                                                               <  1:904891/100‑1 (MQ=255)
                                     gCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCg                                                              >  2:371727/1‑100 (MQ=255)
                                     gCAAGGCCACCCGCGAAATCGACATACCCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTCTACCGGAATAAACTCTGCCGGAGCATaaagc                                                               >  2:2645500/1‑97 (MQ=255)
                                        aTGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCa                                                           >  2:3189835/1‑100 (MQ=255)
                                          gCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTGACCGGAATAAACTCTGCCGGAGCATAAATCGGCAgg                                                         <  1:2943196/100‑1 (MQ=255)
                                            cACCCGCGCATTCGACATACTCACTGCCTTGCTGATCCCAGACTCGGCTGCCCTGACCTTTTACAGGAATAAACTCTGCCGGGGCATAGATCCGCACGAt                                                       >  2:3790243/1‑100 (MQ=255)
                                            cACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGAt                                                       >  2:692447/1‑100 (MQ=255)
                                             acacgcgAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGAACTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATc                                                      <  2:1372682/97‑1 (MQ=255)
                                               ccGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCAc                                                    >  2:4094139/1‑100 (MQ=255)
                                                cgctAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGAt                                                       <  2:22986/96‑1 (MQ=255)
                                                 gcgAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACtt                                                  <  2:2987236/100‑1 (MQ=255)
                                                 gcgAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACtt                                                  <  2:750303/100‑1 (MQ=255)
                                                 gcgAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACtt                                                  >  1:1984957/1‑100 (MQ=255)
                                                   gcaaccGACATACTCCTTGCCTTTCTGATCCCAGAATCTGCTGCCCTGACCTCTTACCGGCCTAACCTCTGCCGGATCATAAACCCACAGGATCACtcc                                                 >  2:336638/6‑97 (MQ=255)
                                                    aaaTCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCAt                                               <  1:4108560/100‑1 (MQ=255)
                                                    aaaTCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCAt                                               >  2:195667/1‑100 (MQ=255)
                                                     aaTCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATc                                              <  1:1761975/100‑1 (MQ=255)
                                                       tCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGa                                            <  1:2062741/100‑1 (MQ=255)
                                                         gACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGaaa                                          >  1:3406551/1‑100 (MQ=255)
                                                         gACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGaaa                                          <  2:3486069/100‑1 (MQ=255)
                                                         gACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGaaa                                          >  2:2279557/1‑100 (MQ=255)
                                                         gACATACTCCTTGCCTTGCTGACCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGaaa                                          >  2:3799526/1‑100 (MQ=255)
                                                              aCTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTcgc                                     >  1:508101/1‑100 (MQ=255)
                                                               cTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTcgcg                                    >  1:1084940/1‑100 (MQ=255)
                                                                  cttgccttgcTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCgt                                 >  1:3425445/1‑100 (MQ=255)
                                                                    tgccttgcTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCgtgt                               <  1:1397660/100‑1 (MQ=255)
                                                                     nccttgcTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTa                              >  1:1858307/2‑100 (MQ=255)
                                                                      ccttgcTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTaa                             <  2:674070/100‑1 (MQ=255)
                                                                      ccttgcTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTaa                             >  1:2283736/1‑100 (MQ=255)
                                                                         tgcTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTg                          <  1:2122652/100‑1 (MQ=255)
                                                                           cTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCt                        >  2:176513/1‑100 (MQ=255)
                                                                            tGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCtg                       >  2:826102/1‑100 (MQ=255)
                                                                             gATCCCAGATTCGGCTTCCCTGGCCTTTTACCGGATTAAACTCTGCCGGCGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGCGTAAGTGc                         >  2:234516/1‑97 (MQ=255)
                                                                             gATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCtgt                      <  2:1135487/100‑1 (MQ=255)
                                                                             caTCCCAGATTCGGTTGCCCGGACCTTTTACCGGAATAAAGTCTGCCGGGGCATAAATCGGCAGGATCTCTTCAGTGAAAGTCGCGCGTGTAATTgtgg                       <  2:2521132/98‑4 (MQ=255)
                                                                              aTCCCAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCtgtt                     >  1:1516211/1‑100 (MQ=255)
                                                                              aTCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCtgtt                     <  1:2367187/100‑1 (MQ=255)
                                                                                cccAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCtgtttg                   <  2:3086358/100‑1 (MQ=255)
                                                                                 ccAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCtgtttgt                  >  1:84161/1‑100 (MQ=255)
                                                                                 ccAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAACTCGGCAGGATCACTTCATCGACAGTCGCGCGCGTAATCGCt                        >  2:2442576/1‑94 (MQ=255)
                                                                                 ccAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCtgtttgt                  >  1:594571/1‑100 (MQ=255)
                                                                                  cAGATTCGGCTGCCCTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCtgtttgtt                 <  2:1521642/100‑1 (MQ=255)
                                                                                   aGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTc                <  2:2828554/100‑1 (MQ=255)
                                                                                    gATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCa               <  1:2016662/100‑1 (MQ=255)
                                                                                      ttCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAAt             >  1:1878634/1‑100 (MQ=255)
                                                                                       tCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAAtt            <  2:1742850/100‑1 (MQ=255)
                                                                                        cGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTg           <  2:877231/100‑1 (MQ=255)
                                                                                        cGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTg           <  2:3673518/100‑1 (MQ=255)
                                                                                        cGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTg           <  1:3572355/100‑1 (MQ=255)
                                                                                          gCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGcc         <  1:1840817/100‑1 (MQ=255)
                                                                                             gCCCTGACCTTTTCCCGGAACAAACTCGGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCAtc      <  2:3925183/100‑1 (MQ=255)
                                                                                             gCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCAtc      <  1:112984/100‑1 (MQ=255)
                                                                                               ccTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCAtct     >  2:2488578/1‑99 (MQ=255)
                                                                                                cTGGCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCa   >  1:641118/1‑100 (MQ=255)
                                                                                                cTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCa   >  1:986905/1‑100 (MQ=255)
                                                                                                   aCCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCAt  <  2:50465/98‑1 (MQ=255)
                                                                                                   |                                                                                                 
CACCAGCGCAGGATGGCAATGGCCCAACGCCGTAACTGCAATGCCACCCGCGAAATCGACATACTCCTTGCCTTGCTGATCCCAGATTCGGCTGCCCTGACCTTTTACCGGAATAAACTCTGCCGGAGCATAAATCGGCAGGATCACTTCATCGAAAGTCGCGCGTGTAATTGCTGTTTGTTCAATTGCCATCTCAT  >  NC_000913/3489989‑3490185

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