Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,255,900 G→T 5.4% A128S (GCC→TCC)  folM → bifunctional dihydrofolate reductase/dihydromonapterin reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,255,9000GT 5.4% 240.4 / 2.2 75A128S (GCC→TCC) folMbifunctional dihydrofolate reductase/dihydromonapterin reductase
Reads supporting (aligned to +/- strand):  ref base G (38/32);  new base T (2/2);  total (41/34)
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.40e-01

TCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTACGCCGCCACCAAGGCCGCGCTGCAGAACATGGTGCTGTCATTCGCCGAAAAATATGCGCCCAAGGTGCGTGT  >  NC_002947/5255752‑5256048
                                                                                                                                                    |                                                                                                                                                    
tcACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGc                                                                                                                                                     >  2:95995/1‑150 (MQ=255)
tcACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTACACGGCGAAGCGCCCCACCACCACAACCTCGCCCTGGGCGACCTGCTGGc                                                                                                                                                     >  2:924142/1‑150 (MQ=255)
    cTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAg                                                                                                                                                 >  1:122644/1‑150 (MQ=255)
     tGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAgg                                                                                                                                                >  1:1537592/1‑150 (MQ=255)
      gCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGt                                                                                                                                               <  2:78833/150‑1 (MQ=255)
               ttATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGTGACCATGTTCCGCCTGCACGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGc                                                                                                                                      <  1:316733/150‑1 (MQ=255)
               ttATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTACCTCATCGATCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGc                                                                                                                                      <  2:1291437/150‑1 (MQ=255)
               ttATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGc                                                                                                                                      >  2:573851/1‑150 (MQ=255)
                  tGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGa                                                                                                                                   <  2:122644/150‑1 (MQ=255)
                      ggACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGAcatc                                                                                                                               >  1:994149/1‑150 (MQ=255)
                           aCGCGCCATCACCCACAACGCCTCGCTATGGCTCGACGACGCCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTACCACCTCACCCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCca                                                                                                                          <  1:997342/150‑1 (MQ=255)
                                         acaACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCGCGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCGAGGTGGAAGAGGCCGACATCATCCACATCTGCGACGaaa                                                                                                            >  2:256549/1‑150 (MQ=255)
                                         acaACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGaaa                                                                                                            >  2:448671/1‑150 (MQ=255)
                                         acaACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCATACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGaaa                                                                                                            >  2:1207815/1‑150 (MQ=255)
                                          caACGCCTCGCTATGGCTCGACGACGGCCTCGAAAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAAc                                                                                                           <  1:500434/150‑1 (MQ=255)
                                            aCGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAAcct                                                                                                         <  2:1304714/150‑1 (MQ=255)
                                                cTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTCCCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACctcctc                                                                                                     <  1:63174/150‑1 (MQ=255)
                                                    cTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTcgcgc                                                                                                 <  2:1537592/150‑1 (MQ=255)
                                                                     ccTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACatcatc                                                                                                                            >  1:700764/1‑106 (MQ=255)
                                                                     ccTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACatcatc                                                                                                                            <  2:700764/106‑1 (MQ=255)
                                                                      cTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATc                                                                               >  2:1629628/1‑150 (MQ=255)
                                                                      cTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCACCCACACCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATc                                                                               <  2:1291615/150‑1 (MQ=255)
                                                                       tCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCg                                                                              >  1:1173952/1‑150 (MQ=255)
                                                                         gACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCgg                                                                             <  1:554507/149‑1 (MQ=255)
                                                                         gACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTACCACCTCAACCTGACCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCgg                                                                             >  2:554507/1‑149 (MQ=255)
                                                                               cTCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTAcgcc                                                                      >  1:1314342/1‑150 (MQ=255)
                                                                                tCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTAcgccg                                                                     >  1:1573564/1‑150 (MQ=255)
                                                                                tCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTAcgccg                                                                     >  2:201040/1‑150 (MQ=255)
                                                                                tCGAGACCATGTTCCGCCCGCACGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTg                                                                                                                                              <  1:845303/77‑1 (MQ=255)
                                                                                tCGAGACCATGTTCCGCCCGCACGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTg                                                                                                                                              >  2:845303/1‑77 (MQ=255)
                                                                                   agaCCATGTTCCGCCTGCACGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGTTACGCCGcca                                                                  <  2:1335465/150‑1 (MQ=255)
                                                                                       cATGTTCCGCCTGCACGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTACGCCGcca                                                                  >  1:1060842/1‑146 (MQ=255)
                                                                                       cATGTTCCGCCTGCACGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTACGCCGcca                                                                  <  2:1060842/146‑1 (MQ=255)
                                                                                        aTGTTCCGCCTGCACGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTACGCCGCCACCAAg                                                             <  1:95995/150‑1 (MQ=255)
                                                                                        aTGTTCCGCCTGCACGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTACGCCGCCACCAAg                                                             <  2:1314342/150‑1 (MQ=255)
                                                                                         tGTTCCGCCTGCACGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTACGCCGCCACCAAgg                                                            <  2:994149/150‑1 (MQ=255)
                                                                                           ttCCGCCTGCACGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTACGCCGcca                                                                  <  1:1416249/142‑1 (MQ=255)
                                                                                           ttCCGCCTGCACGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTACGCCGcca                                                                  >  2:1416249/1‑142 (MQ=255)
                                                                                               gCCTGCACGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTACGCCGCCACCAAGGCCgcgc                                                      >  2:157630/1‑150 (MQ=255)
                                                                                                ccTGCACGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAAc                                                                                                           >  1:139656/1‑96 (MQ=255)
                                                                                                ccTGCACGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAAc                                                                                                           <  2:139656/96‑1 (MQ=255)
                                                                                                ccTGCACGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTACGCCGCCACCAAGGCCGCGCt                                                     <  2:986386/150‑1 (MQ=255)
                                                                                                ccTGCACGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTACGCCGCCACCAAGGCCGCGCt                                                     <  1:1449771/150‑1 (MQ=255)
                                                                                                     aCGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTACGCCGCCACCAAGGCCGCGCTGCAGa                                                >  2:283769/1‑150 (MQ=255)
                                                                                                       gTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGGAGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTACGCCGCCACCAAGGCCGCGCTGCAGAAc                                              >  2:311059/1‑150 (MQ=255)
                                                                                                             gcgcCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTACGCCGCCACCAAGGCCGCGCTGCAGAACATGGTg                                        >  1:228473/1‑150 (MQ=255)
                                                                                                                cccTACCACCTCAACCTGGCCCTGGGCNACCTGCTGGCCANGGTGGANANGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTACGCCGCCACCAAGGCCGCGCTGCAGAACATGGTGCTg                                     >  1:1004505/1‑150 (MQ=255)
                                                                                                                cccTACCACCTCAACCTGGCCCTGGGCGACCTGCTGTCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTAcgc                                                                       >  2:1200214/1‑116 (MQ=255)
                                                                                                                cccTACCACCTCAACCTGGCCCTGGGCGACCTGCTGTCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTAcgc                                                                       <  1:1200214/116‑1 (MQ=255)
                                                                                                                cccTACCACCTCAACCTGGCCCTGGGCGACCTGCTGTCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTAcgc                                                                       >  2:1200338/1‑116 (MQ=255)
                                                                                                                cccTACCACCTCAACCTGGCCCTGGGCGACCTGCTGTCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACCAAACCTCCTCTCGCGGCAGTAAAAGCCACATCGGCTAcgc                                                                       <  1:1200338/116‑1 (MQ=255)
                                                                                                                 ccTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTACGCCGCCACCAAGGCCGCGCTGCAGAACATGGTGCTGt                                    >  1:639688/1‑150 (MQ=255)
                                                                                                                  cTACCACCTCAACCTGGCCCTGGGCGACCTGCTGCCCAAGGTGAAAACGGCCGACATCATCCACATCTGCGACGAAACCGCGTCTCGCGTCGGTAAAAACTACATCGGTTTCGCCGCCACCAAGGCCGCGCTGCAGAACATGGTGCTGTc                                   >  2:184957/1‑150 (MQ=255)
                                                                                                                   tACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTACGCCGCCACCAAGGCCGCGCTGCAGAACATGGTGCTGTCa                                  <  2:1305272/150‑1 (MQ=255)
                                                                                                                   tACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTACGCCGCCACCAAGGCCGCGCTGCAGAACATGGTGCTGTCa                                  <  1:1629628/150‑1 (MQ=255)
                                                                                                                    accaccTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTACGCCGCCACCAAGGCCGCGCTGCAGAACATg                                           >  2:631869/1‑140 (MQ=255)
                                                                                                                    accaccTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTACGCCGCCACCAAGGCCGCGCTGCAGAACATg                                           <  1:631869/140‑1 (MQ=255)
                                                                                                                     ccaccTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTACGCCGCCACCAAGGCCGCGCTGCAGAACATGGTGCTGTCAtt                                <  1:1207815/150‑1 (MQ=255)
                                                                                                                      caccTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCGGTAAAAGCCACATCGGCTACGCCACCACCAAGGCCGCGCTGCAGAACATGGTGCTGTCATTc                               >  1:500603/1‑150 (MQ=255)
                                                                                                                      caccTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTACGCCGCCACCAAGGCCGCGCTGCAGAACATGGTGCTGTCATTc                               <  1:573851/150‑1 (MQ=255)
                                                                                                                      caccTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTACGCCGCCACCAAGGCCGCGCTGCAGAACATGGTGCTGTCATTc                               <  1:201040/150‑1 (MQ=255)
                                                                                                                       accTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTACGCCGCCACCAAGGCCGCGCTGCAGAACATGGTGCTGTCATTCg                              >  2:1040459/1‑150 (MQ=255)
                                                                                                                         cTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCCTCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTACGCCGCCACCAAGGCCGCGCTGCAGAACATGGTGCTGTCATTCGcc                            >  2:1369303/1‑150 (MQ=255)
                                                                                                                          tCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTACGCCGCCACCAAGGCCGCGCTGCAGAACATg                                           <  2:1116601/134‑1 (MQ=255)
                                                                                                                          tCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCAGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTACGCCGCCACCAAGGCCGCGCTGCAGAACATg                                           >  1:1116601/1‑134 (MQ=255)
                                                                                                                            aaCCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAACAGGCCGACATCATCCGCATCTGCGACGAAACCTCCTCCCGCGGCAGTAAAAGCCACATCGGCTACGCCGCCACCAAGGCCGCGCTGCAGAACATGGTGCTGTCATTCGCCGaa                         >  2:771682/1‑150 (MQ=255)
                                                                                                                                   cccTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTACGCCGCCACCAAGGCCAACCAGCGAAACATGGTTTGGTCATTCGCCGAAAAATATg                  >  2:428558/1‑150 (MQ=255)
                                                                                                                                    ccTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGc                                                                            <  2:1524736/91‑1 (MQ=255)
                                                                                                                                    ccTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGc                                                                            >  1:1524736/1‑91 (MQ=255)
                                                                                                                                    ccTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTACGCCGCCACCAAGGcc                                                          >  2:1224535/1‑109 (MQ=255)
                                                                                                                                    ccTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTACGCCGCCACCAAGGcc                                                          <  1:1224535/109‑1 (MQ=255)
                                                                                                                                    ccTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTACGCCGCCACCAAGGCCGCGCTGCAGAACATGGTGCTGTCATTCGCCGAAAAATATgc                 >  1:10805/1‑150 (MQ=255)
                                                                                                                                    ccTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTACGCCGCCACCAAGGCCGCGCTGCAGAACATGGTGCTGTCATTCGCCGAAAAAAATgc                 >  2:176509/1‑150 (MQ=255)
                                                                                                                                        ggCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTACGCCGCCACCAAGGCCGCGCTGCAGAACATGGTGCTGTCATTCGCCGAAAAATATGCGccc             >  1:693669/1‑150 (MQ=255)
                                                                                                                                             ccTGCGGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTACGCCGCCACCAAGGCCGCGCTGCAGAACATGGTGCTGTCATTCGCCGAAAAATATGCGCCCAAGGt        <  1:184957/150‑1 (MQ=255)
                                                                                                                                                   ggCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTACGCCGCCACCAAGGCCGCGCTGCAGAACATGGTGCTGTCATTCCCCGAAAAATATGCGCCCAAGGTGCgtgt  <  1:1040459/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCACCTGCACGATCGTTATGACGGACTACGCGCCATCATCCACAACGCCTCGCTATGGCTCGACGACGGCCTCGACAACCTCGAGACCATGTTCCGCCTGCACGTCGAAGCGCCCTACCACCTCAACCTGGCCCTGGGCGACCTGCTGGCCAAGGTGGAAAAGGCCGACATCATCCACATCTGCGACGAAACCTCCTCGCGCGGCAGTAAAAGCCACATCGGCTACGCCGCCACCAAGGCCGCGCTGCAGAACATGGTGCTGTCATTCGCCGAAAAATATGCGCCCAAGGTGCGTGT  >  NC_002947/5255752‑5256048

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: