Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,570,669 T→C 8.8% Q459Q (CAA→CAG gadB ← glutamate decarboxylase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,570,6690TC 8.8% 145.8 / 4.9 64Q459Q (CAA→CAGgadBglutamate decarboxylase B
Reads supporting (aligned to +/- strand):  ref base T (25/31);  new base C (3/3);  total (29/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.75e-01

CCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGC  >  NC_000913/1570527‑1570811
                                                                                                                                              |                                                                                                                                              
ccTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGctgttctgt                                                                                                                                                >  2:464358/1‑143 (MQ=255)
 cTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGctgttctgtt                                                                                                                                               <  1:941691/143‑1 (MQ=255)
  tGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTg                                                                                                                                              >  1:43570/1‑143 (MQ=255)
  tGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTg                                                                                                                                              >  1:1123168/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:169504‑M1/71‑1 (MQ=255)
                                          aacataacgttgtaaaaaccgaatgcccagcctttaaaaaaacagctgggcattcggtttttacaacgttatGTTATCAGGTGTGTTTAAAGCTGTTCTGCTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTc                                                                                                      <  2:796724‑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)
                                                                                                                                 tAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGAGCATCACCACGATGTCg               >  2:161662/1‑143 (MQ=17)
                                                                                                                                    aGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGAGATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCATCACCACGCTGTCGGTg            >  2:123864/1‑143 (MQ=17)
                                                                                                                                        gttctgttGGGGAATACGCTGCAGTTTCGGGTGAGCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCtt        <  2:377746/143‑1 (MQ=12)
                                                                                                                                        ggtctgttGGGCACTCCCCTGCAGTTGCGGGTGATCGCTGAGATATTGCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCtt        <  1:547097/141‑1 (MQ=255)
                                                                                                                                         ttctgttGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTGCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGCGACAGAATGCGCAGCACCACGATATAGGTGGCGTc       >  2:125127/1‑143 (MQ=11)
                                                                                                                                              ttGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCAccgc  >  1:824536/1‑143 (MQ=11)
                                                                                                                                              |                                                                                                                                              
CCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTTCGTTTCGGGACACCGTTACCGTTAAACGTTATCAGGTATGTTTAAAGCTGTTCTGTTGGGCAATACCCTGCAGTTTCGGGTGATCGCTGAGATATTTCAGGGAGGCTTTGTAGTCTTCCAGCAACAGTTCAGCAAAGTCCATTTCGAAGCCGCGACGACACATAATGCGCATCACCACGATGTCGGTGGCTTCACCGC  >  NC_000913/1570527‑1570811

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: