Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,632,584 A→T 16.7% intergenic (‑49/+20) guaA ← / ← guaB GMP synthetase (glutamine aminotransferase)/IMP dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,632,5840AT16.7% 459.7 / 81.6 192intergenic (‑49/+20)guaA/guaBGMP synthetase (glutamine aminotransferase)/IMP dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (77/81);  new base T (13/19);  total (92/100)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.23e-01

Only 100 of 194 total aligned reads displayed.

GAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGC  >  NC_000913/2632486‑2632683
                                                                                                  |                                                                                                   
gaACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACAt                                                                                                    >  2:4040334/1‑100 (MQ=255)
 aaCCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATg                                                                                                   <  2:664610/100‑1 (MQ=255)
   ccGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCACGAGAAACAGATAATATAAATCGCCCGACCTGaa                                                                                                 >  2:2568462/1‑100 (MQ=255)
   ccGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATAGCCCGACATGa                                                                                                  >  1:4103736/1‑99 (MQ=255)
      aaGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTc                                                                                              <  1:1219487/100‑1 (MQ=255)
        gTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCgg                                                                                            >  2:976921/1‑100 (MQ=255)
        gTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATGCCGCGGCAAGTGAAACAGATAATATACATCGCCAGACATGAAGTggg                                                                                            >  1:4180662/1‑100 (MQ=255)
          ccAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACTTCATGTCGGGc                                                                                          <  1:1075822/100‑1 (MQ=255)
          ccAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACTTCATGTCGGGc                                                                                          <  2:3870962/100‑1 (MQ=255)
          ccAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACTTCATGTCGGGc                                                                                          <  2:3438389/100‑1 (MQ=255)
          ccAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACTTCATGTCGGGc                                                                                          <  2:2681060/100‑1 (MQ=255)
          ccAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACTTCATGTCGGGc                                                                                          <  2:1719161/100‑1 (MQ=255)
          ccAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACTTCATGTCGGGc                                                                                          <  1:708342/100‑1 (MQ=255)
          ccAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACTTCATGTCGGGc                                                                                          <  1:1978203/100‑1 (MQ=255)
          ccAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACTTCATGTCGGGc                                                                                          <  1:1652485/100‑1 (MQ=255)
          ccAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACTTCACGTCGGGc                                                                                          <  2:2324490/100‑1 (MQ=255)
          ccAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGAAGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACTTCATGTCGGGc                                                                                          <  2:186277/100‑1 (MQ=255)
           cAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCg                                                                                         <  1:731658/100‑1 (MQ=255)
           cAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCg                                                                                         >  1:4292065/1‑100 (MQ=255)
            aGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACTTCATGTCGGGCGa                                                                                        >  2:2854404/1‑100 (MQ=255)
            aGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGa                                                                                        >  2:1932907/1‑100 (MQ=255)
             gcatcaGGATGTGATGCTTATGTATGTTTTCTGTCATCGACGCTTATACCGTGGCAAGTGAAACAGATAATATAAATCGCCCGACTTCATGTCGGGCGat                                                                                       <  1:2977098/95‑2 (MQ=255)
             gAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGaa                                                                                       >  2:3185823/1‑100 (MQ=255)
              aaTGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATCAATCGCCCGACATGAAGTCGGGCGaa                                                                                       >  1:280660/1‑99 (MQ=255)
              aaTGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAag                                                                                      <  1:2075319/100‑1 (MQ=255)
              aaTGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAag                                                                                      <  2:3511295/100‑1 (MQ=255)
              aaTGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAag                                                                                      <  2:3253025/100‑1 (MQ=255)
               aTGAGGATGCTATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACTTCATGTCGGGCGat                                                                                       <  1:2520849/98‑2 (MQ=255)
               aTGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAaga                                                                                     >  2:49957/1‑100 (MQ=255)
                tGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATGCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAagag                                                                                    <  1:3205270/100‑1 (MQ=255)
                 gAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAgaga                                                                                   >  2:3653855/1‑100 (MQ=255)
                  aGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGaa                                                                                  >  2:2205339/1‑100 (MQ=255)
                  aGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGaa                                                                                  >  2:419478/1‑100 (MQ=255)
                     atgcgatgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCa                                                                               <  2:2827626/100‑1 (MQ=255)
                      tgcgatgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGAAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAg                                                                              <  1:2910438/100‑1 (MQ=255)
                       gcgatgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAgg                                                                             <  1:478751/100‑1 (MQ=255)
                        cgatgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGa                                                                            >  2:3580348/1‑100 (MQ=255)
                         gatgcTTATGACTGATTGCCGTCAGTGGCGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATcn                                                                               <  2:1656377/96‑2 (MQ=255)
                          atgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGc                                                                          <  2:451823/100‑1 (MQ=255)
                           tgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGcc                                                                         <  2:3000504/100‑1 (MQ=255)
                            gcTTGTGAATGTTTTCCGTCATTGACGCTTAGGCCGAGTCAAGTGAAACAGATAATATAGATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGccc                                                                        <  1:1436904/100‑1 (MQ=255)
                            gcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGccc                                                                        <  1:834304/100‑1 (MQ=255)
                            gcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGccc                                                                        >  1:2120987/1‑100 (MQ=255)
                              ttATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGccc                                                                        <  1:3859969/98‑1 (MQ=255)
                              ttATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAg                                                                      <  2:120564/100‑1 (MQ=255)
                                 tGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACg                                                                   >  1:1277863/1‑100 (MQ=255)
                                 tGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACg                                                                   >  2:3512550/1‑100 (MQ=255)
                                   aaTGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGt                                                                 >  2:596854/1‑100 (MQ=255)
                                   aaTGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGt                                                                 >  2:3141538/1‑100 (MQ=255)
                                      gTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGt                                                              <  1:1930442/100‑1 (MQ=255)
                                          tCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCgg                                                          <  2:1470950/100‑1 (MQ=255)
                                           ccGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCggg                                                         >  1:3213302/1‑100 (MQ=255)
                                           ccGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCggg                                                         <  1:3978411/100‑1 (MQ=255)
                                           ccGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCggg                                                         <  1:4222556/100‑1 (MQ=255)
                                            cGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCggg                                                         >  1:3337615/1‑99 (MQ=255)
                                              tCATTGGCGCGTATGCGGAGTCAAGTGAAACAGATAATATAAATCGCACGACAGGAAGTCGGGCGAACAGAATCAGGAGCCCAGACGGTAGTTCGGCGAc                                                      <  1:3548926/100‑1 (MQ=255)
                                              tCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGAc                                                      <  1:2108479/100‑1 (MQ=255)
                                              tCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGAc                                                      <  1:199795/100‑1 (MQ=255)
                                               cATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGAct                                                     >  2:2112198/1‑100 (MQ=255)
                                                   gACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGActc                                                    >  1:273224/1‑97 (MQ=255)
                                                   gACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCttt                                                 <  1:376372/100‑1 (MQ=255)
                                                     cGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAg                                               <  1:3838064/100‑1 (MQ=255)
                                                     cGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAg                                               <  1:976919/100‑1 (MQ=255)
                                                      gCTTATTGCGTGGCAAGTGGAACAGATAATATAAATCGCCCGGCTTGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAg                                               <  1:3592147/99‑1 (MQ=255)
                                                      gCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGt                                              <  1:1829990/100‑1 (MQ=255)
                                                       cTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGGCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGt                                              <  2:2279515/99‑1 (MQ=255)
                                                             ccGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTc                                       >  2:2318322/1‑100 (MQ=255)
                                                                aGGCAAGTGAAACAGATAATATAAATCGCCAGACATGAGGTCGGGCGAAGCGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGGAATGGTCAcg                                    <  1:4396134/100‑1 (MQ=255)
                                                                  gCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCAcg                                    >  2:1761069/1‑98 (MQ=255)
                                                                   cAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCAcgtcg                                 <  1:3185805/100‑1 (MQ=255)
                                                                     aGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTg                               <  2:550219/100‑1 (MQ=255)
                                                                     aGTGAAACAGATAATATAAATCGCCCGACATGAACTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTg                               <  1:2325487/100‑1 (MQ=255)
                                                                        gAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTCGTAATGGTCACGTCGTGacc                            >  1:1182130/1‑98 (MQ=255)
                                                                          aaCAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGt                          >  1:2145102/1‑100 (MQ=255)
                                                                          aaCAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGt                          >  1:3062752/1‑100 (MQ=255)
                                                                           aCAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCAGGTCGTGAACGTg                         <  2:2905351/100‑1 (MQ=255)
                                                                            cAGATAAGATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTgg                        <  1:4296399/100‑1 (MQ=255)
                                                                              gATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCt                      >  2:1994957/1‑100 (MQ=255)
                                                                              gATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCt                      >  1:740135/1‑100 (MQ=255)
                                                                                taataTAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCttt                    >  2:1185391/1‑100 (MQ=255)
                                                                                 gattctcttCGCCCGACTTCATGTCGGGCGATGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACCCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTc                   >  2:4245209/9‑100 (MQ=25)
                                                                                  attctcttCGCCCGACTTCATGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTcc                  >  2:305761/8‑100 (MQ=38)
                                                                                  atatAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTcc                  >  2:892702/1‑100 (MQ=255)
                                                                                   tataAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCt                 >  2:3805509/1‑100 (MQ=255)
                                                                                     taAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGa               >  1:1886611/1‑100 (MQ=255)
                                                                                     taAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGa               >  2:1146632/1‑100 (MQ=255)
                                                                                      tcttCGCCCGACTTCATGTCGGGCGAAGAGACTCAGGAGCCCAGACGGTCGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGaa              >  1:4252030/4‑100 (MQ=38)
                                                                                      aaaTCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGaa              <  2:24136/100‑1 (MQ=255)
                                                                                       cttCGCCCGACTTCATGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAAt             >  2:130116/3‑100 (MQ=255)
                                                                                         tCGCCCGACTTCATGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGc           >  2:2149636/1‑100 (MQ=255)
                                                                                          cGCCCGACTTCATGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGcc          >  2:3696466/1‑100 (MQ=255)
                                                                                          cGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGcc          >  2:2448171/1‑100 (MQ=255)
                                                                                          cGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGcc          >  1:350358/1‑100 (MQ=255)
                                                                                           gCCCGACTTCATGTCGGGTGAAGAGAATCAGGAGCCCAGACGGTAGTTCGTGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATg            >  1:1036324/1‑97 (MQ=255)
                                                                                            cccGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCg        <  2:1914806/100‑1 (MQ=255)
                                                                                            cccGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCg        <  2:1503116/100‑1 (MQ=255)
                                                                                               gACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCAc     <  2:2030507/100‑1 (MQ=255)
                                                                                                aCATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCAcc    <  2:1416403/100‑1 (MQ=255)
                                                                                                aCATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCAcc    <  2:1109377/100‑1 (MQ=255)
                                                                                                  aTGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGc  >  1:4002071/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
GAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGC  >  NC_000913/2632486‑2632683

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: