Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 207,483 A→C 5.9% V786V (GTA→GTC)  dnaE → DNA polymerase III alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913207,4830AC 5.9% 227.5 / 4.2 68V786V (GTA→GTC) dnaEDNA polymerase III alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (30/34);  new base C (2/2);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

CGCTGAACTGGCGATGAAAATCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTG  >  NC_000913/207336‑207631
                                                                                                                                                   |                                                                                                                                                    
cGCTGAACTGGCGATGAAAATCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTaa                                                                                                                                                     <  4:351115/149‑1 (MQ=255)
       cTGGCGATGAAAATCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgc                                                                                                                                              <  4:369211/149‑1 (MQ=255)
        tGGCGATGAAAATCTTCGACCAGGTGGAGGAAGGCGCTGGTTACGGATTTACCAAATAGCACTCTGCGGCCTATGCTTTGGTGTAATATCAAGCTTTATGGGTGTAAGCGCACTATGCTGCGGAGTGTATGGCGGCGGTAAAGAacgcc                                                                                                                                             >  4:59368/1‑149 (MQ=255)
          gCGATGAAAATCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGa                                                                                                                                           <  1:422206/149‑1 (MQ=255)
             aTGAAAATCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGatat                                                                                                                                        <  2:83337/149‑1 (MQ=255)
                aaaaTCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGa                                                                                                                                     >  1:147334/1‑149 (MQ=255)
                  aaTCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGaca                                                                                                                                   >  3:206374/1‑149 (MQ=255)
                  aaTCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGaca                                                                                                                                   >  2:563586/1‑149 (MQ=255)
                      ttCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAacac                                                                                                                               >  1:154628/1‑149 (MQ=255)
                         gACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCga                                                                                                                            >  2:346760/1‑149 (MQ=255)
                            cTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGaa                                                                                                                         <  2:396590/149‑1 (MQ=255)
                            cTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGaa                                                                                                                         <  1:465929/149‑1 (MQ=255)
                             tggtggAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGTCCGCCGATATGGACAACACCGAGAAg                                                                                                                        <  2:142004/149‑1 (MQ=255)
                             tggtggAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAg                                                                                                                        >  4:43689/1‑149 (MQ=255)
                              ggtggAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAgg                                                                                                                       <  1:563586/149‑1 (MQ=255)
                                tggAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAggtg                                                                                                                     >  2:467679/1‑149 (MQ=255)
                                 ggAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAggtgg                                                                                                                    >  4:555670/1‑149 (MQ=255)
                                  gagaAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAggtggt                                                                                                                   >  2:194245/1‑149 (MQ=255)
                                   agaAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGgtggtg                                                                                                                  <  1:3060/149‑1 (MQ=255)
                                          cGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCtgg                                                                                                           <  2:154628/149‑1 (MQ=255)
                                               gTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGAt                                                                                                      <  3:40763/149‑1 (MQ=255)
                                                 tACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGa                                                                                                    <  4:294201/149‑1 (MQ=255)
                                                   cGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGt                                                                                                  <  1:445902/149‑1 (MQ=255)
                                                    ggATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCtgg                                                                                                           >  4:282896/1‑139 (MQ=255)
                                                    ggATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCtgg                                                                                                           <  3:282896/139‑1 (MQ=255)
                                                     gATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGc                                                                                                >  2:184364/1‑149 (MQ=255)
                                                      aTTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCt                                                                                               <  3:46189/149‑1 (MQ=255)
                                                        ttAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTgg                                                                                             <  3:436033/149‑1 (MQ=255)
                                                             aaaTCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGAt                                                                                        >  2:319695/1‑149 (MQ=255)
                                                             aaaTCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGAt                                                                                        >  3:387348/1‑149 (MQ=255)
                                                               aTCGCACTCTGCGGCCTATGCTTTGATGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATgg                                                                                      <  1:409297/149‑1 (MQ=255)
                                                                 cGCACTCTGCGGCCTATGCTGTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATgggg                                                                                    <  4:107702/149‑1 (MQ=255)
                                                                     ctctGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGa                                                                                >  4:233631/1‑149 (MQ=255)
                                                                      tctGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGaa                                                                               >  4:184457/1‑149 (MQ=255)
                                                                          cGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATc                                                                           <  1:569862/149‑1 (MQ=255)
                                                                            gCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCt                                                                         <  2:147334/149‑1 (MQ=255)
                                                                              cTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGc                                                                       <  3:233631/149‑1 (MQ=255)
                                                                                   cTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGccacca                                                                  <  1:362905/149‑1 (MQ=255)
                                                                                     ttGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGa                                                                <  1:228000/149‑1 (MQ=255)
                                                                                     ttGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGa                                                                <  4:331706/149‑1 (MQ=255)
                                                                                        gtgtCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGatat                                                             >  2:275729/1‑149 (MQ=255)
                                                                                           tCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGTCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGAAATaaa                                                          >  2:239438/1‑149 (MQ=255)
                                                                                               atCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTcc                                                      >  1:63744/1‑149 (MQ=255)
                                                                                                tCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCg                                                     >  1:318255/1‑149 (MQ=255)
                                                                                                   aaCGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGc                                                                       <  1:467315/128‑1 (MQ=255)
                                                                                                   aaCGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGc                                                                       >  2:467315/1‑128 (MQ=255)
                                                                                                      gTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATACAAACTCCGGTCttt                                               <  3:43689/149‑1 (MQ=255)
                                                                                                       ttATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTa                                              >  4:524395/1‑149 (MQ=255)
                                                                                                        tATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCGGATATAAACTCCGGTCTTTAc                                             <  3:555670/149‑1 (MQ=255)
                                                                                                        tATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTAc                                             >  1:57754/1‑149 (MQ=255)
                                                                                                         aTGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTAcc                                            >  2:21791/1‑149 (MQ=255)
                                                                                                             cTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCAttt                                        >  2:232751/1‑149 (MQ=255)
                                                                                                                    cgcACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTc                                 <  2:383752/149‑1 (MQ=255)
                                                                                                                      cACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCaa                               <  1:232751/149‑1 (MQ=255)
                                                                                                                       aCTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAac                              <  2:318255/149‑1 (MQ=255)
                                                                                                                           tCCTGCGGAGTTTATGGCGGCGATAATGACCGCTGATATGGACAACACCGAGAAGGTGGTGGGGCTGGTGGATGAGTGCTGGTGGATGGGGCTGAACATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAacgacg                          >  2:451693/1‑149 (MQ=255)
                                                                                                                                 ggAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCg                    >  1:174974/1‑149 (MQ=255)
                                                                                                                                     tttATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCt                                                                                                             >  1:208186/1‑56 (MQ=255)
                                                                                                                                     tttATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCt                                                                                                             <  2:208186/56‑1 (MQ=255)
                                                                                                                                     tttATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAAt                >  1:597139/1‑149 (MQ=255)
                                                                                                                                       tATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCg              <  3:184457/149‑1 (MQ=255)
                                                                                                                                        atatcggcggTCATTACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTACGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCgt             >  2:117863/5‑149 (MQ=255)
                                                                                                                                         tGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCgtg            <  4:387348/149‑1 (MQ=255)
                                                                                                                                           tcggcggTCATTACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTa          >  4:384596/2‑149 (MQ=255)
                                                                                                                                            cggcggTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCAtt                                         >  4:15460/1‑117 (MQ=255)
                                                                                                                                            cggcggTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCAtt                                         <  3:15460/117‑1 (MQ=255)
                                                                                                                                             ggcggTCATTACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATg        <  3:384596/149‑1 (MQ=255)
                                                                                                                                             ggcggTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATg        <  1:239438/149‑1 (MQ=255)
                                                                                                                                                  tAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTAtt   <  3:200953/149‑1 (MQ=255)
                                                                                                                                                   aaTGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTg  <  2:174974/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGCTGAACTGGCGATGAAAATCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTG  >  NC_000913/207336‑207631

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