Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 236,569 C→T 72.5% A168V (GCC→GTC)  dnaQ → DNA polymerase III epsilon subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913236,5690CT72.5% 277.1 / 128.0 182A168V (GCC→GTC) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base C (19/31);  new base T (73/58);  total (92/90)
Fisher's exact test for biased strand distribution p-value = 4.55e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 185 total aligned reads displayed.

GGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGGAGAGACACAACAGCAACAAGGTGAAGCAAC  >  NC_000913/236470‑236668
                                                                                                   |                                                                                                   
ggAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGc                                                                                                     >  2:2474319/1‑100 (MQ=255)
  aaaaTGTTTCCCGTTAAGCGCAACAGCCTTGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGATCGTTGCACGGGGCATTACTCGATGcc                                                                                                    <  1:2191978/99‑1 (MQ=255)
   aaaTGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCa                                                                                                  <  2:1372345/100‑1 (MQ=255)
   aaaTGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCa                                                                                                  >  2:3863906/1‑100 (MQ=255)
       gTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATc                                                                                              <  1:1566371/100‑1 (MQ=255)
        tttCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGAt                                                                                               <  2:143071/98‑1 (MQ=255)
         gatctggTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCt                                                                                            <  1:365646/95‑1 (MQ=255)
          tCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCtt                                                                                           >  2:908061/1‑100 (MQ=255)
           cccGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTg                                                                                          >  2:1803918/1‑100 (MQ=255)
            ccGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGc                                                                                         >  2:48439/1‑100 (MQ=255)
             cGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCg                                                                                        <  1:3429746/100‑1 (MQ=255)
             cGGTAAGCGCAACAGCCTCCATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCg                                                                                        <  1:3320168/100‑1 (MQ=255)
              ggTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCtggatc                                                                                       >  2:365647/1‑95 (MQ=255)
              ggTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCgg                                                                                       >  2:1027011/1‑100 (MQ=255)
               gTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGa                                                                                      <  1:657631/100‑1 (MQ=255)
                tAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGaa                                                                                     <  2:3188377/100‑1 (MQ=255)
                 aaGCGCAACAGCCTCGATGCGTAAGGTGCTCGCTACGAAATAGAGAACAGTAAACGAACGCTGCACGGGGCATTACTAGATGTCCAGATCCt                                                                                            <  1:2052385/92‑1 (MQ=255)
                  aGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGt                                                                                   >  1:2593117/1‑100 (MQ=255)
                        aCAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCt                                                                             >  1:411071/1‑100 (MQ=255)
                          aGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTgg                                                                           >  2:555281/1‑100 (MQ=255)
                           gCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGc                                                                          <  2:2161395/100‑1 (MQ=255)
                           gCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGc                                                                          <  1:1663656/100‑1 (MQ=255)
                             cgatctGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGa                                                                        <  2:60875/95‑1 (MQ=255)
                             cTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGa                                                                        >  1:3784217/1‑100 (MQ=255)
                              ttgATGCGTTATGTGCTTTTTACGAATTAGTTTATAGTATTTTAACGCTGCACGGTGTATTACTTGATGTTTAGTTTCTTGTGGAAGTTTATCTggt                                                                          <  1:1155154/95‑2 (MQ=255)
                              tCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTCCTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGAt                                                                       >  1:2636545/1‑100 (MQ=255)
                               cGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATg                                                                      <  2:197634/100‑1 (MQ=255)
                                    cGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCgg                                                                 >  1:2425714/1‑100 (MQ=255)
                                     gTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCggt                                                                >  2:1273593/1‑100 (MQ=255)
                                       tATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCggtgg                                                              <  1:832227/100‑1 (MQ=255)
                                        aTGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCggtggt                                                             <  2:2593131/100‑1 (MQ=255)
                                        aTGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCggtggt                                                             <  2:1097311/100‑1 (MQ=255)
                                         tgtgCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCggtgg                                                              <  1:357132/98‑1 (MQ=255)
                                         tgtgCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTc                                                            >  2:1514716/1‑100 (MQ=255)
                                          gtgCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGGCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTc                                                            <  1:3641410/99‑1 (MQ=255)
                                          gtgCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCa                                                           <  1:4135624/100‑1 (MQ=255)
                                            gCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCaaa                                                         >  1:1911148/1‑100 (MQ=255)
                                            gCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCaaa                                                         >  2:3418907/1‑100 (MQ=255)
                                            gCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCaaa                                                         >  2:289583/1‑100 (MQ=255)
                                             cTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAAc                                                        >  2:2470181/1‑100 (MQ=255)
                                              tCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACg                                                       >  2:3905418/1‑100 (MQ=255)
                                              tCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACg                                                       <  2:368881/100‑1 (MQ=255)
                                               cGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCaa                                                          >  1:2853377/1‑96 (MQ=255)
                                               cGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGt                                                      >  1:1849230/1‑100 (MQ=255)
                                               cGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGt                                                      >  2:4272253/1‑100 (MQ=255)
                                               cGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGt                                                      <  2:3982200/100‑1 (MQ=255)
                                                  ttcGAAATAGATAACAATAAACGAGCGCTGCACGGGGCATAACTCGATCTCCAGATCCTTGCGGAAGTTTATCTGGCGATGCCGGGTAGTCAAACg                                                       <  1:3175226/94‑1 (MQ=255)
                                                  tACGAAATAGATAATAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGa                                                   <  1:1421667/100‑1 (MQ=255)
                                                  tACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGa                                                   >  2:4065764/1‑100 (MQ=255)
                                                  tACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGa                                                   <  2:3728629/100‑1 (MQ=255)
                                                    cGAAGTAGATAACAGTATACGAACGCTGCACGAGTCATTACTCGATGTCCAGATCCTGGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACg                                                       <  1:1006602/94‑1 (MQ=255)
                                                    cGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAGGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATg                                                 >  2:3101205/1‑100 (MQ=255)
                                                    cGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATg                                                 <  1:9452/100‑1 (MQ=255)
                                                    cGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATg                                                 <  2:1922449/100‑1 (MQ=255)
                                                    cGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATg                                                 >  1:4273962/1‑100 (MQ=255)
                                                     gAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATgg                                                >  2:511962/1‑100 (MQ=255)
                                                      aaaTAGATAAGAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGc                                               <  2:4194722/100‑1 (MQ=255)
                                                       aaTAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCt                                              >  1:3255026/1‑100 (MQ=255)
                                                       aaTAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCt                                              <  1:3517339/100‑1 (MQ=255)
                                                       aaTAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCCATGACCGGTGGTCAAACGTCGATGGCt                                              >  2:3024822/1‑100 (MQ=255)
                                                         ccgatctcAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCttt                                            <  2:3240082/93‑1 (MQ=255)
                                                            aTAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGc                                         <  2:1584539/100‑1 (MQ=255)
                                                             tAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCg                                        >  2:1610311/1‑100 (MQ=255)
                                                             tAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCg                                        >  2:3904795/1‑100 (MQ=255)
                                                             tAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCg                                        >  1:968293/1‑100 (MQ=255)
                                                              aaCAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGa                                       <  1:3338569/100‑1 (MQ=255)
                                                               aCAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGAt                                      <  1:1044267/100‑1 (MQ=255)
                                                               aCAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGAt                                      <  1:4195525/100‑1 (MQ=255)
                                                                cAGTAAACGAACGCTGCACGGGGCATTACTCGCTGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATg                                     >  2:4181842/1‑100 (MQ=255)
                                                                cAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCtttagatcgg                                     >  1:3240064/1‑93 (MQ=255)
                                                                cAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATg                                     >  2:3032169/1‑100 (MQ=255)
                                                                cAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATg                                     >  2:1735536/1‑100 (MQ=255)
                                                                cAGTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATg                                     >  2:3192138/1‑100 (MQ=255)
                                                                 aGTAAACGAACGCTGCACGGGACATTACTCGATGCACAGAGAATTGTGGACGGTTAACTGGCGATGAGCGGAGGTCAAACGTAGATGGCTTTTGCg                                        >  1:2164012/1‑96 (MQ=255)
                                                                  gtaagggaacgCTGCACGGGGCATTACTCGACGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGa                                   <  1:151335/94‑1 (MQ=255)
                                                                  gTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATgg                                    >  1:1727958/1‑99 (MQ=255)
                                                                  gTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGa                                   <  2:2103338/100‑1 (MQ=255)
                                                                  gTAAACGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGa                                   >  1:3731401/1‑100 (MQ=255)
                                                                  gTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGa                                   <  1:592964/100‑1 (MQ=255)
                                                                    aaaCGAACGCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAg                                 <  1:2802269/100‑1 (MQ=255)
                                                                     aacgaacgCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAgg                                <  2:234763/100‑1 (MQ=255)
                                                                      acgaacgCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGga                               >  2:2496631/1‑100 (MQ=255)
                                                                       cgaacgCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGgag                              <  2:3426091/100‑1 (MQ=255)
                                                                       cgaacgCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGgag                              >  1:257746/1‑100 (MQ=255)
                                                                       cgaacgCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGgag                              <  2:820730/100‑1 (MQ=255)
                                                                          acgCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGgagaga                           >  2:2270769/1‑100 (MQ=255)
                                                                          acgCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGgagaga                           <  2:4205656/100‑1 (MQ=255)
                                                                           cgCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGGAGAGac                          <  2:2980832/100‑1 (MQ=255)
                                                                            gCTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGGAGAGaca                         <  1:2332934/100‑1 (MQ=255)
                                                                             cTGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGGAGAGacac                        >  2:2816422/1‑100 (MQ=255)
                                                                              tGCACGGGGCGTTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGGAGAGAcaca                       <  1:2896346/100‑1 (MQ=255)
                                                                              tGCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGGAGAGAcaca                       >  2:1943881/1‑100 (MQ=255)
                                                                               gCACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGCCCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGGAGAGACacaa                      >  2:1181903/1‑100 (MQ=255)
                                                                                cACGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTGTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGGAGAGACacaac                     <  1:2462255/100‑1 (MQ=255)
                                                                                  cGGGGCATTACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGGAGAGACACAACAg                   <  1:221276/100‑1 (MQ=255)
                                                                                          tACTCGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGGAGAGACACAACAGCAACAAgg           <  1:1933643/100‑1 (MQ=255)
                                                                                              cGATTCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCGTTTGCGATGGAACGAGAGAGACAACAGCAACAAGGTGaa       <  1:733073/100‑1 (MQ=255)
                                                                                              cGATGTCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGGAGAGACACAACAGCAACAAGGTGaa       >  1:2202023/1‑100 (MQ=255)
                                                                                               gATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGGAGAGACACAACAGCAACAAGGTGAAg      <  2:2668304/100‑1 (MQ=255)
                                                                                                   cccAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGGAGAGACACAACAGCAACAAGGTGAAGcaac  >  1:4050299/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAGGAGAGACACAACAGCAACAAGGTGAAGCAAC  >  NC_000913/236470‑236668

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