Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,623,709 T→G I104S (ATT→AGT)  yhhH → putative NTF2 fold immunity protein for polymorphic toxin RhsB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,623,7090TG97.9% 808.3 / 1.6 233I104S (ATT→AGT) yhhHputative NTF2 fold immunity protein for polymorphic toxin RhsB
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base G (109/119);  total (111/122)
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

Only 100 of 233 total aligned reads displayed.

GTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACA  >  NC_000913/3623613‑3623803
                                                                                                |                                                                                              
     ggACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCtggtg                                                                                        >  2:3877854/1‑100 (MQ=255)
      gACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCtggtgg                                                                                       >  1:2000908/1‑100 (MQ=255)
      gACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGAGCtggtgg                                                                                       >  1:2782818/1‑100 (MQ=255)
       aCAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTggg                                                                                      >  1:3706449/1‑100 (MQ=255)
        cAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTgggg                                                                                     >  1:169970/1‑100 (MQ=255)
         aGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGt                                                                                    <  2:3618752/100‑1 (MQ=255)
         aGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGt                                                                                    >  2:2109290/1‑100 (MQ=255)
          ggATGCCGCAGAAGAAGAAAAACCATACGAAATCCCTGAGTTAACAACAAGTTAGGTTGTTGCAGGTACCATTCCCTCAGACCAAAGTGCTGGTGGGGtt                                                                                   >  1:5177219/1‑100 (MQ=255)
          ggATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGtt                                                                                   >  1:385390/1‑100 (MQ=255)
          ggATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGtt                                                                                   >  1:107965/1‑100 (MQ=255)
           gATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGCTACCATTCAATCAGACCAAAGTGCTGGTGGGGttt                                                                                  >  2:5428077/1‑100 (MQ=255)
              gCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTtta                                                                               <  1:2970122/100‑1 (MQ=255)
               ccgccgCAGAAGAAAAACCATACGCAATCCCTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTttan                                                                              <  1:4456086/100‑2 (MQ=255)
               ccGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGGTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTttat                                                                              >  1:3997491/1‑100 (MQ=255)
                cGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTttatt                                                                             >  2:2283444/1‑100 (MQ=255)
                cGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTttatt                                                                             >  2:560589/1‑100 (MQ=255)
                cGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTttatt                                                                             <  1:3400706/100‑1 (MQ=255)
                cGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTttatt                                                                             <  2:906879/100‑1 (MQ=255)
                 gCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTttatta                                                                            <  2:273033/100‑1 (MQ=255)
                 gCAGAAGAACAAAAACCATACGATATCACTGAGTTAACAACACGTTGGGTTGTTCAAGGTACCATTCACGCAGACCAAAGTGCTCGTGGGGTTTttatta                                                                            <  1:4844181/100‑1 (MQ=255)
                  cAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTtattat                                                                           >  1:502436/1‑100 (MQ=255)
                  cAGAAGAAGAAAAAACATACGAAATCACTGAGTTAACAACAAGTTTGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTtattat                                                                           <  2:3971477/100‑1 (MQ=255)
                    gaagaagaaAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAg                                                                         >  1:2604120/1‑100 (MQ=255)
                     aagaagaaAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGa                                                                        <  2:4105353/100‑1 (MQ=255)
                      agaagaaAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGaa                                                                       <  2:3660996/100‑1 (MQ=255)
                      agaagaaAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGaa                                                                       >  1:2820419/1‑100 (MQ=255)
                       gaagaaAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCGAAGTGCTGGTGGGGTTTTTATTATAGaaa                                                                      >  1:5084713/1‑100 (MQ=255)
                       gaagaaAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGaaa                                                                      >  2:5304323/1‑100 (MQ=255)
                         agaaAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATa                                                                    <  1:314508/100‑1 (MQ=255)
                          gaaAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAg                                                                   >  1:1867688/1‑100 (MQ=255)
                           aaaaaCCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAgg                                                                  <  2:3811535/100‑1 (MQ=255)
                           aaaaaCCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAgg                                                                  <  1:5413486/100‑1 (MQ=255)
                           aaaaaCCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAgg                                                                  <  2:2127841/100‑1 (MQ=255)
                               aCCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCaaa                                                              >  2:3536029/1‑100 (MQ=255)
                                  aTACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAat                                                           >  1:1515327/1‑100 (MQ=255)
                                  aTACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAat                                                           <  1:3579413/100‑1 (MQ=255)
                                   taagAAAACAATGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTAAAAAAAAAAAAAAAAatg                                                          <  2:1441889/97‑1 (MQ=255)
                                   tACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAatg                                                          <  2:4118783/100‑1 (MQ=255)
                                     cGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAatgat                                                        <  2:2688248/100‑1 (MQ=255)
                                      gAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAatgatg                                                       <  2:2392980/100‑1 (MQ=255)
                                      gAAATCACTCAGTTAACCACAAGTTGGTTTGGTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAatgatg                                                       <  2:4700865/100‑1 (MQ=255)
                                           cACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGgaga                                                  >  2:2916845/1‑100 (MQ=255)
                                                aGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCt                                             <  1:2515609/100‑1 (MQ=255)
                                                 gTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATATGCAAAAATGATGGGAGAATTCTg                                            >  2:322008/1‑100 (MQ=255)
                                                 gTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTg                                            <  1:2380705/100‑1 (MQ=255)
                                                   tAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGaa                                          <  1:567254/100‑1 (MQ=255)
                                                   tAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGaa                                          <  2:4924923/100‑1 (MQ=255)
                                                    aacaacAAGTTGGGCTGTTGAAGGTACCATTCACTCAGACCAAAGCGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAAt                                         >  2:3081628/1‑100 (MQ=255)
                                                     acaacaAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAAtt                                        >  2:3764084/1‑100 (MQ=255)
                                                      caacaaGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAAttt                                       >  2:720257/1‑100 (MQ=255)
                                                       aacaaGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAAtttt                                      >  2:4857590/1‑100 (MQ=255)
                                                       aacaaGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAAtttt                                      <  1:620769/100‑1 (MQ=255)
                                                       aacaaGTTGGGTTGTTGAAGGTACCATACACTCAGACGAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAAtttt                                      <  2:2855955/100‑1 (MQ=255)
                                                        acaaGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTg                                     >  2:1292870/1‑100 (MQ=255)
                                                        acaaGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTg                                     <  2:2990584/100‑1 (MQ=255)
                                                        acaaGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGAAAAAATGATGGGAGAATTCTGAATTTTg                                     <  2:428708/100‑1 (MQ=255)
                                                         caaGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTgg                                    <  1:3913604/100‑1 (MQ=255)
                                                         caaGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTAAAAAAAAAAAAAAAAATGATGGGAGAATTCTGAATTTTgg                                    <  2:5014419/100‑1 (MQ=255)
                                                          aaGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGc                                   <  1:533507/100‑1 (MQ=255)
                                                           aGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGcc                                  <  2:3156194/100‑1 (MQ=255)
                                                           aGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGcc                                  <  1:227368/100‑1 (MQ=255)
                                                            gttgggttggTGAAGGGACCATTCACTCAGACCAAAGTGCTGGTGTTGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCa                                 <  1:2657814/100‑1 (MQ=255)
                                                            gTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCa                                 >  1:4769345/1‑100 (MQ=255)
                                                               gggTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGGAAAAATGATGGGAGAATTCTGAATTTTGGGCACgg                              >  2:5471651/1‑100 (MQ=255)
                                                               gggTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACgg                              >  2:1680657/1‑100 (MQ=255)
                                                               gggTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACgg                              <  2:4467348/100‑1 (MQ=255)
                                                               gggTTGTTGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACgg                              <  1:3021536/100‑1 (MQ=255)
                                                                ggTTGTTGAAGGTACCATTCACTCAGACCCAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGa                             >  2:2483978/1‑100 (MQ=255)
                                                                 gttgttGAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGaa                            >  1:4868828/1‑100 (MQ=255)
                                                                    gttgAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAAt                         >  2:1797813/1‑100 (MQ=255)
                                                                    gttgAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAAt                         >  1:991541/1‑100 (MQ=255)
                                                                    gttgAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAAt                         <  1:5481959/100‑1 (MQ=255)
                                                                    gttgAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAAt                         >  2:3680279/1‑100 (MQ=255)
                                                                    gttgAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAAGGATGGGACAATTCTGAATTTTGGCCACggccaat                         >  1:4564338/1‑95 (MQ=255)
                                                                      tgAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATaa                       >  2:961843/1‑100 (MQ=255)
                                                                      tgAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATaa                       >  2:3040679/1‑100 (MQ=255)
                                                                      tgAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATaa                       >  2:2416766/1‑100 (MQ=255)
                                                                      tgAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATaa                       >  2:1650994/1‑100 (MQ=255)
                                                                       gAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAg                      <  1:4615142/100‑1 (MQ=255)
                                                                       gAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAACAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAg                      >  2:4249187/1‑100 (MQ=255)
                                                                       gAAGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAg                      >  2:3122483/1‑100 (MQ=255)
                                                                        aaGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGt                     <  1:75876/100‑1 (MQ=255)
                                                                         aGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTa                    >  1:1528378/1‑100 (MQ=255)
                                                                         aGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTa                    <  1:2319113/100‑1 (MQ=255)
                                                                         aGGTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTa                    <  1:5508228/100‑1 (MQ=255)
                                                                          ggTACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAAAAAGTAt                   >  2:3283824/1‑100 (MQ=255)
                                                                            tACCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGt                 <  2:1919297/100‑1 (MQ=255)
                                                                             aCCATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTa                >  1:1087168/1‑100 (MQ=255)
                                                                               cATTCACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAg              <  1:120180/100‑1 (MQ=255)
                                                                                   cACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATa          >  2:4795553/1‑100 (MQ=255)
                                                                                   cACTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATa          >  2:2373424/1‑100 (MQ=255)
                                                                                    aCTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATaa         >  1:3549319/1‑100 (MQ=255)
                                                                                     cTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAAt        >  2:1480732/1‑100 (MQ=255)
                                                                                     cTCAGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGAAAGGATAAt        >  1:853096/1‑100 (MQ=255)
                                                                                     ataaGACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTTAGCATAAt        <  2:846850/97‑1 (MQ=255)
                                                                                      acAGATCAAAGTGTTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATc       <  2:5124797/99‑1 (MQ=255)
                                                                                         gACCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCaaa    >  2:2237699/1‑100 (MQ=255)
                                                                                          aCCAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAAc   <  2:4777080/100‑1 (MQ=255)
                                                                                           ccAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACa  >  1:1570004/1‑100 (MQ=255)
                                                                                           ccAAAGTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACa  >  1:1084132/1‑100 (MQ=255)
                                                                                                |                                                                                              
GTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACA  >  NC_000913/3623613‑3623803

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: ‑