Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,132,227 C→A A798E (GCG→GAG)  metL → Bifunctional aspartokinase/homoserine dehydrogenase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,132,2270CA97.2% 712.4 / 3.7 214A798E (GCG→GAG) metLBifunctional aspartokinase/homoserine dehydrogenase 2
Reads supporting (aligned to +/- strand):  ref base C (4/2);  new base A (105/103);  total (109/105)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 216 total aligned reads displayed.

CGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAAT  >  NC_000913/4132131‑4132322
                                                                                                |                                                                                               
  ttcGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACATCACCGCCGGGGAGAtt                                                                                             <  1:479744/97‑1 (MQ=255)
  tGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGAtt                                                                                             <  2:1743068/99‑1 (MQ=255)
   gCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTc                                                                                            >  1:2614326/1‑99 (MQ=255)
   gCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTc                                                                                            >  1:2328616/1‑99 (MQ=255)
    cGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCa                                                                                           <  1:1435931/99‑1 (MQ=255)
     gATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAg                                                                                          >  2:1589009/1‑99 (MQ=255)
     gATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAg                                                                                          <  1:2692275/99‑1 (MQ=255)
      aTAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGt                                                                                         <  2:2679099/99‑1 (MQ=255)
       tAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTc                                                                                        <  1:928767/99‑1 (MQ=255)
        aaCGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCg                                                                                       <  2:1344738/99‑1 (MQ=255)
         aCGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCgg                                                                                      >  1:4369985/1‑99 (MQ=255)
         aCGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCgg                                                                                      >  1:4634175/1‑99 (MQ=255)
             cTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGatat                                                                                  <  1:4101254/99‑1 (MQ=255)
              tttGCCATCGAAAGCCACTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATc                                                                                 >  2:4873750/1‑99 (MQ=255)
                tGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCaa                                                                               <  2:1388101/99‑1 (MQ=255)
                tGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCaa                                                                               <  1:4582554/99‑1 (MQ=255)
                 gCCATCGAAAGCTGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAAc                                                                              >  1:4964045/1‑99 (MQ=255)
                 gCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAAc                                                                              <  1:2561075/99‑1 (MQ=255)
                 gCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAAc                                                                              <  2:4515049/99‑1 (MQ=255)
                 gCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAAc                                                                              <  1:5606573/99‑1 (MQ=255)
                  ccATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAAcc                                                                             >  1:3421852/1‑99 (MQ=255)
                   cATCGAAAGCCGCTGGTATCGCGATAACCCACTGGTGATCCGCGGCCCTGGCGCTGGGCGCGACGCGACCGCCGGGGAGATTCAGTCGGATATCAACCg                                                                            >  2:3920046/1‑99 (MQ=255)
                    aTCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCCGTCGGATATCAACCgg                                                                           >  2:4144088/1‑99 (MQ=255)
                     tCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGc                                                                          >  1:107480/1‑99 (MQ=255)
                     tCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGc                                                                          >  1:641715/1‑99 (MQ=255)
                     tCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGAGCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGc                                                                          <  2:2870126/99‑1 (MQ=255)
                      cGAAGGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCt                                                                         <  1:659926/99‑1 (MQ=255)
                      cGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCt                                                                         <  1:2574860/99‑1 (MQ=255)
                         aaGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGc                                                                      <  2:3701241/99‑1 (MQ=255)
                           gCCGCTGGTATCGCGATAACCCTGTGGTTATCCGCTGACCTGGCGCGGTGCGCGACGTCACCGCCGTGTAGGTTCTGTCGTATGTCAGCCTGCTGGcac                                                                    >  2:3033083/1‑99 (MQ=255)
                             cGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAg                                                                  >  2:442138/1‑99 (MQ=255)
                             cGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAg                                                                  >  1:517985/1‑99 (MQ=255)
                             cGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAg                                                                  <  1:3490707/99‑1 (MQ=255)
                               cTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAgtt                                                                >  2:2258586/1‑99 (MQ=255)
                               cTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAgtt                                                                >  2:3665737/1‑99 (MQ=255)
                               cTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAgtt                                                                <  2:1900041/99‑1 (MQ=255)
                                tGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAgttg                                                               <  2:4517864/99‑1 (MQ=255)
                                 ggcaTCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAgttgt                                                              <  1:2267813/96‑1 (MQ=255)
                                  gTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAgttgtt                                                             >  1:2168876/1‑99 (MQ=255)
                                  gTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAgttgtt                                                             >  2:3419317/1‑99 (MQ=255)
                                  gTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAgttgtt                                                             >  2:3275934/1‑99 (MQ=255)
                                     tCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTa                                                          >  1:1871174/1‑99 (MQ=255)
                                     tCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGACGGATATCAACCGGCTGGCACAGTTGTTGTa                                                          >  1:1490842/1‑99 (MQ=255)
                                      cgcgATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTaa                                                         <  1:3873179/99‑1 (MQ=255)
                                       gcgATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCCACCGGCTGGCACAGGTGTTGTAAt                                                        >  2:1830999/1‑99 (MQ=255)
                                           tagcccTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCa                                                    <  2:4722661/96‑1 (MQ=255)
                                           gaaCCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCa                                                    <  2:1124543/98‑1 (MQ=255)
                                            aaCCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAg                                                   <  1:449584/99‑1 (MQ=255)
                                             aCCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGa                                                  <  2:2143338/99‑1 (MQ=255)
                                              cccTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGaa                                                 <  2:649595/99‑1 (MQ=255)
                                               ccTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGaaa                                                >  1:2950892/1‑99 (MQ=255)
                                                  ctGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAAttt                                             >  1:3825558/1‑99 (MQ=255)
                                                  ctGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAAttt                                             >  1:3602300/1‑99 (MQ=255)
                                                   tGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTta                                            <  2:3087973/99‑1 (MQ=255)
                                                    ggTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTtaa                                           >  1:3323465/1‑99 (MQ=255)
                                                     gTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTtaat                                          <  2:4370427/99‑1 (MQ=255)
                                                       gATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTtaataa                                        >  1:2479562/1‑99 (MQ=255)
                                                        aTCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTaataat                                       >  1:1211869/1‑99 (MQ=255)
                                                        aTCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACCGTTGTTGTAATTTCAGAAATTTaataat                                       >  2:4584419/1‑99 (MQ=255)
                                                        aTCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTaataat                                       <  1:2032955/99‑1 (MQ=255)
                                                          ccGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGc                                     <  1:4641891/99‑1 (MQ=255)
                                                          ccGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGc                                     <  1:5377781/99‑1 (MQ=255)
                                                           cgcgGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGcc                                    <  2:2254045/99‑1 (MQ=255)
                                                           cgcgGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGcc                                    <  2:5060458/99‑1 (MQ=255)
                                                           cgcgGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGcc                                    >  1:5360885/1‑99 (MQ=255)
                                                             cgGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCg                                  <  1:3381387/99‑1 (MQ=255)
                                                             cgGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCg                                  <  2:4383792/99‑1 (MQ=255)
                                                                aCCTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTa                               >  1:2625116/1‑99 (MQ=255)
                                                                 ccTGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTAc                              <  2:3794719/99‑1 (MQ=255)
                                                                   tGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTc                            >  1:5066188/1‑99 (MQ=255)
                                                                   tGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTc                            >  1:817269/1‑99 (MQ=255)
                                                                   tGGCGCTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTc                            >  1:5315545/1‑99 (MQ=255)
                                                                    ggCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCa                           >  2:3996341/1‑99 (MQ=255)
                                                                     gcgcTGGGCGCGACGTCACCGCCGGGGAGGTTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCAt                          <  2:767191/99‑1 (MQ=255)
                                                                      cgcTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATg                         >  2:129259/1‑99 (MQ=255)
                                                                       gcTGGGCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGt                        >  1:169194/1‑99 (MQ=255)
                                                                          gggCGCGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGtttt                     >  2:4129478/1‑99 (MQ=255)
                                                                           ggCGCGACGTCACCGCCGGGGAGATTCAGTCGGCTATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTc                    <  1:5579367/99‑1 (MQ=255)
                                                                            gcgcGACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCCTGTTTTCg                   >  1:1359060/1‑99 (MQ=255)
                                                                             cgcgACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCgg                  <  2:23668/99‑1 (MQ=255)
                                                                             cgcgACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCgg                  <  1:4908940/99‑1 (MQ=255)
                                                                                gACGTCCCCGCCGGGGAGATTCAGTCCGATATCAACCGACTGGCACAGTTGCTGTAATTTCAGAAATTTAATAATGCCCGGTACGCACGTTTTCGGGtt               >  2:2793030/1‑99 (MQ=255)
                                                                                gACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGtt               >  1:3702983/1‑99 (MQ=255)
                                                                                gACGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGtt               >  1:9378/1‑99 (MQ=255)
                                                                                 aCGTCGCCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAGTAATGCCCGGTACTCATGTTTTCGGGttt              <  1:991237/99‑1 (MQ=255)
                                                                                 aCGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGttt              >  1:509056/1‑99 (MQ=255)
                                                                                  cGTCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTa             <  1:4992414/99‑1 (MQ=255)
                                                                                    tCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATg           >  1:319499/1‑99 (MQ=255)
                                                                                    tCACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATg           >  1:322217/1‑99 (MQ=255)
                                                                                     cACCGCCGGGGAGATTCAGTTGCATATCAACAGGTTGGCAGAGTTGTTGTAATTTCAGTAATTTAATAATGTCCGGTACTCATGTTTTCGGGTTTATgg          <  2:1439669/99‑1 (MQ=255)
                                                                                     cACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATgg          <  1:5462074/99‑1 (MQ=255)
                                                                                     cACCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATgg          <  2:2169421/99‑1 (MQ=255)
                                                                                      aCCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGt         <  2:1258080/99‑1 (MQ=255)
                                                                                      aCCGCCGGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGATGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGCTTATGGt         >  2:2456569/1‑99 (MQ=255)
                                                                                       ccgccgGGGAGATTCAGTCGGATATCGACCGGCTGGCACAGTTGTTGGAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATggat        >  2:5495798/1‑97 (MQ=255)
                                                                                        cgccgGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGttt       <  1:1180699/99‑1 (MQ=255)
                                                                                         gccgGGGAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTc      <  1:4743200/99‑1 (MQ=255)
                                                                                            ggggAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTaa   <  2:2972154/99‑1 (MQ=255)
                                                                                            ggggAGATTCAGTCGGAGATCAACCGGCTGGCACAGTTGTTGTTATTTCAGAAATTTAATAATGCCTGGTACTCATGTTTTCGCGTTTATGGTTTCTaa   <  2:2010935/99‑1 (MQ=255)
                                                                                             gggAGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAAt  >  1:4186680/1‑99 (MQ=255)
                                                                                                |                                                                                               
CGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAAT  >  NC_000913/4132131‑4132322

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: