Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,057,709 G→A 5.5% R82R (CGC→CGT torR ← response regulator in two‑component regulatory system with TorS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,057,7090GA 5.5% 235.0 / 2.4 73R82R (CGC→CGTtorRresponse regulator in two‑component regulatory system with TorS
Reads supporting (aligned to +/- strand):  ref base G (37/32);  new base A (2/2);  total (39/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 = 4.01e-01

TGAGTGTGCGGTTGAGCTTGTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATCAGGCCATTTTCATCGGGTAAGTTGATATCCAGCAGAATTAAATCTACCGACTGATTCTGCATAATTTCCCGTAGCCCAGCACCGCT  >  NC_000913/1057563‑1057855
                                                                                                                                                  |                                                                                                                                                  
tGAGTGTGCGGTTGAGCTTGTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGc                                                                                                                                                   >  5:207011/1‑148 (MQ=255)
 gAGTGTGCGGTTGAGCTTGTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCg                                                                                                                                                  >  2:190936/1‑148 (MQ=255)
    tgtgCGGTTGAGCTTGTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCTGATCGCTGCGTcc                                                                                                                                               <  1:216310/148‑1 (MQ=255)
     gtgCGGTTGAGCTTGTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTTGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCg                                                                                                                                              >  1:42046/1‑148 (MQ=255)
     gtgCGGTTGAGCTTGTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCg                                                                                                                                              >  4:375199/1‑148 (MQ=255)
       gCGGTTGAGCTTGTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAAACCGACCGCTGCGTCCgg                                                                                                                                             >  8:158458/1‑147 (MQ=255)
             gAGCTTGTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCa                                                                                                                                       <  3:375199/147‑1 (MQ=255)
               gCTTGTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGaa                                                                                                                                    <  4:199024/148‑1 (MQ=255)
                 ttGTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAAcc                                                                                                                                        >  1:252016/1‑142 (MQ=255)
                 ttGTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAAcc                                                                                                                                        <  2:252016/142‑1 (MQ=255)
                              tCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGtt                                                                                                                      >  8:137556/1‑147 (MQ=255)
                              tCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTg                                                                                                                     <  6:384995/148‑1 (MQ=255)
                                  ttCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCg                                                                                                                 >  1:65988/1‑148 (MQ=255)
                                      ccAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTtct                                                                                                             >  2:61421/1‑148 (MQ=255)
                                            cAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAgg                                                                                                       >  1:220241/1‑148 (MQ=255)
                                             aGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAggg                                                                                                      <  1:348922/148‑1 (MQ=255)
                                                  tttCACCCGTACTCCCAGTTCGCGCAGTTCCAGCGGTGTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCACAATAATCCCCACCGTTGAGCGTTCTCGCAGGGcgcg                                                                                                  <  1:326663/147‑1 (MQ=255)
                                                   ttCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCgg                                                                                                 <  4:221958/147‑1 (MQ=255)
                                                    tCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGc                                                                                                          <  6:416324/137‑1 (MQ=255)
                                                    tCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGc                                                                                                          >  5:416324/1‑137 (MQ=255)
                                                          gTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACAt                                                                                         <  1:190936/148‑1 (MQ=255)
                                                           tactacCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACAt                                                                                         >  8:395547/1‑147 (MQ=255)
                                                           tactacCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATc                                                                                        <  5:3281/148‑1 (MQ=255)
                                                            actacCAGTTCGCGCAGTTCCAGCGTTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAAc                                                                                           <  1:66512/144‑1 (MQ=255)
                                                            actacCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCgg                                                                                                 >  4:304437/1‑138 (MQ=255)
                                                            actacCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCgg                                                                                                 <  3:304437/138‑1 (MQ=255)
                                                            actacCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAAc                                                                                           >  2:66512/1‑144 (MQ=255)
                                                              tacCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATCAgg                                                                                     >  4:79469/1‑148 (MQ=255)
                                                                ccAGTTCGCGCAGTTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTGCTCGCAGGGCGCGGGTTAACATCAgg                                                                                     <  7:137556/146‑1 (MQ=255)
                                                                 cAGTTCGCGCAGTTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATCAGGCCa                                                                                  >  8:95732/1‑148 (MQ=255)
                                                                     tCGCGCAGTTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATcc                                                                                                                              >  8:269474/1‑100 (MQ=255)
                                                                     tCGCGCAGTTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATcc                                                                                                                              <  7:269474/100‑1 (MQ=255)
                                                                      cgcgCAGTTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATCAGGCCAtttt                                                                              <  2:220241/147‑1 (MQ=255)
                                                                      cgcgCAGTTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATCAGGCCAtttt                                                                              <  7:395547/147‑1 (MQ=255)
                                                                        cgcAGTTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATCTGGCCATTTTCAt                                                                           <  6:207011/148‑1 (MQ=255)
                                                                          cAGTTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATCAGGCCATTTTCATCg                                                                         >  5:249108/1‑148 (MQ=255)
                                                                          cAGCTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATCAGGCCATTTTCATCg                                                                         >  2:242881/1‑148 (MQ=255)
                                                                           aGTTCCAGCGGTTTGGGGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATCAGGCCATTTTCATCgg                                                                        <  8:406758/148‑1 (MQ=255)
                                                                                cAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATCAGGCCATTTTCATCGGGTAAg                                                                   >  2:203388/1‑148 (MQ=255)
                                                                                 aGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATCAGGCCATTTTCATCGGGTAAGt                                                                  >  7:237113/1‑148 (MQ=255)
                                                                                     gTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATCAGGCCATTTTCATCGGGTAAGTTGa                                                               >  2:417820/1‑147 (MQ=255)
                                                                                            gACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATCAGGCCATTTTCATCGGGTAAGTTGATATCcag                                                        >  4:353830/1‑147 (MQ=255)
                                                                                            gACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATCAGGCCATTTTCATCGGGTAAGTTGATATCcag                                                        >  8:145674/1‑147 (MQ=255)
                                                                                               gTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATCAGGCCATTTTCATCGGGTAAGTTGATATCCAGCAGa                                                    <  8:237113/148‑1 (MQ=255)
                                                                                                   tCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATCAGGCCATTTTCATCGGGTAAGTTGATATCCAGCAGAATTa                                                <  2:42046/148‑1 (MQ=255)
                                                                                                    cGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATCAGGCCATTTTCATCGGGTAAGTTGATATCCAGCAGAATTa                                                >  4:254784/1‑147 (MQ=255)
                                                                                                      tCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATCAGGCCATTTTCATCGGGTAAGTTGATATCCAGCAGAATTAAAt                                             <  1:61421/148‑1 (MQ=255)
                                                                                                            cccATTTCCAGCCCAACAATACGGTCAATCCGATCTCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATCAGGCAATTTTCATCGGGTAAGTTGATATCCAGCAGAATTAAATCTagc                                        >  6:132534/1‑145 (MQ=255)
                                                                                                            cccATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAAc                                                                                           <  7:85294/96‑1 (MQ=255)
                                                                                                            cccATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAAc                                                                                           >  8:85294/1‑96 (MQ=255)
                                                                                                                tttCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATCAGGCCATTTTCATCGGGTAAGTTGATATCCAGCAGAATTAAATCTACCGACt                                    >  3:441431/1‑147 (MQ=255)
                                                                                                                 ttCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTg                                                                                                                     <  5:252646/65‑1 (MQ=255)
                                                                                                                 ttCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTg                                                                                                                     >  6:252646/1‑65 (MQ=255)
                                                                                                                   ccAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATCAGGCCATTTTCATCGGGTAAGTTGATATCCAGCAGAATTAAATCTACCGACTGAt                                 <  7:145674/147‑1 (MQ=255)
                                                                                                                           acaaTACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATCAGGCCATTTTCATCGGGTAAGTTGCTATCCAGCAGAATTAAATCTACCGACTGATTCTGCATaa                        >  3:97275/1‑148 (MQ=255)
                                                                                                                              aTACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATCAGGCCATTTTCATCGGGTAAGTTGATATCCAGCAGAATTAAATCTACCGACTGATTCTGCATAAttt                     <  7:95732/148‑1 (MQ=255)
                                                                                                                                aCGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATCAGGCCATTTTCATCGGGTAAGTTGATATCAAGCATAATTAAATCTACCGACTGATTCTGCATAATTTcc                   >  2:339123/1‑148 (MQ=255)
                                                                                                                                 cGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATCAGGCCATTTTCATCGGGTAAGTTGATATCCAGCAGAATTAAATCTACCGACTGATTCTGCATAATTTccc                  <  6:249108/148‑1 (MQ=255)
                                                                                                                                   gTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATCAGGCCATTTTCATCGGGTAAGTTGATATCCAGCAGAATTAAATCTACCGACTGATTCTGCATAATTTc                    >  1:151711/1‑144 (MQ=255)
                                                                                                                                   gTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATCAGGCCATTTTCATCGGGTAAGTTGATATCCAGCAGAATTAAATCTACCGACTGATTCTGCATAATTTc                    <  2:151711/144‑1 (MQ=255)
                                                                                                                                   gTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATCAGGCCATTTTCATCGGGTAAGTTGATATCCAGCAGAATTAAATCTACCGACTGATTCTGCATAATTTCCCGt                <  1:242881/148‑1 (MQ=255)
                                                                                                                                   gTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATCAGGCCATTTTCATCGGGTAAGTTGATATCCAGCAGAATTAAATCTACCGACTGATTCTGCATAATTTCCCGt                <  1:417820/148‑1 (MQ=255)
                                                                                                                                     cAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATCAGGCCATTTTCATCGGGTAAGTTGATATCCAGCAGAATTAAATCTACCGACTGATTCTGCATAATTTCCCGTAg              >  4:143202/1‑148 (MQ=255)
                                                                                                                                     cAATCCGATAGCTGCGTCCGGTAAACAGAATAATCCCCACCGTTGAGCGTTCTCGCATGGCGCGGGTTAACATCTGGCCATTTTCATCGGGTAAGTTGATATCAATCAGAATTAAATCTa                                          <  8:142589/120‑1 (MQ=255)
                                                                                                                                      aaTCCGATCGCTGCGTCCTGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGAAGGGCGCGGGTTAACATCAGGCCATTTTCATCGGGTAAGTTGATATCCAGCAGAATTAAATCTAc                                         >  7:142589/1‑120 (MQ=255)
                                                                                                                                        tCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATCAGGCCATTTTCATCGGGTAAGTTGATATCCAGCAGAATTAAATCTACCGACTGATTCTGCAt                          <  4:203894/133‑1 (MQ=255)
                                                                                                                                        tCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATCAGGCCATTTTCATCGGGTAAGTTGATATCCAGCAGAATTAAATCTACCGACTGATTCTGCAt                          >  3:203894/1‑133 (MQ=255)
                                                                                                                                             tCGCTACGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATCAGGCCATTTTCATCGGGTAAGTTGATATCCAGCAGAATTAAATCTACCGACTGATTCTGCATAATTTcc                   >  6:91859/1‑135 (MQ=255)
                                                                                                                                             tCGCTACGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATCAGGCCATTTTCATCGGGTAAGTTGATATCCAGCAGAATTAAATCTACCGACTGATTCTGCATAATTTcc                   <  1:437557/135‑1 (MQ=255)
                                                                                                                                             tCGCTACGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATCAGGCCATTTTCATCGGGTAAGTTGATATCCAGCAGAATTAAATCTACCGACTGATTCTGCATAATTTcc                   >  2:437557/1‑135 (MQ=255)
                                                                                                                                             tCGCTACGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATCAGGCCATTTTCATCGGGTAAGTTGATATCCAGCAGAATTAAATCTACCGACTGATTCTGCATAATTTcc                   <  5:91859/135‑1 (MQ=255)
                                                                                                                                                cTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATCAGGCCATTTTCATCGGGTAAGTTGATATCCAGCAGAATTAAATCTACCGACTGATTCTGCATAATTTCCCGTAGCCCAGCACCg    >  4:66794/1‑147 (MQ=255)
                                                                                                                                                 tGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATCAGGCCATTTTCATCGGGTAAGTTGATATCCAGCAGAATTAAATCTACCGACTGATTCTGCATAATTTCCCGTAGCCCAGCACCGCt  >  2:426055/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TGAGTGTGCGGTTGAGCTTGTCGCGCGAGGTCGATTCGCCAGAGCAGATTTTTCACCCGTACTACCAGTTCGCGCAGTTCCAGCGGTTTGGTGACGTAATCGTCTGCGCCCATTTCCAGCCCAACAATACGGTCAATCCGATCGCTGCGTCCGGTAACCAGAATAATCCCCACCGTTGAGCGTTCTCGCAGGGCGCGGGTTAACATCAGGCCATTTTCATCGGGTAAGTTGATATCCAGCAGAATTAAATCTACCGACTGATTCTGCATAATTTCCCGTAGCCCAGCACCGCT  >  NC_000913/1057563‑1057855

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