Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,570,651 A→G 8.5% H465H (CAT→CAC gadB ← glutamate decarboxylase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,570,6510AG 8.5% 166.2 / 5.2 70H465H (CAT→CACgadBglutamate decarboxylase B
Reads supporting (aligned to +/- strand):  ref base A (27/35);  new base G (2/4);  total (31/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.54e-01

TTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCATCAC  >  NC_000913/1570510‑1570788
                                                                                                                                             |                                                                                                                                         
ttCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAAAGTTATCAGGTaa                                                                                                                                          >  2:329735/1‑142 (MQ=255)
   cTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGACACGGGTCAGTTTCGGGACGCCGCTACCGTTAAACGTTGTCAGGCTTGtt                                                                                                                                       >  1:172096/1‑143 (MQ=255)
    tGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTGAACGTTATCAGGTATGctt                                                                                                                                      >  2:594730/1‑143 (MQ=255)
       cGGACCCGCACCTGTTTTGTTCTGAAGCTTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGTACAGACTTTGGAAATTGTCCCGAAACGGGCTCCTGGTAGTACTCCTTTACCGTTAAACGTTATCAGGTATGTTTaaa                                                                                                                                   <  2:573221/143‑1 (MQ=255)
          aaCCGCTCCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGct                                                                                                                                <  2:1129248/143‑1 (MQ=255)
           aCCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGctg                                                                                                                               <  1:1252410/143‑1 (MQ=255)
              gCACCTGTTTTGTTCTGAAGCATATACAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACCTTGAAAATTGTCCCGAAACGAGCTCGTTTCGGGACACAGTAACCCGTAGACGGTAACAGGTATGGCAAAAGctggac                                                                                                                            >  2:8082/1‑140 (MQ=255)
                 ccTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGctgttctgt                                                                                                                         >  2:464358/1‑143 (MQ=255)
                  cTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGctgttctgtt                                                                                                                        <  1:941691/143‑1 (MQ=255)
                   tGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTg                                                                                                                       >  1:1123168/1‑143 (MQ=255)
                   tGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTg                                                                                                                       >  1:43570/1‑143 (MQ=255)
                    gTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACGCCGCTACCGTTCAACGTTCTCAGGCAGGTGTAAAGCTGGTCCCTTgg                                                                                                                      >  2:331094/1‑143 (MQ=255)
                     ttttGTTCTGAAGCTTATTCAGAACATTATTTTCCGTTGCTAATGCCAGTTAACATTCTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACAGCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTggg                                                                                                                     <  1:12527/143‑1 (MQ=255)
                     ttttGTTCGGAACCTTATTGAGAACTATCTTTGCCGTTCCCAATGCGAGTGAACAGACTTTGGAGGTTTGCCCGAAACGGGTTCGTTTCGTGACACCGTTATCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTggg                                                                                                                     <  2:343691/143‑1 (MQ=255)
                        tGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGGACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACGGTTACCGTTAAACGTTATCAGGTATGGTGACAGCGGTTCTGTTGGGCaa                                                                                                                  >  2:923021/1‑143 (MQ=255)
                          ttCTGAAGCGTATTCCGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATa                                                                                                                <  1:743509/143‑1 (MQ=255)
                             tGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATAccc                                                                                                             >  1:501470/1‑143 (MQ=255)
                              gAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTACCAGGTATGTTTAACGCTGTTCTGTTGGGCAATACCCt                                                                                                            >  2:1003520/1‑143 (MQ=255)
                               aggCGTACTCAGAAGCAGAGTTTCGGTTGCTCGCGCCACTGAGCGGGCGTTGGAAATTGTCCCGAAGCGGGTGCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTg                                                                                                           <  2:959308/141‑1 (MQ=255)
                               aaGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGCAATTGTCCCGAAACGGGTTCGTTTCGCGACACCGTTACCGTTAAACGTTATCAGGGATGTTTAAAGCTGTTCTGTTGGGCAATACCCTg                                                                                                           >  1:174901/1‑143 (MQ=255)
                                 gCGTATTCAGAACAAGGTTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCa                                                                                                         <  2:517062/143‑1 (MQ=255)
                                  cGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATAACCTGCAg                                                                                                        >  1:851952/1‑143 (MQ=255)
                                      ttCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACATAGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACCCTGTTACCGTTCAACGCTATAAGGTATGTTTAACGCCGTTCCGCTGGGCAATAGCATGCGGTGTc                                                                                                    >  2:896966/1‑143 (MQ=255)
                                       tCAGAACAATATTTTCCGTTGCGAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCg                                                                                                   <  1:1279899/143‑1 (MQ=255)
                                        cAGAAAAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCgg                                                                                                  <  2:1124334/143‑1 (MQ=255)
                                         aGAACGATATTTTGCGTGGCTAGGGGCCGTGGACAGAGTCGGGGCATTGTCCCGAAACGGGTGCGTTTCGGGACACCGGGACCGTTATACGTTTTCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCGGCAGTTTCggg                                                                                                 <  2:1243377/143‑1 (MQ=255)
                                             cagtagtttCCGTTGCGAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGAt                                                                                             <  2:1050912/137‑1 (MQ=255)
                                             cAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGAt                                                                                             >  2:526114/1‑143 (MQ=255)
                                             cAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGAt                                                                                             >  2:1201540/1‑143 (MQ=255)
                                                  ttttCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAATCAAATTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTg                                                                                        <  1:238060/143‑1 (MQ=255)
                                                   tgtCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTga                                                                                       <  1:267388/141‑1 (MQ=255)
                                                    ttCCGTTGCTAACGCCAGTGAACAGACTGTGGCAATTGTCCCCAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTgag                                                                                      <  1:329735/143‑1 (MQ=255)
                                                     tCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTgaga                                                                                     <  1:464358/143‑1 (MQ=255)
                                                           aacataacgttgtaaaaaccgaatgcccagcctttaaaaaaacagctgggcattcggtttttacaacgttatGTTATCAGGTGTGTTTAAAGCTGTTCTGCTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTc                                                                               <  2:796724‑M1/71‑1 (MQ=255)
                                                           aacataacgttgtaaaaaccgaatgcccagcctttaaaaaaacagctgggcattcggtttttacaacgttatGTTATCAGGTGTGTTTAAAGCTGTTCTGCTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTc                                                                               <  2:169504‑M1/71‑1 (MQ=255)
                                                             tAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACGGGTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTGCAg                                                                             >  1:734324/1‑143 (MQ=255)
                                                                     gTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTGAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGGAATACCCTGCAGTTTCGAGTCATCGCTGAGATATTTCACGGAGGcct                                                                     >  2:149977/1‑141 (MQ=255)
                                                                      tGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTCAGCTCTTTCCGGGCGGCttt                                                                    >  1:372828/1‑143 (MQ=37)
                                                                         aCAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTa                                                                 <  2:919465/143‑1 (MQ=255)
                                                                           aGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTCTGTTTAAAGCTGTTCTGTTAGGCACTACCCTGCAGTTCCGGGTGATGACTGACATAATTCAGGGAGGCtgggaagc                                                               >  2:328610/1‑136 (MQ=255)
                                                                            gACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTc                                                              <  1:847141/143‑1 (MQ=255)
                                                                             acgttGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCt                                                             <  1:528658/140‑1 (MQ=255)
                                                                                  ggAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAg                                                        >  1:332780/1‑143 (MQ=255)
                                                                                   gAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGc                                                       <  2:100133/143‑1 (MQ=255)
                                                                                    cccagcctttaaaaaaacagctgggcattcggtttttacaacgttatGTTATCAGGTGTGTTTAAAGCTGTTCTGCTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGGTATTTCAGGGAGGCTTTGTAGTCTTCCAGCa                                                      <  1:647089‑M1/96‑1 (MQ=255)
                                                                                      aTTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTGTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAAc                                                    <  1:29151/143‑1 (MQ=255)
                                                                                       ttGTCCGGAAGGGGGTCCGTGTCGGGAGCGCGTTACCGTTAAACCGGATCAGGTATGTTTAAAGCTGTTCGGGTGGGCAATACCCTGCAGTTTCGGGTGGTCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACa                                                   <  2:1261838/143‑1 (MQ=14)
                                                                                        tGTCCCGAAACGAACCCGTTTCGGGACACCGTTACCGTTAAACGTGATCATGTATGTTTAAAGCTGTTCTGTTGTGCAATACCCCGCAGTTTCGGCCGATCGCTGAGACATTTGAGGGAGGCTTGGTAGTCTTCCCGCacctg                                                  >  2:547097/1‑139 (MQ=255)
                                                                                            ttaaaaaaacagctgggcattcggtttttacaacgttatGTTATCAGGTGTGTTTAAAGCTGTTCTGCTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCa                                              >  1:796724‑M1/40‑143 (MQ=255)
                                                                                            ttaaaaaaacagctgggcattcggtttttacaacgttatGTTATCAGGTGTGTTTAAAGCTGTTCTGCTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCa                                              >  1:169504‑M1/40‑143 (MQ=255)
                                                                                            gccgatcGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCa                                              <  2:255757/137‑1 (MQ=255)
                                                                                            ccGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCACTACCCTGCAGTGTCGGGTGATCGCTGAGATATTTCAGGGAGGCGTTGTAGTCTTCCAGCAACAGGTCa                                              >  2:1166810/1‑143 (MQ=255)
                                                                                            ccGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCa                                              <  1:1108882/143‑1 (MQ=255)
                                                                                            ccGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGATGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCCACAGTTCa                                              >  1:30114/1‑143 (MQ=255)
                                                                                              ggaaggggTGCGCTGCGGGGCACCGTGACCGTTAAACGTGAGCAGGTGTGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGc                                            <  1:1188337/138‑1 (MQ=14)
                                                                                                         gTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAGAGTCCAgtt                                 >  1:1035854/1‑143 (MQ=255)
                                                                                                          tttCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGCTCGCTGAGATCTGTGAGGGGGGAGCCCGCGCCGTCCAGCACCAGATCAGCAAAGTCCATTTc                                >  2:1240135/1‑143 (MQ=17)
                                                                                                            tCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGa                              <  1:997866/143‑1 (MQ=255)
                                                                                                             ctgtacacCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGaa                             <  2:548497/139‑1 (MQ=255)
                                                                                                              gggACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTCCTGTTGGGCAATACTCTTCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAg                            <  1:174338/143‑1 (MQ=255)
                                                                                                                    ccggtacccttaAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGTGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTGCAGCAAAGTCCATTTCGAAGCCGcga                      <  2:1281205/143‑1 (MQ=255)
                                                                                                                      gttaccgttaAACGTTATCAGATATGTTTAAAACTGATCTGTTGTGCAATACACTGCAGTATCAAGAGATCGATGAGATATTTCAACAAAGCTTTGTAGTCTTCCAGCAACAGATCAGCAAAGTACAATTCGAATCAGcgaag                    >  2:238863/1‑141 (MQ=9)
                                                                                                                       ttaccgttaAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGcgacga                   <  1:713635/143‑1 (MQ=255)
                                                                                                                           cgttaAACGTTATCAGGTATGTTTAAAGCTGTTTTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGCcaca               >  1:972357/1‑143 (MQ=255)
                                                                                                                                   gTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGCAGCCGCGACGCCACCTAATgcgg       >  1:128837/1‑142 (MQ=255)
                                                                                                                                   gTTATCAGGTATGTTTAAAGCGGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCGTTTCGAAGCCGCGACGACGCCTTATgcgc       >  1:759929/1‑143 (MQ=11)
                                                                                                                                     ttttgggTATGTTTAAGGGTGTTCTGTTGGGCAATACGCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCAt     <  2:1097910/138‑1 (MQ=255)
                                                                                                                                     tATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCAt     <  1:1201540/143‑1 (MQ=35)
                                                                                                                                     tATCAGGGATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCAt     <  1:526114/143‑1 (MQ=21)
                                                                                                                                      aTCAGGTATGTTTAAAGCTGTTCTGTTAGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCATc    <  1:358839/143‑1 (MQ=25)
                                                                                                                                        cAGGTATGTTTAAAGCTGTTCTGTTGGTCAATCCCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGTGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCATcac  <  2:789745/143‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                         
TTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCATCAC  >  NC_000913/1570510‑1570788

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