Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,960,928 C→G 5.2% G501R (GGC→CGC) ‡ recC ← exonuclease V (RecBCD complex), gamma chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,960,9280CG 5.2% 427.5 / 12.1 134G501R (GGC→CGC) ‡recCexonuclease V (RecBCD complex), gamma chain
Reads supporting (aligned to +/- strand):  ref base C (54/73);  new base G (4/3);  total (58/76)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.77e-01

Only 100 of 137 total aligned reads displayed.

TAAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACGCTTCC  >  NC_000913/2960780‑2961076
                                                                                                                                                    |                                                                                                                                                    
tAAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCACCAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAAGGCCg                                                                                                                                                     >  2:153271/1‑150 (MQ=255)
tAAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCg                                                                                                                                                     <  1:155170/150‑1 (MQ=255)
     tAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTc                                                                                                                                                >  2:530382/1‑150 (MQ=255)
      aGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACACCATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCg                                                                                                                                               <  2:530930/150‑1 (MQ=255)
       gAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCTTCTATGCCCCAACGAATGCCGGATTCGt                                                                                                                                              >  2:1634394/1‑150 (MQ=255)
        aaCCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGtt                                                                                                                                             <  2:352933/150‑1 (MQ=255)
         aCCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTAGTCCGCTGGCGGGGAGTTCCAAATAGATAACGATGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTg                                                                                                                                            >  2:643348/1‑150 (MQ=255)
           cGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGAc                                                                                                                                          <  1:650321/150‑1 (MQ=255)
            gATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGAcc                                                                                                                                         <  1:924979/150‑1 (MQ=255)
            gATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATCCGGCATTCGTTgggg                                                                                                                                         <  1:1439513/150‑4 (MQ=255)
              ttGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATAAGCGTCAGGACAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTCTGCCCCAACGCATGCCGGATTCGTTGAcccc                                                                                                                                       >  1:154771/1‑149 (MQ=255)
               tGCCACTCACCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGGTGACCCAc                                                                                                                                      >  1:155043/1‑150 (MQ=255)
                gCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACt                                                                                                                                     <  2:1225969/150‑1 (MQ=255)
                 ccACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTAGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCCCTg                                                                                                                                    >  2:432630/1‑150 (MQ=255)
                  cACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACATCATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGAGTCGTTGACCCACTgg                                                                                                                                   <  1:645697/150‑1 (MQ=255)
                  cACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGGTGACCCACTgg                                                                                                                                   >  1:1529521/1‑150 (MQ=255)
                    cTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACCACATACCCGTCAGGCCAAAACGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCg                                                                                                                                 <  1:1883981/150‑1 (MQ=255)
                    cTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCg                                                                                                                                 <  1:884790/150‑1 (MQ=255)
                    cTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCg                                                                                                                                 <  1:840936/150‑1 (MQ=255)
                       gCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGAGTGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCTACGAATGCCGGATTCGTTGACCCACTGGCGTaa                                                                                                                              <  2:859893/150‑1 (MQ=255)
                        cccTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTaaa                                                                                                                             <  1:1276951/150‑1 (MQ=255)
                          ctccgcgcTCTCCATCGCGGAGCGCAACAACATACGCGTCAGGCCATCTCGCCAGGGGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGGCGTCATCGATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATa                                                                                                                           <  1:1935456/147‑1 (MQ=255)
                          cTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATa                                                                                                                           <  2:1983297/150‑1 (MQ=255)
                                gcTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAg                                                                                                                     >  2:1745205/1‑150 (MQ=255)
                                 ctctCCCTCGCGTAGCCCAACAACCTACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGCATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGc                                                                                                                    >  1:1584539/1‑150 (MQ=255)
                                 ctctCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCGCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGc                                                                                                                    <  1:413369/150‑1 (MQ=255)
                                 ctctCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGc                                                                                                                    >  2:1728580/1‑150 (MQ=255)
                                   ctcCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGccc                                                                                                                  >  2:373078/1‑150 (MQ=255)
                                         cgcgGAGCCCCACAACATACGCGGCACGCCAAATCGCCATGGGTGTTGGCCGGGGGCCGGGTGTTCCATCTCGCGTACGTTGTCGTCATCAATGCCCCAACGAATGCCGGATTCGGTGACCCACTGGCGTAAATAACGCAGCCCttcttc                                                                                                            <  2:2015692/150‑1 (MQ=255)
                                          gcgTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTAGAACTCTGTCGTCGACTTAAAAATAAAAAAAAATATCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGGAACTACCGCAGCCCTtcgtcg                                                                                                           >  2:657505/1‑150 (MQ=255)
                                            gTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGt                                                                                                         <  2:913974/150‑1 (MQ=255)
                                              aGCCCAACAACATCCGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGa                                                                                                       <  2:1167818/150‑1 (MQ=255)
                                                cccTACAACATACGCGTCAGCCCAAATCGCCAGGCGTTTTGGCCGGTGGAGGGGAGTTCCAGCTCGCGAACGTGGGCGTCATCTATGCCCCACCGAATGCCGTATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATg                                                                                                     <  2:1807717/150‑1 (MQ=255)
                                                cccAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATg                                                                                                     <  1:1750322/150‑1 (MQ=255)
                                                cccAACAACATACGCGTCAGGCAAAATCGCCAGGTGGGGTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCTTAATCTATGCCCCAACTAATGCCGGATTCGTTTACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATg                                                                                                     <  2:1467097/150‑1 (MQ=255)
                                                 ccAACACCATCCGCGTCAGGCCAACTCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGGCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGt                                                                                                    <  1:274707/150‑1 (MQ=255)
                                                  caacaaCATACGCGTCAGGGCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTc                                                                                                   <  2:1549568/150‑1 (MQ=255)
                                                  caacaaCATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCTGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTc                                                                                                   <  1:385423/150‑1 (MQ=255)
                                                  caacaaCATACGCGTCAGGCCAAAACGCCAGGTGTGTTGTCCGGTGGCGTAGAGTTCCAGATATCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCATTCTTCGGTGATGTc                                                                                                   >  2:179021/1‑150 (MQ=255)
                                                         aTACGCGTCAGGCCAACTCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACcgc                                                                                            <  2:1746962/150‑1 (MQ=255)
                                                         aTACGCGTCAGGCCAAATCGCCAGGTGTGTTGACCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACcgc                                                                                            >  1:1098793/1‑150 (MQ=255)
                                                          tACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACcgcg                                                                                           >  1:612169/1‑150 (MQ=255)
                                                          tACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGCTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGGAAATAACGCCGCCCGTCGTCGGTGATGTCAAACcgcg                                                                                           >  2:285110/1‑150 (MQ=255)
                                                            cgcgTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGcgcc                                                                                         >  2:1431333/1‑150 (MQ=255)
                                                             gcgTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGcgccg                                                                                        >  1:1974816/1‑150 (MQ=255)
                                                              cgTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGcgccgc                                                                                       >  2:1702724/1‑150 (MQ=255)
                                                              cgTCAGGCCAAATCGCCAGGTGTGTTGACCGGTGGCGGGGAGTTCCAGCTAGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGcgccgc                                                                                       >  2:42697/1‑150 (MQ=255)
                                                               gTCAGGCCATATCGCCAGGTGTGTTGGCCGGTGGCGGGGAGTTCCAGCTCGCGAACTTTGTCGTCATCTATGCCCCTACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCgccgcc                                                                                      <  1:561227/150‑1 (MQ=255)
                                                               gTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCgccgcc                                                                                      <  2:1605796/150‑1 (MQ=255)
                                                                 cAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAGCCGCGCCGCCAg                                                                                    <  2:1114678/150‑1 (MQ=255)
                                                                 cAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAg                                                                                    <  1:1745205/150‑1 (MQ=255)
                                                                       aaaTCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCgg                                                                              >  2:1845251/1‑150 (MQ=255)
                                                                          tCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGcac                                                                           <  1:230612/150‑1 (MQ=255)
                                                                           cGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGcaca                                                                          >  1:1212600/1‑150 (MQ=255)
                                                                           cGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGcaca                                                                          >  1:1617654/1‑150 (MQ=255)
                                                                              cAGGTGTGTTGTGCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCTGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATcc                                                                       <  2:1333/150‑1 (MQ=255)
                                                                              cAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATcc                                                                       <  1:530382/150‑1 (MQ=255)
                                                                              cAGGTGTGTTGTCCGGTGGCGGGGAGGGCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATcc                                                                       <  2:475479/150‑1 (MQ=255)
                                                                                ggTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAg                                                                     <  2:835648/150‑1 (MQ=255)
                                                                                  tgtgtTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCCGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCa                                                                   >  2:768537/1‑150 (MQ=255)
                                                                                     gttgtCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACg                                                                <  1:414705/150‑1 (MQ=255)
                                                                                       tgtCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGCCCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGAAACGcc                                                              >  1:1736340/1‑150 (MQ=255)
                                                                                       tgtCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATACCGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACCTCCAGCAACGcc                                                              >  1:239050/1‑150 (MQ=255)
                                                                                       tgtCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGcc                                                              <  1:72642/150‑1 (MQ=255)
                                                                                           cGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGca                                                          <  1:153271/150‑1 (MQ=255)
                                                                                               ggCGGGGAGTTCCAGCTCGCGAACGTTGGCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATc                                                      <  2:1339593/150‑1 (MQ=255)
                                                                                                     gAGTTCCAGCTCGCGAACGTTGGCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGa                                                <  1:835943/150‑1 (MQ=255)
                                                                                                      aGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGac                                               <  1:1133355/150‑1 (MQ=255)
                                                                                                          ccAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACaa                                           >  2:1877435/1‑150 (MQ=255)
                                                                                                             gCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAAc                                        <  1:1347221/150‑1 (MQ=255)
                                                                                                              cTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACg                                       >  2:1675300/1‑150 (MQ=255)
                                                                                                                cgcgAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGAc                                     <  1:1056585/150‑1 (MQ=255)
                                                                                                                 gcgAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACt                                    <  2:451856/150‑1 (MQ=255)
                                                                                                                  cgAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACAGGCGTAAAAAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCAACGGCACATCCCGCACCGCCAGCACATCCTCTGACCCCAAACGACTg                                   >  2:1319682/1‑150 (MQ=255)
                                                                                                                     aCGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCTGCAACGCCAGCACATCCTGTGACACAAAACGACTGTCa                                <  2:178165/150‑1 (MQ=255)
                                                                                                                      cGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGGCCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAg                               <  2:1355614/150‑1 (MQ=255)
                                                                                                                        ttGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGc                             <  1:1061142/150‑1 (MQ=255)
                                                                                                                            cgtcATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCGCCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTCTCAGGCAggg                         >  1:693685/1‑148 (MQ=255)
                                                                                                                            cgtcATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTg                         >  1:1460279/1‑150 (MQ=255)
                                                                                                                              tcATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGACAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACCTCCTCTGACACACAACGACTGGCAGGCAGTGAt                       >  2:79986/1‑150 (MQ=255)
                                                                                                                               cATCAATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACAGCGCCGCCAGCACCGGCACATCCAGCAACGACAGCACATCCTCTGACACAAAACGACTGTCACGCAGTGaaa                      >  1:69780/1‑148 (MQ=255)
                                                                                                                                 tCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAAc                    <  2:1858900/150‑1 (MQ=255)
                                                                                                                                   tATGCCCCAACGAATGCCGTATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCTCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAg                  <  2:159814/150‑1 (MQ=255)
                                                                                                                                    aTGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAgg                 <  2:656630/150‑1 (MQ=255)
                                                                                                                                      gCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACAGTCAGGCAGTGATAACAGGCt               <  2:1212600/150‑1 (MQ=255)
                                                                                                                                       ggtcAACGAATCCGGCATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTg              >  2:1634378/4‑150 (MQ=255)
                                                                                                                                       ccccAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTg              <  2:606643/150‑1 (MQ=255)
                                                                                                                                       ccccAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTg              <  2:614731/150‑1 (MQ=255)
                                                                                                                                       ccccAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTg              <  2:154771/150‑1 (MQ=255)
                                                                                                                                        gtcAACGAATCCGGCATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGCCACAAAACGACTGTCGGGCAGTGAGAACAGGCTGa             >  1:1791911/3‑150 (MQ=255)
                                                                                                                                        gtcAACGAATCCGGCATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGa             >  1:802295/3‑150 (MQ=255)
                                                                                                                                        gtcAACGAATCCGGCATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCAAATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGa             >  1:23439/3‑150 (MQ=255)
                                                                                                                                        cccAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCCCATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGAt            >  2:1381978/1‑150 (MQ=255)
                                                                                                                                          cAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATa           <  1:1634394/150‑1 (MQ=255)
                                                                                                                                                aTGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACGCt     <  2:1619079/150‑1 (MQ=255)
                                                                                                                                                aTGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGGGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACGCt     <  1:1728580/150‑1 (MQ=255)
                                                                                                                                                  gCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACGCTTc   <  2:1720917/150‑1 (MQ=255)
                                                                                                                                                   ccGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACGCTTcc  <  2:1116844/150‑1 (MQ=255)
                                                                                                                                                   ccGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACGCTTcc  <  2:1529521/150‑1 (MQ=255)
                                                                                                                                                   ccGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACGCTTcc  <  1:295337/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TAAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACGCTTCC  >  NC_000913/2960780‑2961076

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