Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,155,565 G→A 5.1% A189T (GCA→ACA)  argC → N‑acetyl‑gamma‑glutamylphosphate reductase, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,155,5650GA 5.1% 254.0 / 2.1 78A189T (GCA→ACA) argCN‑acetyl‑gamma‑glutamylphosphate reductase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base G (39/35);  new base A (2/2);  total (41/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.74e-01

CGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATCGCCATCAACCAGAGATCGCCACACACCTCGGTGCTGACGTTATCTTCACCCCACATCTGGGCAATTTCCCGCGCG  >  NC_000913/4155417‑4155712
                                                                                                                                                    |                                                                                                                                                   
cGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGCCCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTg                                                                                                                                                     >  1:250571/1‑149 (MQ=255)
   taataaATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAg                                                                                                                                                  <  7:403456/149‑1 (MQ=255)
   taataaATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAg                                                                                                                                                  >  2:148280/1‑149 (MQ=255)
     ataaATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTg                                                                                                                                                     >  6:66845/1‑144 (MQ=255)
     ataaATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTg                                                                                                                                                     <  5:66845/144‑1 (MQ=255)
     ataaATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAgg                                                                                                                                                 <  1:148280/148‑1 (MQ=255)
        aaTTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGt                                                                                                                                             <  6:47998/149‑1 (MQ=255)
         aTTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGt                                                                                                                                             >  3:118743/1‑148 (MQ=255)
            aaaaGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGTTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTaaa                                                                                                                                          <  4:181292/148‑1 (MQ=255)
            aaaaGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAg                                                                                                                                         >  8:15322/1‑149 (MQ=255)
             aaaGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAgg                                                                                                                                                 <  8:290035/140‑1 (MQ=255)
             aaaGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAgg                                                                                                                                                 >  7:290035/1‑140 (MQ=255)
               agaagCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCgg                                                                                                                                      >  2:291136/1‑149 (MQ=255)
                gaagCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGc                                                                                                                                     >  2:322036/1‑149 (MQ=255)
                 aagCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGc                                                                                                                                     >  8:153396/1‑148 (MQ=255)
                   gCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCAt                                                                                                                                  <  2:353679/149‑1 (MQ=255)
                   gCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCAt                                                                                                                                  >  1:61694/1‑149 (MQ=255)
                   gCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCAt                                                                                                                                  <  7:153396/149‑1 (MQ=255)
                     gaatgggATGGCGGTGGCCGGCCGATCTCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGTTGCCGATCTTCTTCACCTCAATGAGTTGCTTGTGAATAACGGCAACACCGGCGCGAGCGGTGCAGTGCGTAACGCGGCCAttt                                                                                                                                <  2:250571/143‑1 (MQ=255)
                      aaTTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTc                                                                                                                               <  7:372566/149‑1 (MQ=255)
                       aTTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTc                                                                                                                               >  3:207102/1‑148 (MQ=255)
                          tGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCaaa                                                                                                                            <  1:291136/148‑1 (MQ=255)
                             ttGCGGTGCCGGGCTGTTATCCGCCGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCTGTGATCAACGCCACCAGCGGCGTGCGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGc                                                                                                                        >  3:40466/1‑149 (MQ=255)
                             ttGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGc                                                                                                                        >  4:298852/1‑149 (MQ=255)
                                cGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCttt                                                                                                                     <  1:322036/149‑1 (MQ=255)
                                  gggCCGGTCTGTTATACGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTTATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTtg                                                                                                                   <  4:215315/147‑1 (MQ=255)
                                  gTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTtg                                                                                                                   <  2:380177/149‑1 (MQ=255)
                                   tGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTtgt                                                                                                                  <  4:118743/149‑1 (MQ=255)
                                             ttATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGc                                                                                                         <  7:118082/148‑1 (MQ=255)
                                             ttATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGc                                                                                                         <  7:112263/148‑1 (MQ=255)
                                               aTCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTg                                                                                                      >  2:118077/1‑149 (MQ=255)
                                               aTCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATCCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGACACCAGCGTCGTGAGAGCTGCAGGGCGTAAAGCGGCCATTTCAAACACCTTTTGTGAAGTTAGccag                                                                                                      >  2:55707/1‑147 (MQ=255)
                                                   gACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAACGCGGCCATTTCAAACAGCATTTGTGAAGTTAGCCTGCAAc                                                                                                  >  6:13398/1‑149 (MQ=255)
                                                        cggcACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTAt                                                                                             <  2:61694/149‑1 (MQ=255)
                                                          gcACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATgg                                                                                           <  1:317280/149‑1 (MQ=255)
                                                             caGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGt                                                                                        >  1:346364/1‑149 (MQ=255)
                                                               gCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCt                                                                                      <  3:367256/149‑1 (MQ=255)
                                                                cTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCtt                                                                                     <  2:369132/149‑1 (MQ=255)
                                                                  ggCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCaa                                                                                                   <  3:200140/133‑1 (MQ=255)
                                                                  ggCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCaa                                                                                                   >  4:200140/1‑133 (MQ=255)
                                                                   gcgcTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTAc                                                                                  <  3:263519/149‑1 (MQ=255)
                                                                   gcgcTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTAc                                                                                  <  4:207102/149‑1 (MQ=255)
                                                                   gcgcTGAAACAGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGGGATCAACGCCTCCATCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCGATTTCAAACAGCTTTTGTTAAGTTAGCCTGCAACCGTATGGCGTCTTTa                                                                                   <  2:142334/148‑1 (MQ=255)
                                                                        gAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCACACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTAATc                                                                             >  1:152184/1‑149 (MQ=255)
                                                                            ccGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATCGCCa                                                                         >  6:162483/1‑149 (MQ=255)
                                                                               ttgattgaTGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGaa                                                                                                               <  6:315712/108‑1 (MQ=255)
                                                                               ttgattgaTGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGaa                                                                                                               >  5:315712/1‑108 (MQ=255)
                                                                                 gattgatGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATCGCCATCAAc                                                                    >  7:212355/1‑149 (MQ=255)
                                                                                  attgatGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATCGCCATCAAc                                                                    >  1:118359/1‑148 (MQ=255)
                                                                                   ttgatGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATAGCGTCTTTACTCATCGCCATCAAc                                                                    >  3:22043/1‑147 (MQ=255)
                                                                                      atGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATCGCCATCAACCagag                                                               <  5:162483/149‑1 (MQ=255)
                                                                                        gCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATCGCCATCAACCa                                                                  >  5:97171/1‑144 (MQ=255)
                                                                                         ccGATCTTCTTTACCTCAATCAGTGGCCGGAGATCAACGCCAACAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGAGTCTTTACTCATCGCCATCAACGAGAGAt                                                             <  6:15024/148‑1 (MQ=255)
                                                                                           gATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATCGCCATCAACCAGAGATCGc                                                          >  8:204876/1‑149 (MQ=255)
                                                                                           gATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATCGCCATCAACCAGAGATCGc                                                          >  6:302221/1‑149 (MQ=255)
                                                                                                    gACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATCGCCATCAACCAGAGATCGCCACACACCt                                                 >  8:380596/1‑149 (MQ=255)
                                                                                                      ccTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTTAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATCGCCATCAACCAGAGATCGCCACACACCTCg                                               >  2:95603/1‑149 (MQ=255)
                                                                                                      ccTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATCGCCATCAACCAGAGATCGc                                                          <  2:309408/138‑1 (MQ=255)
                                                                                                      ccTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATCGCCATCAACCAGAGATCGc                                                          >  1:309408/1‑138 (MQ=255)
                                                                                                        tctaTCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATCGCCATCAACCAGAGATCGCCACACACCTCGGt                                             <  6:73902/146‑1 (MQ=255)
                                                                                                        tCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATCGCCATCAACCAGAGATCGCCACACACCTCGGt                                             >  7:165719/1‑149 (MQ=255)
                                                                                                         cAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTACAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATCGCCATCAACCag                                                                 <  5:242062/128‑1 (MQ=255)
                                                                                                         cAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTACAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATCGCCATCAACCag                                                                 >  6:242062/1‑128 (MQ=255)
                                                                                                          aaTCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATCGCCATCAACCAGAGATCGCCACACACCTCGGTGc                                           >  7:342690/1‑149 (MQ=255)
                                                                                                                    cGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATCGCCATCAACCAGAGATCGCCACACACCTCGGTGCTGACGTTATc                                 <  3:307591/149‑1 (MQ=255)
                                                                                                                        gATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATCGCCATCAACCAGAGATCGCCACACACCTCGGTGCTGACGTTATCTTCa                             >  1:284450/1‑149 (MQ=255)
                                                                                                                         aTCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATCGCCATCAACCAGAGATCGCCACACACCTCGGTGCTGACGTTATCTTCAc                            <  2:346364/149‑1 (MQ=255)
                                                                                                                          tCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATCGCCATCAACCAGAGATCGCCACACACCTCGGTGCTGACGTTATCTTCAcc                           >  3:310898/1‑149 (MQ=255)
                                                                                                                                 caccaGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATCGCCATCAACCagag                                                               >  6:115771/1‑106 (MQ=255)
                                                                                                                                 caccaGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATCGCCATCAACCAGAGATCGCCACACACCTCGGTGCTGACGTTATCTTCACCCCACATc                    >  4:68639/1‑149 (MQ=255)
                                                                                                                                 caccaGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATCGCCATCAACCAGAGATCGCCACACACCTCGGTGCTGACGTTATCTTCACCCCACATc                    >  4:73646/1‑149 (MQ=255)
                                                                                                                                 caccaGCGGCGTGAGCGGTACAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGaa                                                                                                               >  4:116521/1‑58 (MQ=255)
                                                                                                                                 caccaGCGGCGTGAGCGGTACAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGaa                                                                                                               <  3:116521/58‑1 (MQ=255)
                                                                                                                                  accaGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATCGCCATCAACCAGAGATCGCCACACACCTCGGTGCTGACGTTATCTTCACCCCACATCt                   >  4:209269/1‑149 (MQ=255)
                                                                                                                                    caGCGGCGTGATCTGTGGAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTACCCTGCAACCGTATGGCGTCTTTACTCATCGCCATCAACCAGAGATc                                                            <  8:86897/106‑1 (MQ=255)
                                                                                                                                    caGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTCAAGTTAGCCTGCACCCGTATGGCGTCTTTACTCATCGCCATCAACCAGAGATc                                                            >  7:86897/1‑106 (MQ=255)
                                                                                                                                      gcggcgTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCACCCGTATGGCGTCTTTACTCATCGCCATCAACCAGAGATCGCCACCCACCTCGGTGCTGACGTTATCTTCACCCCACATCTggg                >  6:185310/1‑148 (MQ=255)
                                                                                                                                          cgTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATCGCCATCAACCAGAGATCGCCACACACCTCGGTGCTGACGTTATCTTCACCCCACATCTGGGCAAtt           <  7:380596/149‑1 (MQ=255)
                                                                                                                                             gAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATCGCCATCAACCAGAGATCGCCACACACCTCGGTGCTGACGtt                                    <  6:95207/121‑1 (MQ=255)
                                                                                                                                                 ggTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATCGCCATCAACCAGAGATCGCCACACACCTCGGTGCTGACGTTATCTTCACCCCACATCTGGGCAATTTCCcgcg    >  4:277299/1‑149 (MQ=255)
                                                                                                                                                   tGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATCGCCATCAACCAGAGATCGCCACACACCTCGGTGCTGACGTTATCTTCACCCCACATCTGGGCAATTTCCcgcgcg  <  5:302221/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGGTAATAAATTAAAAGAAGCGAATTTGATTGCGGTGCCGGGCTGTTATCCGACGGCGGCACAGCTGGCGCTGAAACCGTTGATTGATGCCGATCTTCTTGACCTCAATCAGTGGCCGGTGATCAACGCCACCAGCGGCGTGAGCGGTGCAGGGCGTAAAGCGGCCATTTCAAACAGCTTTTGTGAAGTTAGCCTGCAACCGTATGGCGTCTTTACTCATCGCCATCAACCAGAGATCGCCACACACCTCGGTGCTGACGTTATCTTCACCCCACATCTGGGCAATTTCCCGCGCG  >  NC_000913/4155417‑4155712

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: