Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,907,100 T→C 6.0% E160G (GAA→GGA)  phoU ← negative regulator of PhoR/PhoB two‑component regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,907,1000TC 6.0% 208.8 / 3.6 67E160G (GAA→GGA) phoUnegative regulator of PhoR/PhoB two‑component regulator
Reads supporting (aligned to +/- strand):  ref base T (26/37);  new base C (2/2);  total (28/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.34e-01

TATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTACG  >  NC_000913/3906954‑3907248
                                                                                                                                                  |                                                                                                                                                    
gatTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCa                                                                                                                                                    <  2:340447/148‑1 (MQ=255)
 atTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCAc                                                                                                                                                   <  2:138413/149‑1 (MQ=255)
    cTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGAt                                                                                                                                                >  1:327779/1‑149 (MQ=255)
       gCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATaaa                                                                                                                                             >  4:140902/1‑149 (MQ=255)
          gCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATAc                                                                                                                                          <  2:265566/149‑1 (MQ=255)
          gCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCGTGAAGATaaaaaa                                                                                                                                          <  1:573584/149‑4 (MQ=255)
            ggTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGt                                                                                                                                        >  2:472616/1‑149 (MQ=255)
              tCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTAc                                                                                                                                      <  1:251617/149‑1 (MQ=255)
              tCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTAc                                                                                                                                      <  4:530325/149‑1 (MQ=255)
              tCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTAc                                                                                                                                      <  2:25223/149‑1 (MQ=255)
                  cAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCt                                                                                                                                  >  4:193694/1‑149 (MQ=255)
                        gTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCa                                                                                                                            >  2:388987/1‑149 (MQ=255)
                          tCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAAt                                                                                                                          <  3:140902/149‑1 (MQ=255)
                              tAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTcc                                                                                                                      <  3:158931/149‑1 (MQ=255)
                                gAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCAt                                                                                                                    >  1:347845/1‑149 (MQ=255)
                                    gcgcgcAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCcgc                                                                                                                >  3:57711/1‑149 (MQ=255)
                                    gcgcgcAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCcgc                                                                                                                <  3:532174/149‑1 (MQ=255)
                                      gcgcAGAACAGCGCAGTAAGTACGCTCGGAATCGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGACCATCcgcgc                                                                                                              >  4:361879/1‑149 (MQ=255)
                                        gcAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGa                                                                                                            >  4:155281/1‑149 (MQ=255)
                                         cAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGaa                                                                                                           <  4:454921/149‑1 (MQ=255)
                                          aGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAAc                                                                                                          <  3:5086/149‑1 (MQ=255)
                                                          taCGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCAcgtcgt                                                                                          <  1:472616/149‑1 (MQ=255)
                                                             gCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTgca                                                                                       <  3:212598/149‑1 (MQ=255)
                                                             gCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTgca                                                                                       <  3:179917/149‑1 (MQ=255)
                                                                       gTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGAt                                                                             <  3:268706/149‑1 (MQ=255)
                                                                       gTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGAt                                                                             <  2:105029/149‑1 (MQ=255)
                                                                        tACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATg                                                                            <  3:285031/149‑1 (MQ=255)
                                                                          cgcgAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGTGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGt                                                                          <  2:13929/149‑1 (MQ=255)
                                                                           cgaatttttCCATCATTTGGGTAATAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTAt                                                                        <  2:291680/142‑1 (MQ=255)
                                                                                tCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGc                                                                                                 >  3:285925/1‑120 (MQ=255)
                                                                                tCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGc                                                                                                 <  4:285925/120‑1 (MQ=255)
                                                                                  ttCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGc                                                                  <  3:361846/149‑1 (MQ=255)
                                                                                  ttCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGc                                                                  <  3:193694/149‑1 (MQ=255)
                                                                                     catcatGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGt                                                                          <  3:195684/138‑1 (MQ=255)
                                                                                     catcatGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGt                                                                          >  4:195684/1‑138 (MQ=255)
                                                                                          tGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATAAGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGAc                                                          >  1:213802/1‑149 (MQ=255)
                                                                                             aGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTcc                                                       <  1:248638/149‑1 (MQ=255)
                                                                                                tcatcaGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGa                                                    <  2:327779/149‑1 (MQ=255)
                                                                                                  atcaGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACt                                                  >  3:19070/1‑149 (MQ=255)
                                                                                                       ttGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCa                                             <  4:307319/149‑1 (MQ=255)
                                                                                                         gACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGc                                           >  4:37264/1‑149 (MQ=255)
                                                                                                            gAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAAc                                        <  4:245442/149‑1 (MQ=255)
                                                                                                            gAACAATACCTTCGTATTCCTGATCGACTTCTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGCTATGGCGGCCCAGCGACTCAAGACTTACCAGCAAc                                        <  4:350489/149‑1 (MQ=255)
                                                                                                              aCAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATGCGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACgg                                      >  2:341455/1‑149 (MQ=255)
                                                                                                                     cTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGAtgc                               >  3:164538/1‑149 (MQ=255)
                                                                                                                     cTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAAGACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGCATGGTATGACGGCCCAGCGACTCCGGACTTACCAGCAACGGCTGAtgc                               <  2:92638/149‑1 (MQ=255)
                                                                                                                          tATTCCTGATCGACTTTTTTATCTCCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATc                                                                                  >  2:228043/1‑93 (MQ=255)
                                                                                                                          tATTCCTGATCGACTTTTTTATCTCCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATc                                                                                  <  1:228043/93‑1 (MQ=255)
                                                                                                                          tATTCCTGATCGACTTTTTTATCTCCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATc                                                                                  <  3:162817/93‑1 (MQ=255)
                                                                                                                          tATTCCTGATCGACTTTTTTATCTCCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATc                                                                                  >  4:162817/1‑93 (MQ=255)
                                                                                                                           aTTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCa                                                               <  3:68586/111‑1 (MQ=255)
                                                                                                                           aTTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCa                                                               >  4:68586/1‑111 (MQ=255)
                                                                                                                           aTTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTgg                         >  2:530692/1‑149 (MQ=255)
                                                                                                                            ttCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACtt                                                 >  2:374136/1‑124 (MQ=255)
                                                                                                                            ttCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACtt                                                 <  1:374136/124‑1 (MQ=255)
                                                                                                                            ttCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTggg                        <  4:57711/149‑1 (MQ=255)
                                                                                                                               cTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGgaga                     >  1:14811/1‑149 (MQ=255)
                                                                                                                               cTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGgaga                     >  2:218157/1‑149 (MQ=255)
                                                                                                                                tGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCCGCACGTCGTGCAGCAGCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGaa                    >  2:207904/1‑149 (MQ=255)
                                                                                                                                tGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGaa                    <  1:384618/149‑1 (MQ=255)
                                                                                                                                tGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGACGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGaa                    >  1:123966/1‑149 (MQ=255)
                                                                                                                                   tCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAAttt                 <  4:253888/149‑1 (MQ=255)
                                                                                                                                     gACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTtct               <  4:262675/149‑1 (MQ=255)
                                                                                                                                        ttttttATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGga                       <  2:373487/138‑1 (MQ=255)
                                                                                                                                        ttttttATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGga                       >  1:373487/1‑138 (MQ=255)
                                                                                                                                        ttttttATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCa            >  4:228369/1‑149 (MQ=255)
                                                                                                                                         tttttATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAg           >  1:324124/1‑149 (MQ=255)
                                                                                                                                                 ttCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTAc   >  3:375884/1‑149 (MQ=255)
                                                                                                                                                 ttCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTAc   >  4:97043/1‑149 (MQ=255)
                                                                                                                                                  tCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTACg  >  4:162011/1‑149 (MQ=255)
                                                                                                                                                  tCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTACg  >  2:526150/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TATTCTGGCAGCGGTCGCCAATACGTTCGATAGAACGCGCGCAGAACAGCGCAGTAAGTACGCTCGGAATGGTACGCGAATCTTCCATCATGTAGGTCATCAGTTGACGAACAATACCTTCGTATTCCTGATCGACTTTTTTATCTTCACGATAAATACGTACCGCTTCGTCAATGTCCATCCGCGCGAACGCGTCCAGCACGTCGTGCAGCATCTGGATGGTATGACGGCCCAGCGACTCCAGACTTACCAGCAACGGCTGATGCTGCTGGGAGAATTTCTCCAGCGCAGTACG  >  NC_000913/3906954‑3907248

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: