Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 205,220 C→T 5.9% A32V (GCG→GTG)  dnaE → DNA polymerase III alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913205,2200CT 5.9% 212.2 / 2.4 67A32V (GCG→GTG) dnaEDNA polymerase III alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (29/34);  new base T (2/2);  total (31/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 = 1.00e+00

TGGGACTTGCGTCCTGATTCTTGTGTCGAGATTAAGTAAACCGGAATCTGAAGATGTCTGAACCACGTTTCGTACACCTGCGGGTGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCGGGACATGGCGCAGGGATTAAGCCTATCGTCGGGGCAGATTTTAACGTCCAGTGCGACCTGCTGGGTGATGAGT  >  NC_000913/205073‑205366
                                                                                                                                                   |                                                                                                                                                  
tGGGACTTGCGTCCTGATTCTTGTGTCGAGATTAAGTAAACCGGAATCTGAAGATGTCTGAACCACGATTCGTACACCTGCGCGTGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAgccggcg                                                                                                                                                   >  6:313792/1‑149 (MQ=255)
   gACTTGCGTCCTGATTCTTGTGTCGAGATTAAGTAAACCGGAATCTGAAGATGTCTGAACCACGTTTCGTACACCTGCGGGTGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGgcggcggcg                                                                                                                                                >  2:190174/1‑149 (MQ=255)
        gCGTCCTGATTCTTGTGTCGAGATTAAGTAAACCGGAATCTGCAGATGTCTGAACCACGTTTCGTACACCTGCGGGTGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTgg                                                                                                                                            <  7:379851/148‑1 (MQ=255)
        gCGTCCTGATTCTTGTGTCGAGATTAAGTAAACCGGAATCTGAAGATGTCTGAACCACGTTTCGTACACCTGCGGGTGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTgg                                                                                                                                            <  8:197359/148‑1 (MQ=255)
        gCGTCCTGATTCTTGTGTCGAGATTAAGTAAACCGGAATCTGAAGATGTCTGAACCACGTTTCGTACACCTGCGGGTGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTg                                                                                                                                             <  2:360889/147‑1 (MQ=255)
         cGTCCTGATTCTTGTGTCGAGATTAAGTAAACCGGAATCTGAAGATGTCTGAACCACGTTTCGTACACCTGCGGGTGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGt                                                                                                                                          <  1:265738/149‑1 (MQ=255)
          gTCCTGATTCTTGTGTCGAGATTAAGTAAACCGGAATCTGAAGATGTCTGAACCACGTTTCGTACACCTGCGGGTGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTa                                                                                                                                         <  5:297102/149‑1 (MQ=255)
             cTGATTCTTGTGTCGAGATTAAGTAAACCGGAATCTGAAGATGTCTGAACCACGTTTCGTACACCTGCGGGTGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATg                                                                                                                                       >  3:293010/1‑148 (MQ=255)
             cTGATTCTTGTGTCGAGATTAAGTAAACCGGAATCTGAAGATGTCTGAACCACGTTTCGTACACCTGCGGGTGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGc                                                                                                                                      <  8:204828/149‑1 (MQ=255)
             cTGATTCTTGTGTCGAGATTAAGTAAACCGGAATCTGAAGATGTCTGAACCACGTTTCGTACACCTGCGGGTGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGc                                                                                                                                      <  4:293010/149‑1 (MQ=255)
                      gtgtCGAGATTAAGTAAACCGGAATCTGAAGATGTCTGAACCACGTTTCGTACACCTGCGGGTGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTg                                                                                                                              <  4:219384/148‑1 (MQ=255)
                       tgtCGAGATTAAGTAAACCGGAATCTGAAGATGTCTGAACCACGTTTCGTACACCTGCGGGTGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGc                                                                                                                            >  4:395257/1‑149 (MQ=255)
                        gtCGAGATTAAGTAAACCGGAATCTGAAGATGTCTGACCCACGTTTCGTACACCTGCGGGTGCACAGCGCCTACTCGATGATCGATGGCCTGGCCACACCCGCACCGTTGGGAAAAAAGGCGGCGGCGGTGGGTATGCCAGCACTGGc                                                                                                                            >  5:329801/1‑148 (MQ=255)
                                   gTAAACCGGAATCTGAAGATGTCTGAACCACGTTTCATACACCTGCGGGTGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGAttt                                                                                                                >  8:305119/1‑149 (MQ=255)
                                    tAAACCGGAATCTGAAGATGTCTGAACCACGTTTCGTACACCTGCGGGTGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGAttt                                                                                                                >  2:42834/1‑148 (MQ=255)
                                      aaCCGGAATCTGAAGCTGGCTGCACCCCGTTTCGTACACCTGCGGGGGCCCAGCGCCTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTGGGGTATGCCAGCACTGGCGATCACCGACTTCa                                                                                                              <  1:36421/148‑1 (MQ=255)
                                       aCCGGAATCTGAAGATGTCTGAACCACGTTTCGTACACCTGCGGGTGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAATCCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCAcc                                                                                                            <  6:301081/149‑1 (MQ=255)
                                         cGGAATCTGAAGATGTCTGAACCACGTTTCGTACACCTGCGGGTGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCaa                                                                                                          >  1:271507/1‑149 (MQ=255)
                                                 gAAGATGTCTGAACCACGTTTCGTACACCTGCGGGTGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTtgtg                                                                                                  >  4:354851/1‑149 (MQ=255)
                                                  aaGATGTCTGAACCACGTTTCGTACACCTGCGGGTGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTtgtg                                                                                                  <  4:426213/148‑1 (MQ=255)
                                                    gATGTCTGAACCACGTTTCGTACACCTGCGGGTGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTc                                                                                               >  3:115401/1‑149 (MQ=255)
                                                     aTGTCTGAACCACGTTTCGTACACCTGCGGGTGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCt                                                                                              <  7:53819/149‑1 (MQ=255)
                                                          tGAACCACGTTTCGTACACCTGCGGGTGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTg                                                                                          >  8:263314/1‑148 (MQ=255)
                                                           gAACCACGTTTCGTACACCTGCGGGTGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTg                                                                                          >  8:95003/1‑147 (MQ=255)
                                                               cACGTTTCGTACACCTGCGGGTGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTc                                                                                    <  4:158190/149‑1 (MQ=255)
                                                                  gTTTCGTACACCTGCGGGTGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTAc                                                                                 <  2:271507/149‑1 (MQ=255)
                                                                  gTTTCATACACCTGCGGGTGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTAc                                                                                 <  7:305119/149‑1 (MQ=255)
                                                                   tttCGTACACCTGCGGGTGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACg                                                                                >  8:9493/1‑149 (MQ=255)
                                                                        tACACCTGCGGGTGCACAGCGTCTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCg                                                                           <  6:95543/149‑1 (MQ=255)
                                                                        tACACCTGCGGGTGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGTGGCGTTgg                                                                                                                                            <  7:190372/84‑1 (MQ=255)
                                                                        tACACCTGCGGGTGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGTGGCGTTgg                                                                                                                                            >  8:190372/1‑84 (MQ=255)
                                                                        tACACCTGCGGGTGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGTGGCGTTgg                                                                                                                                            >  2:446064/1‑84 (MQ=255)
                                                                        tACACCTGCGGGTGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGTGGCGTTgg                                                                                                                                            <  1:446064/84‑1 (MQ=255)
                                                                           acCTGCGGGTGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCggg                                                                         >  4:405260/1‑148 (MQ=255)
                                                                           acCTGCGGGTGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCggg                                                                         <  3:405260/148‑1 (MQ=255)
                                                                           acCTGCGGGTGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCGGGa                                                                        <  4:115401/149‑1 (MQ=255)
                                                                             cTGCGGGTGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCGGGACa                                                                      >  5:191860/1‑149 (MQ=255)
                                                                             cTGCGGGTGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCGGGACa                                                                      >  4:186333/1‑149 (MQ=255)
                                                                                cGGGTGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCGGGACAt                                                                     <  3:186333/147‑1 (MQ=255)
                                                                                   gTGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCGGGACATGgcg                                                                 >  2:244766/1‑148 (MQ=255)
                                                                                    tGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCGGGACATGGCGCa                                                               <  7:95003/149‑1 (MQ=255)
                                                                                     gCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGTGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCGGGACATGGCGCAg                                                              >  3:386392/1‑149 (MQ=255)
                                                                                     gCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCGGGACATGGCGCAg                                                              <  3:22140/149‑1 (MQ=255)
                                                                                     gCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCGGGACATGGCGCAg                                                              >  2:128071/1‑149 (MQ=255)
                                                                                      cacaGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGTGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCGGGACATGGCGCAg                                                              <  4:386392/148‑1 (MQ=255)
                                                                                      cacaGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTACGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCGGGACATGGCGCAg                                                              <  1:190174/148‑1 (MQ=255)
                                                                                          gCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGc                                                                            >  8:275921/1‑130 (MQ=255)
                                                                                          gCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGc                                                                            <  7:275921/130‑1 (MQ=255)
                                                                                          gCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCGGGACATGGCGCAGGGAtt                                                         <  7:9493/149‑1 (MQ=255)
                                                                                            gACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCGGGACATGGCGCAGGGATTaa                                                       <  7:263314/149‑1 (MQ=255)
                                                                                               tactCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCGGGACATGGCGCAGGGATTAAGcc                                                    >  8:293744/1‑149 (MQ=255)
                                                                                                    gatgatCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCGGGACATGGCGCAGGGATTAAGCCTAtcg                                               >  1:401490/1‑149 (MQ=255)
                                                                                                       gatCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCGGGACATGGCGCAGGGATTAAGCCTAtcgtcg                                            >  2:460380/1‑149 (MQ=255)
                                                                                                           gATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCGGGACATGGCGCAGGGATTAAGCCTATCGTCgggg                                         >  1:433735/1‑148 (MQ=255)
                                                                                                                 ctggccAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCGGGACATGGCGCAGGGATTAAGCCTATCGTCGGGGCAGAttt                                  <  3:395257/149‑1 (MQ=255)
                                                                                                                    gccAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCGGGACATGGCGCAGGGATTAAGCCTATCGTCGGGGCAGATTTTaa                               <  1:128071/149‑1 (MQ=255)
                                                                                                                     ccAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCGGGACATGGCGCAGGGATTAAGCCTATCGTCGGGGCAGATTTTAAc                              <  7:293744/149‑1 (MQ=255)
                                                                                                                      cAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCGGGACATGGCGCAggg                                                            <  7:12808/118‑1 (MQ=255)
                                                                                                                      cAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCGGGACATGGCGCAggg                                                            >  8:12808/1‑118 (MQ=255)
                                                                                                                      cAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCGGGACATGGCGCAGGGATTAAGCCTATCGTCGGGGCAGATTTTAACg                             >  4:260862/1‑149 (MQ=255)
                                                                                                                      cAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCGGGACATGGCGCAGGGATTAAGCCTATCGTCGGGGCAGATTTTAACg                             >  1:174112/1‑149 (MQ=255)
                                                                                                                            cGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCGGGACATGGCGCAGGGATTAAGCCTATCGTCGGGGCAGATTTTAACGTCCAGt                       >  3:124045/1‑149 (MQ=255)
                                                                                                                             gCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCGGGACATGGCGCAGGGATTAAGCCTATCGTCGGGGCAGATTTTAACGTCCAGTg                      >  1:49957/1‑149 (MQ=255)
                                                                                                                              cACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCGGGACATGGCGCAGGGATTAAGCCTATCGTCGGGGCAGATTTTAACGTCCAGTGc                     <  1:244766/149‑1 (MQ=255)
                                                                                                                                ccGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCGGGACATGGCGCAGGGATTAAGCGTATCGTCGGGGCAGATTTTAACGTCCAGTGCGa                   <  3:354851/149‑1 (MQ=255)
                                                                                                                                   ttGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCGGGACATGGCGCAGGGATTAAGCCTAtcgtcg                                            >  8:201574/1‑121 (MQ=255)
                                                                                                                                   ttGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCGGGACATGGCGCAGGGATTAAGCCTAtcgtcg                                            <  7:201574/121‑1 (MQ=255)
                                                                                                                                     ggTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCGGGACATGGCGCAGGGATTAAGCCTATCGTCGGGGCAGATTTTAACGTCCAGTGCGACctgc              >  3:36102/1‑149 (MQ=255)
                                                                                                                                             aGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCGGGACATGGCGCAGGGATTAAGCCTATCGTCGTGGCAGATTTTAACGTCCAGTGCGACCTGCTGGGggat      <  5:313792/149‑1 (MQ=255)
                                                                                                                                               gcggcggcgTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCGGGACATGGCGCAGGGATTAAGCCTATCGTCGGGGCAGATTTTAACGTCCAGTGCGACCTGCTGGGtgatga    <  4:36102/149‑1 (MQ=255)
                                                                                                                                               gcggcggcgTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCGGGACATGGCGCAGGGATTAAGCCTATCGTCGGGGCAGATTTTAACGTCCAGTGCGACCTGCTGGGtgatga    <  1:42834/149‑1 (MQ=255)
                                                                                                                                               gcggcggcgTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTCGTCTGGTGAAGGTCTACGGAGCGGGACATGGCGCAGGGATTAAGCCTATCGTCGGGGCAGATTTTAACGTCCAGTGCGACCTGCTGGGggatga    <  6:329801/149‑1 (MQ=255)
                                                                                                                                                 ggcggcgTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCGGGACATGGCGCAGGGATTAAGCCTATCGTCGGGGCAGATTTTAACGTCCAGTGCGACCTGCTGGGTGATGAGt  >  1:207696/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGGGACTTGCGTCCTGATTCTTGTGTCGAGATTAAGTAAACCGGAATCTGAAGATGTCTGAACCACGTTTCGTACACCTGCGGGTGCACAGCGACTACTCGATGATCGATGGCCTGGCCAAAACCGCACCGTTGGTAAAAAAGGCGGCGGCGTTGGGTATGCCAGCACTGGCGATCACCGATTTCACCAACCTTTGTGGTCTGGTGAAGTTCTACGGAGCGGGACATGGCGCAGGGATTAAGCCTATCGTCGGGGCAGATTTTAACGTCCAGTGCGACCTGCTGGGTGATGAGT  >  NC_000913/205073‑205366

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