Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,891,841 A→G 15.2% G465G (GGT→GGC) ‡ yoaA ← putative ATP‑dependent helicase, DinG family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,891,8410AG15.2% 280.9 / 36.4 105G465G (GGT→GGC) ‡yoaAputative ATP‑dependent helicase, DinG family
Reads supporting (aligned to +/- strand):  ref base A (45/44);  new base G (9/7);  total (54/51)
Fisher's exact test for biased strand distribution p-value = 7.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.76e-02

Only 100 of 107 total aligned reads displayed.

AGCGCATTACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAACAGCAACGAC  >  NC_000913/1891695‑1891982
                                                                                                                                                  |                                                                                                                                             
agCGCATTACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGAc                                                                                                                                              <  3:997570/148‑1 (MQ=255)
     aTTACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGGTCgttg                                                                                                                                         >  3:997127/1‑148 (MQ=255)
      ttACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCgttg                                                                                                                                         >  1:679511/1‑147 (MQ=255)
      ttACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCgttg                                                                                                                                         <  2:679511/147‑1 (MQ=255)
           ggCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCt                                                                                                                                   <  1:752376/148‑1 (MQ=255)
               cTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGAc                                                                                                                                              <  3:240309/133‑1 (MQ=255)
               cTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGAc                                                                                                                                              >  4:240309/1‑133 (MQ=255)
                                       ccTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGGTCgttgtttt                                                                                                                                     <  4:1117715/118‑3 (MQ=255)
                                       ccTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGGTCgttgtttt                                                                                                                                     >  3:1117715/1‑116 (MQ=255)
                                       ccTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCga                                                                                                                               <  2:257328/124‑1 (MQ=255)
                                       ccTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCga                                                                                                                               >  1:257328/1‑124 (MQ=255)
                                       ccTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCgcagca                                                                                                                <  4:913301/139‑1 (MQ=255)
                                       ccTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCgcagca                                                                                                                >  3:913301/1‑139 (MQ=255)
                                        cTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTaa                                                                                                      >  1:428007/1‑148 (MQ=255)
                                        cTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTaa                                                                                                      >  2:388604/1‑148 (MQ=255)
                                        cTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCAGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGGTCgttgttt                                                                                                                                      >  3:736689/1‑115 (MQ=255)
                                        cTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCAGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGGTCgttgttt                                                                                                                                      <  4:736689/116‑2 (MQ=255)
                                          ttGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGGTCgttgtttta                                                                                                                                    <  1:50761/116‑4 (MQ=255)
                                          ttGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGGTCgttgtttta                                                                                                                                    >  2:50761/1‑113 (MQ=255)
                                            gCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCtt                                                                                                                                  <  2:314867/116‑1 (MQ=255)
                                            gCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCtt                                                                                                                                  >  1:314867/1‑116 (MQ=255)
                                                 tttCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCgatgat                                                                                                                           >  4:311302/1‑118 (MQ=255)
                                                 tttCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCgatgat                                                                                                                           <  3:311302/118‑1 (MQ=255)
                                                  ttCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTgg                                                                                                  >  2:623512/1‑142 (MQ=255)
                                                       ccTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCt                                                                                                                                   >  4:1119012/1‑104 (MQ=255)
                                                       ccTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCt                                                                                                                                   <  3:1119012/104‑1 (MQ=255)
                                                       ccTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAAcc                                                                                       >  1:126618/1‑148 (MQ=255)
                                                         tGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCgatgat                                                                                                                           <  1:1073850/110‑1 (MQ=255)
                                                         tGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCgatgat                                                                                                                           >  2:1073850/1‑110 (MQ=255)
                                                         tGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACctg                                                                                     <  2:932426/148‑1 (MQ=255)
                                                           caacaaTACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACctggc                                                                                   >  1:805635/1‑148 (MQ=255)
                                                             acaaTACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGa                                                                                                                                               >  3:506178/1‑86 (MQ=255)
                                                             acaaTACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGa                                                                                                                                               <  4:506178/86‑1 (MQ=255)
                                                               aaTACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCtggt                                                                               <  1:388604/148‑1 (MQ=255)
                                                                aTACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCgttgt                                                                                                                                        <  2:771620/90‑1 (MQ=255)
                                                                aTACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCgttgt                                                                                                                                        >  1:771620/1‑90 (MQ=255)
                                                                 tACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGGTCgttgtt                                                                                                                                       <  3:28917/90‑1 (MQ=255)
                                                                 tACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGGTCgttgtt                                                                                                                                       >  4:28917/1‑90 (MQ=255)
                                                                    aGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCtggt                                                                               <  1:1113057/143‑1 (MQ=255)
                                                                    aGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCtggt                                                                               >  2:1113057/1‑143 (MQ=255)
                                                                     ggAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTa                                                                                                       <  4:651567/118‑1 (MQ=255)
                                                                     ggAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGttt                                                                         <  2:805635/148‑1 (MQ=255)
                                                                     ggAAGCGTCAGGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTa                                                                                                       >  3:651567/1‑118 (MQ=255)
                                                                      gAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCgc                                                                                                                    <  1:802852/104‑1 (MQ=255)
                                                                      gAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCgc                                                                                                                    >  2:802852/1‑104 (MQ=255)
                                                                      gAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTg                                                                        <  2:428007/148‑1 (MQ=255)
                                                                       aaGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTg                                                                        >  1:484108/1‑147 (MQ=255)
                                                                            tCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGGTCgttgtttta                                                                                                                                    <  1:752452/82‑4 (MQ=255)
                                                                            tCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGGTCgttgtttta                                                                                                                                    >  2:752452/1‑79 (MQ=255)
                                                                            tCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAAc                                                                                        >  2:898139/1‑126 (MQ=255)
                                                                            tCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAAc                                                                                        <  1:898139/126‑1 (MQ=255)
                                                                            tCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCtggt                                                                               <  4:296107/135‑1 (MQ=255)
                                                                            tCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCtggt                                                                               >  3:296107/1‑135 (MQ=255)
                                                                                              gCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACagag                                                >  1:705849/1‑148 (MQ=255)
                                                                                                tCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACctggctgg                                                                                <  2:65408/114‑1 (MQ=255)
                                                                                                tCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACctggctgg                                                                                >  1:65408/1‑114 (MQ=255)
                                                                                                tCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTa                                              >  4:191134/1‑148 (MQ=255)
                                                                                                 cGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTgg                                                                                                  >  3:1133679/1‑95 (MQ=255)
                                                                                                 cGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTgg                                                                                                  <  4:1133679/95‑1 (MQ=255)
                                                                                                    ccAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACctggct                                                                                  >  1:816847/1‑108 (MQ=255)
                                                                                                    ccAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACctggct                                                                                  <  2:816847/108‑1 (MQ=255)
                                                                                                     cAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCa                                                                   >  4:635838/1‑122 (MQ=255)
                                                                                                        aTCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCgatga                                                                                                                            >  1:29361/1‑62 (MQ=255)
                                                                                                        aTCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCgatga                                                                                                                            <  2:29361/62‑1 (MQ=255)
                                                                                                         tCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTggc                                     >  3:704693/1‑148 (MQ=255)
                                                                                                           gcgcATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTggcgg                                   <  4:704693/148‑1 (MQ=255)
                                                                                                            cgcATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTggcggc                                  >  1:180849/1‑148 (MQ=255)
                                                                                                                tcatGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTa                              <  2:180849/148‑1 (MQ=255)
                                                                                                                 catGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGAt                                                                >  4:1125314/1‑113 (MQ=255)
                                                                                                                 catGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGAt                                                                <  3:1125314/113‑1 (MQ=255)
                                                                                                                   tGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGaa                                                                                         >  2:34072/1‑86 (MQ=255)
                                                                                                                   tGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGaa                                                                                         <  1:34072/86‑1 (MQ=255)
                                                                                                                      cGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCaaa                        <  2:126618/148‑1 (MQ=255)
                                                                                                                         gCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCg                                                                                                <  2:1108180/73‑1 (MQ=255)
                                                                                                                         gCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCg                                                                                                >  1:1108180/1‑73 (MQ=255)
                                                                                                                              gTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTg                                                                        <  2:848842/92‑1 (MQ=255)
                                                                                                                              gTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTg                                                                        >  1:848842/1‑92 (MQ=255)
                                                                                                                                  aaaGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTg                                       <  3:819653/121‑1 (MQ=255)
                                                                                                                                  aaaGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTg                                       >  4:819653/1‑121 (MQ=255)
                                                                                                                                    aGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCt                                 >  3:276015/1‑125 (MQ=255)
                                                                                                                                    aGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAACa          >  3:1097281/1‑148 (MQ=255)
                                                                                                                                         ctaacaacGGTCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGTCTGTAATCAAAtgggct                  >  3:933325/3‑135 (MQ=255)
                                                                                                                                         ctaacaacGGTCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAAtgggct                  <  4:933325/133‑1 (MQ=255)
                                                                                                                                         aaaaCAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAATAGCAAc     >  1:838915/1‑148 (MQ=255)
                                                                                                                                          aaaCAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCaa                         <  4:603918/127‑1 (MQ=255)
                                                                                                                                          aaaCAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCaa                         >  3:603918/1‑127 (MQ=255)
                                                                                                                                          aaaCAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAATAGCAACg    <  2:838915/148‑1 (MQ=255)
                                                                                                                                           aacaacGGTCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAACAGCAACGa   >  1:722618/1‑148 (MQ=255)
                                                                                                                                           aacaacGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAACAGCAACGa   <  2:172207/148‑1 (MQ=255)
                                                                                                                                           aacaacGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAACAGCAACGa   >  1:172207/1‑148 (MQ=255)
                                                                                                                                            acaacGGTCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAACAGCAACGAc  >  4:304751/1‑148 (MQ=255)
                                                                                                                                            acaacGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCgca                                                          >  2:797395/1‑92 (MQ=255)
                                                                                                                                            acaacGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCgca                                                          <  1:797395/92‑1 (MQ=255)
                                                                                                                                            acaacGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCaa                         >  2:511316/1‑125 (MQ=255)
                                                                                                                                            acaacGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCaa                         <  1:511316/125‑1 (MQ=255)
                                                                                                                                             caacGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTa                                                                                                       >  1:961292/1‑46 (MQ=255)
                                                                                                                                             caacGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTa                                                                                                       <  2:961292/46‑1 (MQ=255)
                                                                                                                                             caacGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTggcggc                                  <  1:957761/115‑1 (MQ=255)
                                                                                                                                               acGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCaa                         >  4:302880/1‑122 (MQ=255)
                                                                                                                                               acGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCaa                         <  3:302880/122‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
AGCGCATTACCGGCGCTGACAAATTGCTGCAACAGTTGCCCTTTGCTGGTTTCCCCCTGCAACAATACAGGAAGCGTCATGGTAGCGCGGAACTGCTCGGCCAGATCGCGCATCATGGCGTGCGAGGTACAAAGCATAAAACAACGACCGTTGTTAGCTTCGATGATCGGTCGCAGCATTGCCGCTAACTGGCGAGCAGAACCTGGCTGGTTGGTTTGCGGCAGATTGCGCAGCACACAGAGTAACGCCTGGCGGCTGTAATCAAATGGGCTGGGCAACAGCAACGAC  >  NC_000913/1891695‑1891982

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