Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,623,820 C→T 5.5% intergenic (+38/‑67) yhhH → / → yrhC putative NTF2 fold immunity protein for polymorphic toxin RhsB/pseudogene fragment

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,623,8200CT 5.5% 185.0 / 2.9 74intergenic (+38/‑67)yhhH/yrhCputative NTF2 fold immunity protein for polymorphic toxin RhsB/pseudogene fragment
Reads supporting (aligned to +/- strand):  ref base C (40/27);  new base T (2/2);  total (44/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

GGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGACTCCAGACAGTCAAACGATAGAAAAAGATAGCCTTTATGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAAT  >  NC_000913/3623677‑3623968
                                                                                                                                               |                                                                                                                                                    
ggTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTCTAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGCAGAGt                                                                                                                                                 <  1:10921/149‑1 (MQ=255)
 gttgttGAAGGTACCATTCCCTCAGACCAAATTGCTGGTGGGGTTTTTTTTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGTCCACGGAAAATAAGTATGTATTGATAATCAAACATCGTGTGCGCTGATGGCAGAGtt                                                                                                                                                >  2:196122/1‑148 (MQ=255)
 gttgttGAAGGTACCATTCACTCAGGCCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGCAGAGTa                                                                                                                                                >  1:518611/1‑149 (MQ=255)
  tagatgAAGGTACAATTCACTCAGAACAAATTGCTGGTGGGGATTTTATTATAGAAAAAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATAGTGTGCGCTGATGTCAGAGTac                                                                                                                                               >  2:557018/5‑148 (MQ=255)
    gttgAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGt                                                                                                                                             >  2:375053/1‑149 (MQ=255)
    gttgAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGt                                                                                                                                             <  2:635809/149‑1 (MQ=255)
       gAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGTgga                                                                                                                                          >  1:847307/1‑149 (MQ=255)
            tACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGTGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCTAACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGACt                                                                                                                                     <  1:583750/149‑1 (MQ=255)
             aCCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGACTc                                                                                                                                    >  2:862478/1‑149 (MQ=255)
              cgaTTCACTCAGACCAAATTGCTGGTTGGGTTTTTATAATAGAATTAGGCAAATATGATGGGAGAAATCTGAATTTTGGCGACGGTAAATTAGTATGTAAGCATATTCATACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGACTcc                                                                                                                                   <  1:522520/147‑1 (MQ=255)
              ccATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGACTcc                                                                                                                                   <  1:342598/149‑1 (MQ=255)
                aTTCACTCAGACCAAATTGCTGGGGGGGTTAATATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAAATTTGGCCACGGATAATAAGTATCTAAGGAAAATCAAACATCGTGGGCGCTGATGGAAGAGTAGGTGTAGGACACCaa                                                                                                                                 >  2:946299/1‑148 (MQ=255)
                 ttCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGACTCCAGa                                                                                                                                <  1:32374/149‑1 (MQ=255)
                 ttCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGACTCCAGa                                                                                                                                >  2:54554/1‑149 (MQ=255)
                   cACTCAGCCCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGACTCCAGACa                                                                                                                              <  2:895034/149‑1 (MQ=255)
                     cTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGCTGGCACAGTAGGTGGAGGACTCCAGACAGt                                                                                                                            >  2:917552/1‑149 (MQ=255)
                       cAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGACTCCAGACAGTCa                                                                                                                          <  2:183553/149‑1 (MQ=255)
                         gACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGACTCCAGACAGTCaaa                                                                                                                        <  1:823783/149‑1 (MQ=255)
                          aCCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGACTCCAGACAGTCAAAc                                                                                                                       >  1:948745/1‑149 (MQ=255)
                            cAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGTATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGAAGAGTAGGTGGAGGACTCCAGACAGTAAAActa                                                                                                                     >  1:993344/1‑147 (MQ=255)
                             aaaTTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCTTGTGCACTGATGGCAGAGTAGGTGGAGGACTCCAGACAGTCAAACGAt                                                                                                                    >  2:19744/1‑149 (MQ=255)
                              aaTTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGTAGAGTAGGTGGAGTACTCCAGACAGTCAAACGATa                                                                                                                   >  1:812538/1‑149 (MQ=255)
                              aaTTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGTAGAGTAGGTGGAGGACTCCAGACAGTCAAACGATa                                                                                                                   >  1:373468/1‑149 (MQ=255)
                                  gCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGACTCCAGACAGTCAAACGATAGaaa                                                                                                               >  2:505079/1‑149 (MQ=255)
                                     tgtggtgTTTTTATTATATAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGTATAATCAAACATCGTTTGCGCTGATGGCAGAGTAGGTGGAGGACTCCAGACAGTCAAACGATAGAAAAAg                                                                                                            <  1:175069/148‑1 (MQ=255)
                                     ggtggGGTTTTTATTATAGAAATAGGAAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGACTCCAGACAGTCAAACGATAGAAAAAg                                                                                                            <  1:144642/149‑1 (MQ=255)
                                         gggTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGACTCCAGACAGTCAAACGATAGAAAAAGATac                                                                                                        >  2:408093/1‑148 (MQ=255)
                                         gggTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGACTCAAGAAAGTCAAACGATAGAAAAAGATAg                                                                                                        >  1:457149/1‑149 (MQ=255)
                                          ggTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTTGCCACGGAAAATAAGTATGTAAGTATAATCAAACATCGTGTGCGCTGATGGCATAGTAGGTGGAGGACTCCAGACAGTCAAACGATAGAAAAAGATAGc                                                                                                       <  2:190946/149‑1 (MQ=255)
                                           gTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGTTGGAGGACTCCAGACAGTCAAACGATAGAAAAAGATAGca                                                                                                      >  1:823101/1‑148 (MQ=255)
                                                   tataGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAATTAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGACTCCAGAAAGTCAAACGATAGAAAAAGATAGCCTTTATGGa                                                                                              >  1:803694/1‑149 (MQ=255)
                                                   taaaGAAATAGGCAAAAATGCTGGCAGAAGTCTGAATTTTGCCCACGGAAAATAAGTATGTAAGGATAATAAAACATCGTGTGCGCTGAGGGCAGAGTAGGTGGCGGACTCCAGACAGTCAAACGATAGAAAAAGATAGCCTTTATGGa                                                                                              <  2:110361/146‑1 (MQ=255)
                                                        aaaTAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCGACGGAAAATAAGTATGTTAGGATAATCAAACATAGTGTGCGCGGATGGCAGAGTAGGTGGAGGACTCGAGACAGTCAAACGATAGAAAAAGATTGCCTTGATGGAGGTTc                                                                                         >  2:257682/1‑149 (MQ=255)
                                                        aaaTAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGACTCCAGACAGTCAAACGATAGAAAAAGATAGCCTTTATGGAGGTTc                                                                                         <  1:753612/149‑1 (MQ=255)
                                                           tAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGGGGAGGACTCCAGACAGTCAAACGATAGAAAAAGATAGCCTTTATGGAGGTTCCTg                                                                                      <  1:54554/149‑1 (MQ=255)
                                                             ggCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGCCTCCAGACAGTCAACCCATAGAAAAAGATAGCCTTTATGGAGGTTCCTGCa                                                                                    >  2:673738/1‑149 (MQ=255)
                                                              gCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGAGTCCAGACAGTCAAACGATAGTAATAGATCGCCTGGTTGGAGGTTCCTGCaa                                                                                   >  2:99043/1‑149 (MQ=255)
                                                                aaaaaTGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGCATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGGGGAGGACTCCAGACAGTCAAACGAGAGAAAAAGATAGCCTTTATGGAGGTTCCTGCAATg                                                                                 >  1:916657/1‑149 (MQ=255)
                                                                  aaaTGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGAAGAGTAGGTGGAGGACTCCAGACAGTCAAACGATAGAAAATGATAGCCTTTATGGAGGTTCCTGCAATGTc                                                                               <  1:594728/149‑1 (MQ=255)
                                                                     tgatgGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGACTCCAGACAGTCAAACGATAGAAAAAGATAGCCTTTATGGAGGTTCCTGCAATGTCaaa                                                                            <  1:155094/149‑1 (MQ=255)
                                                                     tgatgGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGACTCCAGACAGTCAAACGATAGAAAAAGATAGCCTTTATGGAGGTTCCTGCAATGTCaaa                                                                            <  1:546175/149‑1 (MQ=255)
                                                                     tgatgGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAAGCCAACATCGTGTGCGCGGATGGCAGAGTAGGTGGAGGACTCCAGACAGTCAAACGATAGAAAAAGATAGCCTTTATGGAGGTTCCTGCAATGTCaaa                                                                            <  1:896294/149‑1 (MQ=255)
                                                                         gggAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGACTCCAGACAGTCAAACGATAGAAAAAGATAGCCTTTAGGGAGGTTCCTGCAATGTCAAataca                                                                        >  2:521010/1‑149 (MQ=255)
                                                                           gagaATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGTTAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGACTCCAGACAGTCAAACGATAGAAAAAGATAGCCTTTATGGAGGTTCCTGCAATGTCAAatacata                                                                      >  2:276362/1‑149 (MQ=255)
                                                                           gagaATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGACTCCAGACAGTCAAACGATAGAAAAAGATAGCCTTTATGGAGGTTCCTGCAATGTCAAatacata                                                                      >  2:825720/1‑149 (MQ=255)
                                                                            agaATTCTGAATTTTGGCCACGTAAACTAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGACTCCAGACAGTCAAACGATAGAAAAAGATAGCCTTTATGGAGGTTCCTGCAATGTCAAatacatac                                                                     <  1:696125/149‑1 (MQ=255)
                                                                                 tCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGACTCCAGACAGTCAAACGATAGAAAAAGATAGCCTTTATGGAGGTTCCTGCAATGTCAAATACATACCagaa                                                                >  1:370568/1‑149 (MQ=255)
                                                                                   tGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGACTCCAGACAGTCAAACGATAGAAAAAGATAGCCTTTATGGAGGTTCCTGCAATGTCAAATACATACCagaaaa                                                              <  2:518611/149‑1 (MQ=255)
                                                                                    gAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGACTCCAGACAGTCAAACGATAGAAAAAGATAGCCTTTATGGAGGTTCCTGCAATGTCAAATACAAACCagaaaac                                                             >  2:146990/1‑148 (MQ=255)
                                                                                      aTTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGACGCCAGACAGTCAAACGATAGAAAAAGATAGCCTTTATGGAGGTTCCTGCAATGTCAAATACATACCagaaaagaa                                                           <  2:582592/149‑1 (MQ=255)
                                                                                       ttttGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGACTCCAGACAGTCAAACGATAGAAAAAGATAGCCTTTATGGAGGTTCCTGCAATGTCAAATACATACCagaaaagaaa                                                          >  1:113191/1‑149 (MQ=255)
                                                                                                   aaaaTAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGACTCCAGACAGTCAAACGATAGAAAAAGATAGCCTTTATGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAaga                                              >  1:626286/1‑149 (MQ=255)
                                                                                                      aTAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGACTCCAGACAGTCAAACGATAGAAAAAGATAGCCTTTATGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTa                                           >  2:49953/1‑149 (MQ=255)
                                                                                                      aTAAGTATGTAAGGATAATCAAACATAGTGTGCGATGATGGCAGAGTAGGTGGAGGACTCCAGACAGTCAAACGATAGAAAAAGATAGCCTTTATGGAGGTTACTGCAAAGACAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTa                                           >  2:835187/1‑149 (MQ=25)
                                                                                                       tAAGTATGTAAGGATAATCAAACTTCGTGTGCGCTGATGTCAGAGTAGGTGGAGTACTCCAGACAGTCAAACGTTATAAAAAGATAGGCTTTATGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGTGTAt                                          <  1:557018/149‑1 (MQ=21)
                                                                                                       tAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGTCAGAGTAGGTTGAGGACTCCAGACAGTCAAACGATAGAAAAAGATAGCCTTTATGGAGGTTCCTGCAAGGTCAAATACAAACCACAAAACCAAGGCAAGTAAagctgat                                          >  2:671468/1‑144 (MQ=25)
                                                                                                          gtaagtaaGGATAATCAATCATCGTGTGCGCTGAAGGCACAGTAGGTGGAGGACTCCAGACAGACAAACCATAGAAAAACAAAGCCTATATAGAGGTTCCTCCAAAGACAAATACATACCAGAAAAGAAAGGCAAGTAAAGATTaacgt                                       >  1:327703/1‑145 (MQ=21)
                                                                                                          gTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGACTCCAGACAGTCAAACGATAGAAAAAGATAGACTTTATGGAGTTTCCTGCAATGTCAAAAACATACCAGAAAATAAAGCCAAGTAATGAGAATgtt                                       >  1:514084/1‑147 (MQ=25)
                                                                                                              gTAAGGATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGACTCCAGACAGTCAAACGATAGAAAAAGATAGCCTTTATGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTtat                                   >  2:104853/1‑149 (MQ=255)
                                                                                                                  ggATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGACTCCAGACAGTCAAACGATAGAAAAAGATAGCCTTTATGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATataa                               <  1:375053/149‑1 (MQ=255)
                                                                                                                   gATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGACTCCAGACAGTCAAACGATAGAAAAAGATAGCCTTTATGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATataaa                              >  2:523622/1‑148 (MQ=255)
                                                                                                                       aTCAAACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGACTCCAGACAGTCAAACGATAGAAAAAGATAGCCTTTATGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATAAAt                          <  2:111369/149‑1 (MQ=255)
                                                                                                                        taaaaCATCGTGGGCGCCGATGCCAGTGTAGGTGGAGGACACCAGACACTAAAACGATAGAAAAAGATAGCCTTACAGGAGGTTCCTGCAATGTCAAATACATACCAGACAACAAAGGCAAGAAAAGAGTTTGGTTTATATAATAAAtt                         <  2:266812/147‑2 (MQ=17)
                                                                                                                        tCAAACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGACTACAGACAGTCAAACGATAGAAAAAGATAGCCTTTATGGAGGTTCCTGCAATGTCAAATACAAACAAGAAAAGAAAGGCAAGTAAAGAGTATCGTTTATATAATAAATg                         >  1:1004972/1‑149 (MQ=255)
                                                                                                                         cAAACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGACTCCAGACAGTCAAACGATAGAAAAAGATAGCCTTTATGGAGGTTCCTGCAATGGCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATAAATGt                        <  1:570908/149‑1 (MQ=255)
                                                                                                                          aaaCATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGACTCCAGACAGTCAAACGATAGAAAAAGATAGCCTTTATGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATCGTTTATATAATAAATGTa                       >  1:298528/1‑149 (MQ=255)
                                                                                                                               tCGTGTGCGCTGATGGCAGAGTAGGGTGAGGACTCCAGACAGTCAAACGATAGAAAAAGATAGCCTTTATGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGaaa                  <  2:546735/149‑1 (MQ=255)
                                                                                                                                 gtgtGCGCTGATGGCAGATTAGGTGGAGGACTCCAGACAGTCAAACGATAGAAAAAGATAGCCTTTATGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAATTAAAGAGTATGGTTTATATAATAAATTTAAGAAAca                >  2:771791/1‑148 (MQ=255)
                                                                                                                                     gcgcTGATGGCAGAGTAGGTGGAGGACTCCAGACAGTCAAACGATAGAAAAAGATAGCCTTTATGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGATTATATAATAAATGTAAGAAACTAAat            >  1:601825/1‑149 (MQ=255)
                                                                                                                                      cgcTGATGGTAGAGTAGGTGGAGGACTCCAGACAGTCAAACGATAGAAAAAGATAGCCTTTATGGAGGTTCCTGCAATGTCAATTACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAatg           <  2:812538/149‑1 (MQ=255)
                                                                                                                                      cgcTGATGGTAGAGTAGGTGGAGGACTCCAGACAGTCAAACGATAGAAAAAGATAGCCTTTATGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAatg           <  2:373468/149‑1 (MQ=255)
                                                                                                                                      cgcTGATGGCAGAGTAGGTGGAGGACTCCAGACAGTCAAACGATAGAAAAAGATAGCCTTTATGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAatg           >  1:321111/1‑149 (MQ=255)
                                                                                                                                       gcTGATGGCAGAGTAGGTGGAGGACTCCAGACAGTCAAACGATAGAAAAAGATAGCCTTTATGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAatga          >  1:772969/1‑149 (MQ=255)
                                                                                                                                         tGATGGCAGAGTAGGTGGAGGACTCCAGACAGTCAAACCATAGATAAAGATAGCCTTTATGGAGTTTCCTGCAATGGCAAATACATACCAGAAAAGAAAGGCAATTTAAGTGTATGGTTTATATAATAAATGTAAGATACTAAattatt        >  2:744185/1‑148 (MQ=17)
                                                                                                                                               cAGAGTAGGTTGAGGACTCCAGACAGTCAAACGATAGAAAAAGATAGCCTTTATGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAACGGCATGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATTATGAat  <  2:948745/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
GGTTGTTGAAGGTACCATTCACTCAGACCAAATTGCTGGTGGGGTTTTTATTATAGAAATAGGCAAAAATGATGGGAGAATTCTGAATTTTGGCCACGGAAAATAAGTATGTAAGGATAATCAAACATCGTGTGCGCTGATGGCAGAGTAGGTGGAGGACTCCAGACAGTCAAACGATAGAAAAAGATAGCCTTTATGGAGGTTCCTGCAATGTCAAATACATACCAGAAAAGAAAGGCAAGTAAAGAGTATGGTTTATATAATAAATGTAAGAAACTAAATGATGATGAAT  >  NC_000913/3623677‑3623968

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