Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,997,981 +G 96.6% intergenic (‑46/‑119) RCAP_rcc01844 ← / → RCAP_rcc01845 sensor histidine kinase/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,997,9811.→G96.6% 13.6 0.0 ‑0.5207intergenic (‑46/‑119)RCAP_rcc01844/RCAP_rcc01845sensor histidine kinase/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (100/100):  ref base (4/3):  total (104/103)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.12e-01

Only 100 of 220 total aligned reads displayed.

GCAGAATCGGAAACGGTAACGTCTCGGCCAAGGAGATTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTC.GGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCTTTTGCTGAATTCGCGGATGTTGTTGTCATGAGCGATGGGTTT  >  NC_014034/1997881‑1998076
                                                                                                     |                                                                                               
gCAGAATCGGAAACGGTAACGTCTCGGCCAAGGAGATTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATTCTTGACCCGAGCAAGCTTc                                                                                                  <  2:2442848/101‑1 (MQ=255)
gCAGAATCGGAAACGGTAACGTCTCGGCCAAGGAGATTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTc                                                                                                  <  2:4284825/101‑1 (MQ=255)
gCAGAATCGGAAACGGTAACGTCTCGGCCAAGGAGATTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTc                                                                                                  >  2:3996605/1‑101 (MQ=255)
   gAATCGGAAACGGTAAGGTCTCGGCCAAGGAGATTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTC.ggg                                                                                              >  1:3896279/1‑101 (MQ=255)
     aTCGGAAACGGTAACGTCTCGGCCAAGGAGATTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTC.gggga                                                                                            >  2:4435138/1‑99 (MQ=255)
       cGGAAACGGTAACGTCTCGGCCAAGGAGATTGTTGCCGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTC.GGGgaaa                                                                                          <  2:864618/101‑2 (MQ=255)
        ggAAACGGTAACGTCTCGGCCAAGGAGATTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTC.GGGgaaac                                                                                         >  1:3976941/1‑99 (MQ=255)
        ggAAACGGTAACGTCTCGGCCAAGGAGATTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTC.GGGgaaac                                                                                         >  2:3523169/1‑99 (MQ=255)
          aaaCGGTAACGTCTCGGCCAAGGAGATTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGa                                                                                        <  1:364712/101‑1 (MQ=255)
           aaCGGTAACGTCTCGGCCAAGGAGATTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGAt                                                                                       >  2:1645687/1‑101 (MQ=255)
              ggTAACGTCTCGGCCAAGGAGATTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATAcc                                                                                    >  1:2341087/1‑101 (MQ=255)
              ggTAACGTCTCGGCCAAGGAGATTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATAcc                                                                                    <  1:2139505/101‑1 (MQ=255)
                tAACGTCTCGGCCAAGGAGATTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCa                                                                                  <  1:1115951/101‑1 (MQ=255)
                tAACGTCTCGGCCAAGGAGATTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCa                                                                                  >  1:1302747/1‑101 (MQ=255)
                tAACGTCTCGGCCAAGGAGATTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCa                                                                                  <  1:2407986/101‑1 (MQ=255)
                 aaCGTCTCGGCCAAGGAGATTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCaa                                                                                 <  1:3103644/101‑1 (MQ=255)
                 aaCGTCTCGGCCAAGGAGATTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCaa                                                                                 >  2:3669122/1‑101 (MQ=255)
                  aCGTCTCGGCCAAGGAGATTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAg                                                                                >  2:3586887/1‑101 (MQ=255)
                  aCGTCTCGGCCAAGGAGATTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAg                                                                                <  1:3523417/101‑1 (MQ=255)
                  aCGTCTCGGCCAAGGAGATTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAg                                                                                <  1:575591/101‑1 (MQ=255)
                   cGTCTCGGCCAAGGAGATTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGt                                                                               >  2:3332407/1‑101 (MQ=255)
                    gTCTCGGCCAAGGAGATTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGtt                                                                              <  1:4068629/101‑1 (MQ=255)
                    gTCTCGGCCAAGGAGATTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGtt                                                                              <  2:635950/101‑1 (MQ=255)
                      ctcGGCCAAGGAGATTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCa                                                                            <  1:4189859/101‑1 (MQ=255)
                       tcGGCCAAGGAGATTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCaa                                                                           >  1:1616346/1‑101 (MQ=255)
                        cGGCCAAGGAGATTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAg                                                                          <  1:644163/101‑1 (MQ=255)
                        cGGCCAAGGAGATTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAg                                                                          <  2:1302650/101‑1 (MQ=255)
                         ggCCAAGGAGATTGTTGACGAATCCGCCATTTTTCGGGCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAgg                                                                         <  2:2981115/101‑1 (MQ=255)
                          gCCAAGGAGATTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGc                                                                        >  1:4370392/1‑101 (MQ=255)
                              aGGAGATTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCa                                                                    >  2:1335594/1‑101 (MQ=255)
                              aGGAGATTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCa                                                                    <  2:3434749/101‑1 (MQ=255)
                              aGGAGATTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCa                                                                    <  1:1053241/101‑1 (MQ=255)
                              aGGAGATTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCa                                                                    <  1:396947/101‑1 (MQ=255)
                               ggAGATTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCAt                                                                   <  1:997807/101‑1 (MQ=255)
                               ggAGATTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCAt                                                                   <  2:1662270/101‑1 (MQ=255)
                               ggAGATTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCAt                                                                   >  2:2920221/1‑101 (MQ=255)
                                 agaTTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGc                                                                 >  1:87262/1‑101 (MQ=255)
                                 agaTTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGc                                                                 <  1:2223854/101‑1 (MQ=255)
                                 agaTTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGc                                                                 <  1:210155/101‑1 (MQ=255)
                                   aTTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCaa                                                               <  1:1007079/101‑1 (MQ=255)
                                     tgttgACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCaaaa                                                             <  2:3255582/101‑1 (MQ=255)
                                     tgttgACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCaaaa                                                             <  1:2920442/101‑1 (MQ=255)
                                     tgttgACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCaaaa                                                             >  1:814151/1‑101 (MQ=255)
                                     tgttgACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCaaaa                                                             <  1:3954919/101‑1 (MQ=255)
                                        tgACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGcc                                                          <  2:771562/101‑1 (MQ=255)
                                        tgACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGcc                                                          <  1:3183665/101‑1 (MQ=255)
                                            gAATCCGCCATTTTTCGGTCCCCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAggg                                                      <  2:1308090/101‑1 (MQ=255)
                                            gAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAggg                                                      >  1:1374879/1‑101 (MQ=255)
                                            gAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAggg                                                      >  2:3157699/1‑101 (MQ=255)
                                             aaTCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAgggg                                                     >  1:1764505/1‑101 (MQ=255)
                                              aTCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAggggg                                                    >  2:4323176/1‑101 (MQ=255)
                                               tCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAgggggg                                                   >  2:175567/1‑101 (MQ=255)
                                                ccGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGc                                                  <  2:1025003/101‑1 (MQ=255)
                                                 cGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCt                                                 <  2:235191/101‑1 (MQ=255)
                                                  gCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTa                                                <  1:2033796/101‑1 (MQ=255)
                                                     aTTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGc                                             <  1:3405828/101‑1 (MQ=255)
                                                     aTTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGc                                             >  2:1097600/1‑101 (MQ=255)
                                                      tttttCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCt                                            >  1:4212544/1‑101 (MQ=255)
                                                        tttCGGTCACCCTGATATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCttt                                          >  2:3758794/1‑101 (MQ=255)
                                                        tttCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCttt                                          <  1:1851669/101‑1 (MQ=255)
                                                        tttCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCttt                                          >  2:2523907/1‑101 (MQ=255)
                                                         ttCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCtttt                                         >  1:3244034/1‑101 (MQ=255)
                                                         ttCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCtgtt                                         <  1:2559916/101‑1 (MQ=255)
                                                           cGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCTTTTGc                                       <  2:4345191/101‑1 (MQ=255)
                                                           cGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCTTTTGc                                       <  1:91273/101‑1 (MQ=255)
                                                            ggTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCTTTTGCt                                      <  2:1853032/101‑1 (MQ=255)
                                                            ggTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCTTTTGCt                                      <  2:2340937/101‑1 (MQ=255)
                                                             gTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCTTTTGCTg                                     >  2:4123239/1‑101 (MQ=255)
                                                              tCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCTTTTGCTGa                                    <  1:1444082/101‑1 (MQ=255)
                                                              tCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCTTTTGCTGa                                    >  1:542030/1‑101 (MQ=255)
                                                              tCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCTTTTGCTGa                                    >  2:3794382/1‑101 (MQ=255)
                                                               cACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCTTTTGCTGaa                                   <  2:152787/101‑1 (MQ=255)
                                                                aCCCTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCTTTTGCTGAAt                                  <  1:2422889/101‑1 (MQ=255)
                                                                 cccTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCTTTTGCTGAAtt                                 <  1:2954939/101‑1 (MQ=255)
                                                                 cccTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCTTTTGCTGAAtt                                 <  2:802628/101‑1 (MQ=255)
                                                                  ccTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCTTTTGCTGAATTc                                <  1:3378001/101‑1 (MQ=255)
                                                                   cTGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCTTTTGCTGAATTcg                               <  1:2621334/101‑1 (MQ=255)
                                                                    tGACATATCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCTTTTGCTGAATTcgc                              >  2:2792882/1‑101 (MQ=255)
                                                                         tatCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCTTTTGCTGAATTCGCGGAtg                         <  2:2582988/101‑1 (MQ=255)
                                                                           tCCAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCTTTTGCTGAATTCGCGGAtgtt                       >  2:2080980/1‑101 (MQ=255)
                                                                            ccAATGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCTTTTGCTGAATTCGCGGAtgctc                      >  1:4139051/1‑98 (MQ=255)
                                                                              aaTGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCTTTTGCTGAATTCGCGGAtgttgtt                    >  2:707742/1‑101 (MQ=255)
                                                                               aTGCTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCTTTTGCTGAATTCGCGGAtgttgttg                   <  2:1603463/101‑1 (MQ=255)
                                                                                  cTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCTTTTGCTGAATTCGCGGATGTTGTTGTCa                >  1:684275/1‑101 (MQ=255)
                                                                                  cTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCTTTTGCTGAATTCGCGGATGTTGTTGTCa                >  2:4209899/1‑101 (MQ=255)
                                                                                  cTTGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCTTTTGCTGAATTCGCGGATGTTGTTGTCa                <  1:1645793/101‑1 (MQ=255)
                                                                                    tGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAGGCTTTTGCTGAATTCGCGGATGTTGTTGTCATg              >  1:3787307/1‑101 (MQ=255)
                                                                                    tGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCTTTTGCTGAATTCGCGGATGTTGTTGTCATg              >  2:4247602/1‑101 (MQ=255)
                                                                                    tGACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCTTTTGCTGAATTCGCGGATGTTGTTGTCATg              >  1:2527351/1‑101 (MQ=255)
                                                                                     gccccGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCTTTTGCTGAATTCGCGGATGTTGTTGTCATGa             <  1:513870/99‑1 (MQ=255)
                                                                                     gACCCGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCTTTTGCTGAATTCGCGGATGTTGTTGTCATGa             >  1:171497/1‑101 (MQ=255)
                                                                                       cccGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCTTTTGCTGAATTCGCGGATGTTGTTGTCATGAGc           >  2:1205040/1‑101 (MQ=255)
                                                                                         cGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCTTTTGCTGAATTCGCGGATGTTGTTGTCATGAGCGa         <  1:2744248/101‑1 (MQ=255)
                                                                                         cGAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCTTTTGCTGAATTCGCGGATGTTGTTGTCATGAGCGa         >  1:2801399/1‑101 (MQ=255)
                                                                                          gAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCTTTTGCTGAATTCGCGGATGTTGTTGTCATGAGCGAt        >  2:3464969/1‑101 (MQ=255)
                                                                                          gAGCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGATTTTGCTGAATTCGCGGATGTTGTTGTCATGAGCGAt        >  1:3813259/1‑101 (MQ=255)
                                                                                            gCAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGCCTAAGCTTTTGCTGAATCGCCGGATGTTGTTGTAATGGGCGATgg      >  1:1567766/1‑101 (MQ=255)
                                                                                             cAAGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCTTTTGCTGAATTCGCGGATGTTGTTGTCATGAGCGATggg     >  1:1908058/1‑101 (MQ=255)
                                                                                              aaGCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCTTTTGCTGAATTCGCGGATGTTGTTGTCATGAGCGATGGGt    >  2:2700144/1‑101 (MQ=255)
                                                                                                gCTTCGGGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCTTTTGCTGAATTCGCGGATGTTGTTGTCATGAGCGATGGGttt  <  2:3294936/101‑1 (MQ=255)
                                                                                                     |                                                                                               
GCAGAATCGGAAACGGTAACGTCTCGGCCAAGGAGATTGTTGACGAATCCGCCATTTTTCGGTCACCCTGACATATCCAATGCTTGACCCGAGCAAGCTTC.GGGAAACGATACCCAAGTTCAAGGCCGCATGCAAAAGCCAGGGGGGCTAAGCTTTTGCTGAATTCGCGGATGTTGTTGTCATGAGCGATGGGTTT  >  NC_014034/1997881‑1998076

Base quality scores:  ATCG  < 3 ≤  ATCG  < 24 ≤  ATCG  < 31 ≤  ATCG  < 35 ≤  ATCG  < 41 ≤  ATCG