Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,737,998 T→C 20.1% V85A (GTA→GCA)  pheA → chorismate mutase and prephenate dehydratase, P‑protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,737,9980TC20.1% 237.7 / 54.7 114V85A (GTA→GCA) pheAchorismate mutase and prephenate dehydratase, P‑protein
Reads supporting (aligned to +/- strand):  ref base T (43/48);  new base C (9/14);  total (52/62)
Fisher's exact test for biased strand distribution p-value = 6.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 116 total aligned reads displayed.

AACGCGATTTGCTGGAAAGATTAATTACGCTCGGTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATCCGCACTCAGCACGCATCGCTTTTCTCGGCCCCAAAGGTTCTTATTCCCATC  >  NC_000913/2737899‑2738096
                                                                                                   |                                                                                                  
aaCGCGATTTGCTGGAAAGATTAATTACGCTCGGTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGt                                                                                                    <  2:72866/100‑1 (MQ=255)
aaCGCGATTTGCTGGAAAGATTAATTACGCTCGGTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGt                                                                                                    <  2:2322503/100‑1 (MQ=255)
 aCGCGATTTGCTGGAAAGATTAATTACGCTCGGTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTa                                                                                                   <  2:1414946/100‑1 (MQ=255)
 aCGCGATTTGCTGGAAAGATTAATTACGCTCGGTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTa                                                                                                   <  2:438900/100‑1 (MQ=255)
  cgcgATTTGCTGGAAAGATTAATTACGCTCGGTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTAt                                                                                                  >  2:3145083/1‑100 (MQ=255)
  cgcgATTTGCTGGAAAGATTAATTACGCTCGGTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTAt                                                                                                  >  1:1544552/1‑100 (MQ=255)
   gcgATTTGCTGGAAAGATTAATTACGCTCGGTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTAtt                                                                                                 <  2:1639076/100‑1 (MQ=255)
    cgATTTGCTGGAAAGATTAATTACGCTCGGTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTa                                                                                                >  1:1301423/1‑100 (MQ=255)
      aTTTGCTGGAAAGATTAATTACGCTCGGTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAAc                                                                                              <  1:1107862/100‑1 (MQ=255)
      aTTTGCTGGAAAGATTAATTACGCTCGGTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAAc                                                                                              >  2:861164/1‑100 (MQ=255)
           cTGGAAAGATTAATTACGCTCGGTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTcagc                                                                                         >  2:460009/1‑100 (MQ=255)
            tGGAAAGATTAATTACGCTCGGTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTcagca                                                                                        >  2:885578/1‑100 (MQ=255)
             ggAAAGATTAATTACGCTCGGTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTcagcag                                                                                       <  2:1647000/100‑1 (MQ=255)
             ggAAAGATTAATTACGCTCGGTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTcagcag                                                                                       <  1:71812/100‑1 (MQ=255)
              gAAAGATTAATTACGCTCGGTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAgg                                                                                      >  2:1894141/1‑100 (MQ=255)
                 aGATTAATTACGCTCGGTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGCATTAACTCAGCAGGCtt                                                                                   >  2:281495/1‑100 (MQ=255)
                 aGATTAATTACGCTCGGTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGCATTAACTCAGCAGGCtt                                                                                   <  1:2913597/100‑1 (MQ=255)
                  gATTAATTACGCTCGGTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCttt                                                                                  <  1:2433190/100‑1 (MQ=255)
                   attaattaCGCTCGGTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTg                                                                                 >  1:2747703/1‑100 (MQ=255)
                    ttaattaCGCTCGGTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGc                                                                                <  2:589510/100‑1 (MQ=255)
                    ttaattaCGCTCGGTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGCATTAACTCAGCAGGCTTTGc                                                                                <  2:3267150/100‑1 (MQ=255)
                    ttaattaCGCTCGGTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGCATTAACTCAGCAGGCTTTGc                                                                                <  1:2480537/100‑1 (MQ=255)
                      aattaCGCTCGGTAAAGCGCACCATCTGGACGCCCATTGCATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTc                                                                              <  2:1788194/100‑1 (MQ=255)
                      aattaCGCTCGGTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTc                                                                              >  1:1385976/1‑100 (MQ=255)
                       attaCGCTCGGTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTcc                                                                             >  1:2665301/1‑100 (MQ=255)
                          aCGCTCGGTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGCATTAACTCAGCAGGCTTTGCTCcaac                                                                          <  1:281497/100‑1 (MQ=255)
                           cGCTCGGTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCcaaca                                                                         >  2:3043769/1‑100 (MQ=255)
                           cGCTCGGTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGCATTAACTCAGCAGGCTTTGCTCcaaca                                                                         >  1:2089608/1‑100 (MQ=255)
                            gCTCGGTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCcaacaa                                                                        <  2:1557906/100‑1 (MQ=255)
                            gCTCGGTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCcaacaa                                                                        >  2:3020362/1‑100 (MQ=255)
                            gCTCGGTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCcaacaa                                                                        >  2:2189992/1‑100 (MQ=255)
                             cTCGGTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGCATTAACTCAGCAGGCTTTGCTCCaacaac                                                                       >  2:464094/1‑100 (MQ=255)
                              tCGGTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCGGCAGGCTTTGCTCCAacaaca                                                                      <  1:2348404/100‑1 (MQ=255)
                              tCGGTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAacaaca                                                                      <  2:344274/100‑1 (MQ=255)
                               cGGTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACAt                                                                     <  2:3048646/100‑1 (MQ=255)
                                ggTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACAtc                                                                    >  1:1836098/1‑100 (MQ=255)
                                ggTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACAtc                                                                    <  1:2910691/100‑1 (MQ=255)
                                ggTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACAt                                                                     >  2:2943014/1‑99 (MQ=255)
                                 gTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACAtct                                                                   >  2:88539/1‑100 (MQ=255)
                                 gTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGCATTAACTCAGCAGGCTTTGCTCCAACAACAtct                                                                   >  1:83869/1‑100 (MQ=255)
                                  tAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACAtctc                                                                  <  1:1434281/100‑1 (MQ=255)
                                  tAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACAtctc                                                                  >  1:792614/1‑100 (MQ=255)
                                   aaaGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCa                                                                 >  1:1308065/1‑100 (MQ=255)
                                   aaaGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGCATTAACTCAGCAGGCTTTGCTCCAACAACATCTCa                                                                 >  2:2649515/1‑100 (MQ=255)
                                    aaGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCaa                                                                <  2:2497586/100‑1 (MQ=255)
                                    aaGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGCATTAACTCAGCAGGCTTTGCTCCAACAACATCTCaa                                                                <  1:1038684/100‑1 (MQ=255)
                                     aGCGAACCATCTGGACGCCCATTACATTAATCGCCTGTTCAAGCTCATCCTTGAAGATTCCGTATTAACTCAGAAGGCTTTCCTCCACCAACATCTAAAt                                                               >  2:634807/1‑100 (MQ=255)
                                       cgcACCATCTGGACGACCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATaa                                                             <  2:1213824/100‑1 (MQ=255)
                                          aCCATCTGGACGCCCATTACATTACTCTCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAAt                                                          <  2:1950470/100‑1 (MQ=255)
                                          aCCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAAt                                                          <  1:903479/100‑1 (MQ=255)
                                          aCCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAAt                                                          <  1:1898440/100‑1 (MQ=255)
                                          aCCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGCATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAAt                                                          <  2:1374453/100‑1 (MQ=255)
                                             aTCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTaa                                                       <  1:2336396/100‑1 (MQ=255)
                                               ctttaCGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGCATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATc                                                     <  2:2877267/96‑1 (MQ=255)
                                               cTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATc                                                     <  2:802129/100‑1 (MQ=255)
                                               cTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATc                                                     <  2:1062043/100‑1 (MQ=255)
                                               cTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATc                                                     <  1:1749432/100‑1 (MQ=255)
                                               cTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGCATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATc                                                     <  1:3080149/100‑1 (MQ=255)
                                                tGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATcc                                                    <  2:2665288/100‑1 (MQ=255)
                                                tGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATcc                                                    >  1:218126/1‑100 (MQ=255)
                                                tGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATcc                                                    <  2:2957294/100‑1 (MQ=255)
                                                 ggACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATCCg                                                   <  2:1446208/100‑1 (MQ=255)
                                                 ggACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATCCg                                                   >  1:3086753/1‑100 (MQ=255)
                                                  gACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATCCGc                                                  <  2:792612/100‑1 (MQ=255)
                                                  gACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGCATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATCCGc                                                  <  1:2649527/100‑1 (MQ=255)
                                                   aCGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATCCGc                                                  >  2:2877939/1‑99 (MQ=255)
                                                   aCGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATCCGCa                                                 <  2:1544548/100‑1 (MQ=255)
                                                     gccggttacattacTCGCCTGTACCCGCTCAGCACTGAAGATTCCGCATTAAGTCAGCAGGCTTTGCTTCAACGACATCTCAAGAAAAt                                                          <  2:3101048/84‑1 (MQ=255)
                                                       ccATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATCCGCACTCa                                             >  1:2950365/1‑100 (MQ=255)
                                                        cattacattaCTCGCCTGTTCCAGCTCATCATTGAAGATTCCGCATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATCCGCACTCAg                                            >  2:868594/1‑100 (MQ=255)
                                                           tacattacTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATCCGCACTCAGCAc                                         >  1:403075/1‑100 (MQ=255)
                                                           tacattacTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATCCGCACTCAGCAc                                         <  1:2092402/100‑1 (MQ=255)
                                                                  cTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATCCGCACTCAGCACGCATCGc                                  <  1:1632095/100‑1 (MQ=255)
                                                                     gCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATCCGCACTCAGCACGCATCGCttt                               <  2:207591/100‑1 (MQ=255)
                                                                       cTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATCCGCACTCAGCACGCATCGCTTTtc                             >  2:1145068/1‑100 (MQ=255)
                                                                        tGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATCCGCACTCAGCACGCATCGCTTTtct                            <  1:790999/100‑1 (MQ=255)
                                                                        tGTTCCAGCTCATCATTGAAGATTCCGCATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATCCGCACTCAGCACGCATCGCTTTtct                            >  1:474765/1‑100 (MQ=255)
                                                                          ttCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACCTCTCAATAAAATTAATCCGCACTCAGCACGCATCGCTTTTCTCg                          <  1:207360/100‑1 (MQ=255)
                                                                          ttCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATCCGCACTCAGCACGCATCGCTTTTCTCg                          >  2:771271/1‑100 (MQ=255)
                                                                          ttCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATCCGCACTCAGCACGCATCGCTTTTCTCg                          <  1:1774030/100‑1 (MQ=255)
                                                                          ttCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATCCGCACTCAGCACGCATCGCTTTTCTCg                          >  2:2139867/1‑100 (MQ=255)
                                                                          ttCCAGCTCATCATTGAAGATTCCGCATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATCCGCACTCAGCACGCATCGCTTTTCTCg                          <  1:2767456/100‑1 (MQ=255)
                                                                            ccAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATCCGCACTCAGCACGCATCGCTTTTCTCGGc                        >  2:1100425/1‑100 (MQ=255)
                                                                               gCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATCCGCACTCAGCACGCATCGCTTTTCTCGGcccc                     >  1:2243612/1‑100 (MQ=255)
                                                                               gCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATCCGCACTCAGCACGCATCGCTTTTCTCGGcccc                     <  1:2538787/100‑1 (MQ=255)
                                                                               gCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATCCGCACTCAGCACGCATCGCTTTTCTCGGcccc                     >  2:865650/1‑100 (MQ=255)
                                                                                  catcatTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATCCGCACTCAGCACGCATCGCTTTTCTCGGCCCCaaa                  >  2:1183761/1‑100 (MQ=255)
                                                                                   atcatTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATCCGCACTCAGTACGCATCGCTTTTCTCGGCCCCAAAg                 <  1:1989921/100‑1 (MQ=255)
                                                                                     catTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATCCGCACTCAGCACGCATCGCTTTTCTCGGCCCCAAAGGt               <  1:992583/100‑1 (MQ=255)
                                                                                        tGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATCCGCACTCAGCACGCATCGCTTTTCTCGGCCCCAAAGGTTCt            <  1:531272/100‑1 (MQ=255)
                                                                                        tGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATCCGCACTCAGCACGCATCGCTTTTCTCGGCCCCAAAGGTTCt            <  1:2115192/100‑1 (MQ=255)
                                                                                         gAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATCCGCACTCAGCACGCATCGCTTTTCTCGGCCCCAAAGGTTCtt           <  2:418029/100‑1 (MQ=255)
                                                                                          aaGATTCCGCATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATCCGCACTCAGCACGCATCGCTTTTCTCGGCCCCAAAGGTTCTTa          <  2:2908375/100‑1 (MQ=255)
                                                                                            gATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATCCGCACTCAGCACGCATCGCTTTTCTCGGCCCCAAAGGTTCTTAtt        >  1:2131363/1‑100 (MQ=255)
                                                                                             aTTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATCCGCACTCAGCACGCATCGCTTTTCTCGGCCCCAAAGGTTCTTATTc       <  1:2190001/100‑1 (MQ=255)
                                                                                             aTTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATCCGCACTCAGCACGCATCGCTTTTCTCGGCCCCAAAGGTTCTTATTc       <  1:1895253/100‑1 (MQ=255)
                                                                                              ttCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATCCGCACTCAGCACGCATCGCTTTTCTCGGCCCCAAAGGTTCTTATTcc      <  2:661859/100‑1 (MQ=255)
                                                                                              ttCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATCCGCACTCAGCACGCATCGCTTTTCTCGGCCCCAAAGGTTCTTATTcc      >  2:8865/1‑100 (MQ=255)
                                                                                                ccGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATCCGCACTCAGCACGCATCGCTTTTCTCGGCCCCAAAGGTTCTTATTCCCa    >  1:2254006/1‑100 (MQ=255)
                                                                                                  gTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATCCGCACTCAGCACGCATCGCTTTTCTCGGCCCCAAAGGTTCTTATTCCCATc  >  2:526624/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
AACGCGATTTGCTGGAAAGATTAATTACGCTCGGTAAAGCGCACCATCTGGACGCCCATTACATTACTCGCCTGTTCCAGCTCATCATTGAAGATTCCGTATTAACTCAGCAGGCTTTGCTCCAACAACATCTCAATAAAATTAATCCGCACTCAGCACGCATCGCTTTTCTCGGCCCCAAAGGTTCTTATTCCCATC  >  NC_000913/2737899‑2738096

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