Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,295,744 G→A 69.3% A422T (GCT→ACT)  gltP → glutamate/aspartate:proton symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,295,7440GA69.3% 107.6 / 80.9 98A422T (GCT→ACT) gltPglutamate/aspartate:proton symporter
Reads supporting (aligned to +/- strand):  ref base G (21/9);  new base A (40/28);  total (61/37)
Fisher's exact test for biased strand distribution p-value = 3.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.66e-01

ACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGA  >  NC_000913/4295649‑4295840
                                                                                               |                                                                                                
aCATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTa                                                                                              <  1:1594671/100‑1 (MQ=255)
   tGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGa                                                                                           >  2:1717607/1‑100 (MQ=255)
    ggCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAg                                                                                          >  2:1481704/1‑100 (MQ=255)
      cgcgTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCg                                                                                        <  1:1253082/100‑1 (MQ=255)
      cgcgTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCg                                                                                        >  2:1593917/1‑100 (MQ=255)
      cgcgTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCg                                                                                        >  1:468623/1‑100 (MQ=255)
      cgcgTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCg                                                                                        >  1:466815/1‑100 (MQ=255)
         gTACTGCGCTGAACGCGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGa                                                                                     <  1:2214385/100‑1 (MQ=255)
           aCTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAg                                                                                   <  2:407124/100‑1 (MQ=255)
           aCTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAg                                                                                   <  2:1213170/100‑1 (MQ=255)
            cTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGt                                                                                  >  2:853795/1‑100 (MQ=255)
            cTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGt                                                                                  <  2:666372/100‑1 (MQ=255)
            cTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGt                                                                                  <  1:697120/100‑1 (MQ=255)
            cTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGt                                                                                  >  1:2157613/1‑100 (MQ=255)
              gcgcTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGc                                                                                >  2:976883/1‑100 (MQ=255)
                 cTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTgg                                                                             >  2:1543899/1‑100 (MQ=255)
                 cTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGCTgg                                                                             >  2:553330/1‑100 (MQ=255)
                  tGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGCTggg                                                                            <  2:686172/100‑1 (MQ=255)
                   gAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGc                                                                           <  1:228284/100‑1 (MQ=255)
                     aCGTGGTGGGTAATGCGCTTGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGCTgg                                                                             <  2:504512/96‑1 (MQ=255)
                        tggtggGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGCTGGGCAAAtt                                                                      >  1:302821/1‑100 (MQ=255)
                         ggtggGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGCTGGGCAAAttt                                                                     >  1:1739733/1‑100 (MQ=255)
                         ggtggGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGCTGGGCAAAttt                                                                     >  2:1956028/1‑100 (MQ=255)
                          gtggGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGCTGGGCAAATTTg                                                                    >  1:1343783/1‑100 (MQ=255)
                            gggTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGAt                                                                  >  1:1688557/1‑100 (MQ=255)
                              gTAAGGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGATaa                                                                <  2:141750/100‑1 (MQ=255)
                                aaTGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGATaaaa                                                              >  1:840129/1‑100 (MQ=255)
                                 aTGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAAc                                                             <  2:449713/100‑1 (MQ=255)
                                 aTGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAAc                                                             <  2:2187750/100‑1 (MQ=255)
                                  tGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACt                                                            >  2:1301346/1‑100 (MQ=255)
                                  tGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACt                                                            >  1:479717/1‑100 (MQ=255)
                                   gcgcTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGGAAGAAAGAGATGACTACTGAGCTTGAAATGATGGGCAAAGTTGATacaa                                                              >  2:2163652/1‑97 (MQ=255)
                                    cgcTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGc                                                          >  1:1168598/1‑100 (MQ=255)
                                    cgcTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGc                                                          >  2:11422/1‑100 (MQ=255)
                                      cTGGCGGTGATGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGTTGACTTATGAGCGTGAAGTGTTGTGCAAATTTGATAAAACTGCg                                                         <  2:1945082/99‑1 (MQ=255)
                                        ggCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGAt                                                      <  2:590219/100‑1 (MQ=255)
                                         gCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATc                                                     <  2:1248207/100‑1 (MQ=255)
                                          cGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCa                                                    >  1:2149301/1‑100 (MQ=255)
                                          cGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCa                                                    >  2:264365/1‑100 (MQ=255)
                                          cGGTGCTGGTGATTACCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCa                                                    >  2:1224075/1‑100 (MQ=255)
                                           ggTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCaa                                                   <  2:1188409/100‑1 (MQ=255)
                                              gCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCaataa                                                >  2:1647443/1‑100 (MQ=255)
                                              gCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCaataa                                                >  1:920453/1‑100 (MQ=255)
                                              gCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCaataa                                                >  1:1234592/1‑100 (MQ=255)
                                                tGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAAtt                                              >  1:31542/1‑100 (MQ=255)
                                                tGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAACCGCGGAATTATGAGCGTGAAGCACTATGCAAAATTGATAAAAATGCGGATCAATAAtt                                              >  2:126649/1‑100 (MQ=255)
                                                 ggTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTg                                             <  1:2184398/100‑1 (MQ=255)
                                                 ggTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTg                                             <  1:728343/100‑1 (MQ=255)
                                                 ggTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAAATGCGGATCAATAATTg                                             >  2:191067/1‑100 (MQ=255)
                                                 ggTGATTGCCAAGAGGGAACACAAATTTGACCGTAAGAGAGCGCTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTg                                             >  2:1383925/1‑100 (MQ=255)
                                                   tGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGaa                                           <  2:209705/100‑1 (MQ=255)
                                                   tGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGaa                                           >  2:1657533/1‑100 (MQ=255)
                                                    gATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAg                                          >  1:326672/1‑100 (MQ=255)
                                                     aTTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGAATTATGAGCGTGAAGTGCTAGGCAAATCTGATAAAACTGCGGATCAATAAttta                                            >  2:1914397/1‑95 (MQ=255)
                                                      ttGCCAAGTGGGTACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGAt                                        >  2:1005141/1‑100 (MQ=255)
                                                      ttGCCAAGTGGGTACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGAt                                        >  2:242049/1‑100 (MQ=255)
                                                      ttGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGGGCGTGAAGGGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGAt                                        <  1:2124494/100‑1 (MQ=255)
                                                       tGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGATGAATTATAAGCGTGAAGTGCTGGGCAAATTTGATAACACTGCGGATCAATAATTGAAGAtt                                       <  1:660751/100‑1 (MQ=255)
                                                        gCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGaat                                       >  2:72631/1‑97 (MQ=255)
                                                        gCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCAGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTg                                      >  2:1194977/1‑100 (MQ=255)
                                                         ccAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGc                                     >  2:1974422/1‑100 (MQ=255)
                                                         ccAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGc                                     >  1:1058451/1‑100 (MQ=255)
                                                          cAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGAGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGcc                                    >  2:931553/1‑100 (MQ=255)
                                                             gTGGGAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCgg                                  <  1:170377/99‑1 (MQ=255)
                                                               gggAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGa                               <  1:996472/100‑1 (MQ=255)
                                                               gggAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGa                               <  2:1800824/100‑1 (MQ=255)
                                                                ggAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGat                              <  1:1490268/100‑1 (MQ=255)
                                                                ggAACACAAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGat                              <  1:1723182/100‑1 (MQ=255)
                                                                    cacaAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATAAATAATTGAAGATTGCCGGGGATATcc                          >  1:2030258/1‑100 (MQ=255)
                                                                    cacaAATTTGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATcc                          >  1:227655/1‑100 (MQ=255)
                                                                     acaAATTTGACCGTAAGAAAGAGCTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAAATGAAGATTGCCGGGGATATCCa                         >  2:314142/1‑100 (MQ=255)
                                                                      caAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCAc                        >  2:822676/1‑100 (MQ=255)
                                                                      caAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCAc                        >  2:2003239/1‑100 (MQ=255)
                                                                          tttGACCGTAAGAGAGCGCTGACTGATGAGCGTGAAGGGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCgg                    <  1:1278208/100‑1 (MQ=255)
                                                                          tttGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCg                     >  1:300953/1‑99 (MQ=255)
                                                                           ttGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGc                   <  2:2157598/100‑1 (MQ=255)
                                                                           ttGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGc                   <  2:841384/100‑1 (MQ=255)
                                                                            tGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCa                  >  1:1687847/1‑100 (MQ=255)
                                                                            tGACCGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCa                  >  1:1309631/1‑100 (MQ=255)
                                                                              aCCGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAAt                <  1:1889912/100‑1 (MQ=255)
                                                                               ccGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGCTATCCACCCGGCAAtg               >  1:1023839/1‑100 (MQ=255)
                                                                               ccGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAAtg               >  2:274058/1‑100 (MQ=255)
                                                                               ccGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAAtg               <  2:920452/100‑1 (MQ=255)
                                                                               ccGTAAGAAAGCGCTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAAt                >  2:1868002/1‑99 (MQ=255)
                                                                                cGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAAtgt              >  2:258024/1‑100 (MQ=255)
                                                                                   aaGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAAtgtgtg           <  1:568559/100‑1 (MQ=255)
                                                                                     gAAAGCGCTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGaa         >  1:758209/1‑100 (MQ=255)
                                                                                       aaGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATg       >  2:706807/1‑100 (MQ=255)
                                                                                         gcgcTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGcc     >  2:1291520/1‑100 (MQ=255)
                                                                                         gcgcTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGcc     <  1:1755788/100‑1 (MQ=255)
                                                                                         gcgcTGACTTATGAGCGTGAAGTGCTGGGCAAATGTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGcc     <  1:1641437/100‑1 (MQ=255)
                                                                                          cgcTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGAAAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCt    >  2:1199674/1‑100 (MQ=255)
                                                                                           gcTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTg   <  1:1832504/100‑1 (MQ=255)
                                                                                           gcTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTg   >  2:536075/1‑100 (MQ=255)
                                                                                            cTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTg   <  1:1830789/99‑1 (MQ=255)
                                                                                            cTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGa  >  2:1667799/1‑100 (MQ=255)
                                                                                            cTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGa  >  2:1393336/1‑100 (MQ=255)
                                                                                            cTGACTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGa  >  1:1847851/1‑100 (MQ=255)
                                                                                               |                                                                                                
ACATGGCGCGTACTGCGCTGAACGTGGTGGGTAATGCGCTGGCGGTGCTGGTGATTGCCAAGTGGGAACACAAATTTGACCGTAAGAAAGCGCTGGCTTATGAGCGTGAAGTGCTGGGCAAATTTGATAAAACTGCGGATCAATAATTGAAGATTGCCGGGGATATCCACCCGGCAATGTGTGAATGCCTGA  >  NC_000913/4295649‑4295840

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