Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,960,926 G→C 5.7% G501G (GGC→GGG) ‡ 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,9260GC 5.7% 422.6 / 13.1 136G501G (GGC→GGG) ‡recCexonuclease V (RecBCD complex), gamma chain
Reads supporting (aligned to +/- strand):  ref base G (55/73);  new base C (5/3);  total (60/76)
Fisher's exact test for biased strand distribution p-value = 3.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.46e-01

Only 100 of 137 total aligned reads displayed.

CATAAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACGCTTC  >  NC_000913/2960778‑2961075
                                                                                                                                                    |                                                                                                                                                     
catAATGTAGAAGCGTTTGCCACTCGCGCTGCGCGCTCTGCATGGCGTAGGCCATCAAGATACGCGTCTGGCCAAATCGCCAGGTGTGTTGTCCGGTGGAGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGc                                                                                                                                                      <  2:1738982/150‑1 (MQ=255)
  tAAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCACCAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAAGGCCg                                                                                                                                                    >  2:153271/1‑150 (MQ=255)
  tAAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACACCATACGCGTCAGGCCAAATCGCCAGGTGGGTTGTGCGGAGGCGGGGAGTTCCCGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGcgg                                                                                                                                                    >  1:136809/1‑148 (MQ=255)
       tAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTc                                                                                                                                               >  2:530382/1‑150 (MQ=255)
        aGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACACCATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCg                                                                                                                                              <  2:530930/150‑1 (MQ=255)
        aGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATCCGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCg                                                                                                                                              <  2:16052/150‑1 (MQ=255)
          aaCCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGtt                                                                                                                                            <  1:1740789/150‑1 (MQ=255)
          aaCCGATTGCCACTCGCCCTGAGCGCTCTCCATCTCGTAGCCCAACAACATACGCGTCAGGCCAAATCGACAGGTGTGTTGTCCGGTGTCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGtt                                                                                                                                            <  1:221492/150‑1 (MQ=255)
             cGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTTAc                                                                                                                                         >  2:1462267/1‑150 (MQ=255)
             cGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGAc                                                                                                                                         <  1:650321/150‑1 (MQ=255)
             cGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGAc                                                                                                                                         <  1:242206/150‑1 (MQ=255)
             cGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATCCGGCATTCGTTggg                                                                                                                                         <  1:1726293/150‑3 (MQ=255)
              gATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGAcc                                                                                                                                        <  1:924979/150‑1 (MQ=255)
              gATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATCCGGCATTCGTTgggg                                                                                                                                        <  1:1439513/150‑4 (MQ=255)
              gATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATCCGGCATTCGTTgggg                                                                                                                                        <  2:472351/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)
                        cGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCCCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTa                                                                                                                              >  1:1537192/1‑150 (MQ=255)
                         gCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGAGTGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCTACGAATGCCGGATTCGTTGACCCACTGGCGTaa                                                                                                                             <  2:859893/150‑1 (MQ=255)
                          cccTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTaaa                                                                                                                            <  1:1276951/150‑1 (MQ=255)
                           ccTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATCCCGGATTCGTTGACCCACTGGCGTaact                                                                                                                           >  1:1581157/1‑148 (MQ=255)
                            ctccgcgcTCTCCATCGCGGAGCGCAACAACATACGCGTCAGGCCATCTCGCCAGGGGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGGCGTCATCGATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATa                                                                                                                          <  1:1935456/147‑1 (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)
                                       ccATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGCACGTTGTCGTCATCTATGCCCCAACGAATGCCGGAATCGTTGCCCCACTGGCGTAACTACCGCCGcccgt                                                                                                               >  2:1207422/1‑148 (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)
                                                aGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATACCGCAGCCCTTCTTCGGTGa                                                                                                      >  2:1892276/1‑150 (MQ=255)
                                                aGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTAGAACGCTGTCGTCGACTTAACGCTAACAAACATAGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGa                                                                                                      >  2:1380153/1‑150 (MQ=255)
                                                  cccTACAACATACGCGTCAGCCCAAATCGCCAGGCGTTTTGGCCGGTGGAGGGGAGTTCCAGCTCGCGAACGTGGGCGTCATCTATGCCCCACCGAATGCCGTATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATg                                                                                                    <  2:1807717/150‑1 (MQ=255)
                                                  cccAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATg                                                                                                    <  1:1750322/150‑1 (MQ=255)
                                                  cccAACAACATACGCGTCAGGCCAAATCCCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATg                                                                                                    <  2:1537192/150‑1 (MQ=255)
                                                   ccAACACCATCCGCGTCAGGCCAACTCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGGCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGt                                                                                                   <  1:274707/150‑1 (MQ=255)
                                                    caacaaCATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCTGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTc                                                                                                  <  1:385423/150‑1 (MQ=255)
                                                    caacaaCATACGCGTCAGGCCAAAACGCCAGGTGTGTTGTCCGGTGGCGTAGAGTTCCAGATATCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCATTCTTCGGTGATGTc                                                                                                  >  2:179021/1‑150 (MQ=255)
                                                       caacaTACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCAATCTCGCGAAAGTTGTCGTCATCTATGCCCCAACGAATGCCCGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCaaa                                                                                               >  2:163129/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)
                                                                cgTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGcgccgc                                                                                      >  2:1702724/1‑150 (MQ=255)
                                                                 gTCAGGCCATATCGCCAGGTGTGTTGGCCGGTGGCGGGGAGTTCCAGCTCGCGAACTTTGTCGTCATCTATGCCCCTACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCgccgcc                                                                                     <  1:561227/150‑1 (MQ=255)
                                                                 gTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCgccgcc                                                                                     <  2:632775/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)
                                                                       ccAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCAcc                                                                               <  2:1841339/150‑1 (MQ=255)
                                                                            tCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGcac                                                                          <  1:230612/150‑1 (MQ=255)
                                                                             cGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGcaca                                                                         >  1:1617654/1‑150 (MQ=255)
                                                                             cGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGcaca                                                                         >  1:1212600/1‑150 (MQ=255)
                                                                                cAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATATCGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATcc                                                                      >  2:599241/1‑150 (MQ=255)
                                                                                cAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATcc                                                                      <  1:530382/150‑1 (MQ=255)
                                                                                  ggTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAg                                                                    <  2:835648/150‑1 (MQ=255)
                                                                                       gttgtCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACg                                                               <  1:414705/150‑1 (MQ=255)
                                                                                         tgtCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGcc                                                             <  1:72642/150‑1 (MQ=255)
                                                                                          gtCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCa                                                            >  2:741067/1‑150 (MQ=255)
                                                                                          gtCCGGTGGCGGGGAGTTCCAGCTAGCGAACGTTGTCGTCATCTATCCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATACAGCAACCCCa                                                            >  2:836704/1‑150 (MQ=255)
                                                                                                 ggCGGGGAGTTCCAGCTCGCGAACGTTGGCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATc                                                     <  2:1339593/150‑1 (MQ=255)
                                                                                                      ggAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTg                                                <  2:1322294/150‑1 (MQ=255)
                                                                                                       gAGTTCCAGCTCGCGAACGTTGGCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGa                                               <  1:835943/150‑1 (MQ=255)
                                                                                                        aGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGac                                              <  1:1133355/150‑1 (MQ=255)
                                                                                                        aGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGac                                              <  1:1462267/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)
                                                                                                                        cGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGGCCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAg                              <  2:1355614/150‑1 (MQ=255)
                                                                                                                          ttGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGc                            <  1:1061142/150‑1 (MQ=255)
                                                                                                                              cgtcATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCGCCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTCTCAGGCAggg                        >  1:693685/1‑148 (MQ=255)
                                                                                                                                tcATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGACAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACCTCCTCTGACACACAACGACTGGCAGGCAGTGAt                      >  2:79986/1‑150 (MQ=255)
                                                                                                                                 cATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATa                     >  2:37462/1‑150 (MQ=255)
                                                                                                                                 cATCAATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACAGCGCCGCCAGCACCGGCACATCCAGCAACGACAGCACATCCTCTGACACAAAACGACTGTCACGCAGTGaaa                     >  1:69780/1‑148 (MQ=255)
                                                                                                                                 aaTCTATGCCCCACCGAATGCCGGATGCGTGGACCCACTGGCGTAAATAAGGCACCCCTTCTTCGTTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATa                     <  2:1769449/149‑1 (MQ=255)
                                                                                                                                   tCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAAc                   <  2:1858900/150‑1 (MQ=255)
                                                                                                                                    cTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCTCATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACa                  >  2:295337/1‑150 (MQ=255)
                                                                                                                                     tATGCCCCAACGAATGCCGTATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCTCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAg                 <  2:159814/150‑1 (MQ=255)
                                                                                                                                      aTGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAgg                <  2:656630/150‑1 (MQ=255)
                                                                                                                                       tGCCCCAACGAATGCCGGATTCGTTGACCCACTAGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGc               <  2:1735841/150‑1 (MQ=255)
                                                                                                                                        gCCCCCACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCt              <  1:1803775/150‑1 (MQ=255)
                                                                                                                                         ggtcAACGAATCCGGCATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTg             >  2:1634378/4‑150 (MQ=255)
                                                                                                                                         ccccAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTg             <  2:154771/150‑1 (MQ=255)
                                                                                                                                         ccccAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTg             <  2:614731/150‑1 (MQ=255)
                                                                                                                                          gtcAACGAATCCGGCATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGa            >  1:802295/3‑150 (MQ=255)
                                                                                                                                          gtcAACGAATCCGGCATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCAAATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGa            >  1:23439/3‑150 (MQ=255)
                                                                                                                                                gAATGCCGGATTCGTCGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGTCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACg      >  2:16792/1‑150 (MQ=255)
                                                                                                                                                  aTGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGGGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACGCt    <  1:1728580/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CATAAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACGCTTC  >  NC_000913/2960778‑2961075

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: