Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,648,319 C→T 5.2% noncoding (26/84 nt) agrB → antisense sRNA antitoxin for the DinQ toxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,648,3190CT 5.2% 343.3 / 3.4 115noncoding (26/84 nt)agrBantisense sRNA antitoxin for the DinQ toxin
Reads supporting (aligned to +/- strand):  ref base C (69/38);  new base T (3/3);  total (73/42)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 125 total aligned reads displayed.

TATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATA  >  NC_000913/3648170‑3648467
                                                                                                                                                     |                                                                                                                                                    
tatGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTccagc                                                                                                                                                      >  1:1382991/1‑150 (MQ=255)
tatGACGATGACCCAGTGAAAGTATATAAATCGTCACTGAGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTccagc                                                                                                                                                      >  2:1387357/1‑150 (MQ=255)
 atGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCATGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAATAACTATCTCACTccagcc                                                                                                                                                     <  2:362838/150‑1 (MQ=255)
   gACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTccagccag                                                                                                                                                   >  2:150112/1‑150 (MQ=255)
   gACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTccagccag                                                                                                                                                   >  2:1589809/1‑150 (MQ=255)
           cccAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGccacca                                                                                                                                           >  2:666975/1‑150 (MQ=255)
           cccAGTGAAAGTATATAAATCGTCACTGAGATATATACAGAAGTGCTCCCTACTACAGCTGAAGAAATCTCTAATTCTTGAAAAGTTAGCCACTGGCTTATAGTATTGAGCTTTTAGATAAGAACTCTCTCACTCCATCAATAGccaaca                                                                                                                                           >  2:101179/1‑150 (MQ=255)
           accAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAATAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGccacca                                                                                                                                           <  2:686235/149‑1 (MQ=255)
                    aGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGc                                                                                                                                  >  1:1812997/1‑150 (MQ=255)
                          taAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGATAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGACACCAACTCAGGGCTGGaaa                                                                                                                            >  1:1968851/1‑150 (MQ=255)
                            aaTCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGt                                                                                                                          <  2:688819/150‑1 (MQ=255)
                                  cACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGGTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTaaaaaa                                                                                                                    >  2:1273588/1‑150 (MQ=255)
                                    cTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAAcc                                                                                                                  >  1:1670058/1‑150 (MQ=255)
                                    cTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGCTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAAcc                                                                                                                  >  2:470795/1‑150 (MQ=255)
                                      gtgATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGa                                                                                                                >  1:851404/3‑150 (MQ=255)
                                      gCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGa                                                                                                                >  2:1160597/1‑150 (MQ=255)
                                       cGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGAc                                                                                                               >  1:3241/1‑150 (MQ=255)
                                          tatataCCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCa                                                                                                            >  1:262759/1‑150 (MQ=255)
                                           atataCCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCaa                                                                                                           <  1:1533942/150‑1 (MQ=255)
                                              taCCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGCGCTGTTTGATAAGAACTCTCTCACTCCAGCCAGAGCCACGAGCTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGt                                                                                                        >  1:111572/1‑150 (MQ=255)
                                              taCCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATTAGAACTCTCTCACTCCCGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGt                                                                                                        >  1:106827/1‑150 (MQ=255)
                                              taCCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAACCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGt                                                                                                        >  1:494392/1‑150 (MQ=255)
                                               aCCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGCGCCACCCACTCCGGGCTGGAAAGTAAAAAACCGACGCAAAGTc                                                                                                       >  1:135705/1‑150 (MQ=255)
                                                ccGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCg                                                                                                      >  2:498227/1‑150 (MQ=255)
                                                ccGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCg                                                                                                      >  2:368463/1‑150 (MQ=255)
                                                 cGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGTCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGATCCACCAACTCAGGTCTGGAAAGTTAAAAACCGACGCAAAGTCgg                                                                                                     <  2:902503/150‑1 (MQ=255)
                                                      tGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCCCTGGCTAATAGTATTGAGATGTTAGATAAGCACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttt                                                                                                >  1:238756/1‑150 (MQ=255)
                                                      tGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttt                                                                                                <  1:1612441/150‑1 (MQ=255)
                                                       gCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                               <  1:215636/150‑1 (MQ=255)
                                                        ctccctccGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttttt                                                                                              <  2:647794/150‑1 (MQ=255)
                                                           cctccGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACg                                                                                           <  1:1866580/150‑1 (MQ=255)
                                                              ccGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTcc                                                                                        >  2:338345/1‑150 (MQ=255)
                                                               cGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCt                                                                                       >  2:739605/1‑150 (MQ=255)
                                                                 ccAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGa                                                                                     <  2:1935312/150‑1 (MQ=255)
                                                                     cTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGAttaa                                                                                 >  1:265668/1‑148 (MQ=255)
                                                                      tGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGATAGATAAGAACTCTCTCACTCCTGCCAGAGCCACCAACTCAGGGCTGCAAAGTAAAAAACCGACGCAAAGACGGTTTTTTTACGTCCTGCTTAAg                                                                                >  1:313129/1‑150 (MQ=37)
                                                                       ttagaaATAGCTATTTCTTGCAATTTTATCAATTGGGTAATAGTATTGATGCGTTAGATCATTCCTATCTCACTCCAGGCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGa                                                                               <  2:2000296/148‑1 (MQ=11)
                                                                       gaagaaATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGa                                                                               >  1:1257437/1‑150 (MQ=255)
                                                                        aagaaATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGAc                                                                              >  2:638774/1‑150 (MQ=255)
                                                                        aagaaATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGAc                                                                              >  2:397603/1‑150 (MQ=255)
                                                                                tatAATTCTTGAATTGTTATGGAGTTGATAATTGTTTTGTTCTGTTATATAATATCTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTc                                                                      <  1:210506/148‑1 (MQ=17)
                                                                                gcttttgcTTGCAATTTTAGCATTTTGGTATTAGTATTGATGTGTTAGATTATATCTATCTCACTCCAGCGAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTc                                                                      <  1:257947/143‑1 (MQ=14)
                                                                                gCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTc                                                                      >  1:458349/1‑150 (MQ=255)
                                                                                gCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTAACTCCAGCCAGAGACACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTc                                                                      >  1:37065/1‑150 (MQ=255)
                                                                                 cTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCa                                                                     <  2:357160/150‑1 (MQ=255)
                                                                                  tAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGCACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCaa                                                                    >  2:1459237/1‑150 (MQ=255)
                                                                                  tAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCaa                                                                    >  2:1446307/1‑150 (MQ=255)
                                                                                  tAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGGCGGTTTTTTTACGTCCTGATTCAGACCGCCTTTCaa                                                                    >  2:1011890/1‑150 (MQ=255)
                                                                                    aTTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAAt                                                                  <  2:452054/150‑1 (MQ=255)
                                                                                        ttGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAAt                                                              >  1:1659777/1‑150 (MQ=255)
                                                                                            aaTGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGcc                                                          <  1:338345/150‑1 (MQ=255)
                                                                                                  aGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCa                                                    >  1:707410/1‑150 (MQ=255)
                                                                                                   gCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAGCTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTGTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCGAACTCaa                                                   >  2:97013/1‑150 (MQ=255)
                                                                                                   gCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCaa                                                   >  1:1306036/1‑150 (MQ=255)
                                                                                                   gCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCaa                                                   >  1:1996383/1‑150 (MQ=255)
                                                                                                      aCTGGCTAATAGTATTGAGCTGTTAGATAATAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAAt                                                <  1:1893509/150‑1 (MQ=255)
                                                                                                         ggCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTca                                             >  1:1949534/1‑150 (MQ=255)
                                                                                                          tcTAATAGTATTTATCTGTTAGATAAGAACTCTCTCACTCCAGCCTTAGACACCAACTCAGGGCTGGAATGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTTATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTcac                                            <  2:1879801/149‑1 (MQ=37)
                                                                                                               tAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTa                                       >  1:708615/1‑150 (MQ=255)
                                                                                                                 gTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTAtt                                     >  1:1954030/1‑150 (MQ=255)
                                                                                                                  tATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTa                                    <  2:1142744/150‑1 (MQ=255)
                                                                                                                   aTTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTAc                                   >  1:1325246/1‑150 (MQ=255)
                                                                                                                     tgatgtGTTAGATTATTTCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCt                                 <  2:460571/145‑1 (MQ=18)
                                                                                                                      gAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCtt                                <  1:150112/150‑1 (MQ=255)
                                                                                                                          tGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCctc                            <  1:1879119/150‑1 (MQ=255)
                                                                                                                          tGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCctc                            >  2:1363788/1‑150 (MQ=255)
                                                                                                                           gTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCctct                           >  2:1755253/1‑150 (MQ=255)
                                                                                                                           gTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCctct                           >  2:1305014/1‑150 (MQ=255)
                                                                                                                           gTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTg                          <  1:1233315/150‑1 (MQ=255)
                                                                                                                             gaGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGCGATGGAAAGTAAAAAACCGACGCAAAGACTGGTATTTTACGACCTGCTTCACCCCTCCTTTCAAATGAAAAGCCAACTCACCATTCACAACAATAACCTTCCTCTGc                         >  1:1384227/2‑150 (MQ=21)
                                                                                                                              aGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGCCGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACATATTCCCTTCCTCTGCa                        >  2:612042/1‑150 (MQ=255)
                                                                                                                              aGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCACTGCa                        >  2:772051/1‑150 (MQ=255)
                                                                                                                               gATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACCCCTATTACCTTCCTCTGCAc                       >  1:771521/1‑150 (MQ=255)
                                                                                                                               gATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATCAATAGCCAACTCAAAATTCCCACCTATTACCTTCCTCTGccc                       >  1:1869795/1‑148 (MQ=255)
                                                                                                                                 tAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACtt                     >  1:1559448/1‑150 (MQ=255)
                                                                                                                                 tAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTAAGACCTCCTTTCAAATGAATAGCCAACGCAAAATTCCCACCTATTCCCTTCCTCTGCACtt                     >  2:2012957/1‑150 (MQ=255)
                                                                                                                                  aaGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTa                    >  1:228576/1‑150 (MQ=255)
                                                                                                                                  aaGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTa                    <  1:1160597/150‑1 (MQ=255)
                                                                                                                                   aGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTac                   <  2:851404/150‑1 (MQ=255)
                                                                                                                                   aGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTac                   <  1:1507835/150‑1 (MQ=255)
                                                                                                                                      aCTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTAcaca                <  1:666975/150‑1 (MQ=255)
                                                                                                                                       ctctctCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACAt               <  2:1776336/150‑1 (MQ=255)
                                                                                                                                         ctctcACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTTCCTTCCTCTGCACTTACACATTc             >  1:2005391/1‑150 (MQ=255)
                                                                                                                                         ctctcACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTc             >  2:670104/1‑150 (MQ=255)
                                                                                                                                         ctctcACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTc             >  2:892202/1‑150 (MQ=255)
                                                                                                                                         ctctcACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTc             >  2:884493/1‑150 (MQ=255)
                                                                                                                                          tctcACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTATTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCg            <  1:108594/150‑1 (MQ=255)
                                                                                                                                           ctcACTCCAGCCAGAGCCCCCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTCCCTTCCTCTGCACTTACACATTCGt           >  2:362015/1‑150 (MQ=255)
                                                                                                                                           ctcACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGt           <  2:1670058/150‑1 (MQ=255)
                                                                                                                                           ctcACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTACAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGt           >  2:556633/1‑150 (MQ=255)
                                                                                                                                            tcACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGtt          <  2:1316422/150‑1 (MQ=255)
                                                                                                                                            tcACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGtt          <  2:733615/150‑1 (MQ=255)
                                                                                                                                             cACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAAt                                                              >  1:1060018/1‑97 (MQ=255)
                                                                                                                                             cACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAAt                                                              <  2:1060018/97‑1 (MQ=255)
                                                                                                                                             cACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACATCCTTTCATATGAATAGCCAACTCAAAAATAACACCTATGAGCTTGCTCTTCACTTACTAATTCGtca         >  2:108594/1‑148 (MQ=18)
                                                                                                                                              aCTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTaa        >  2:331311/1‑150 (MQ=255)
                                                                                                                                               cTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCa                                                                                 <  2:2001267/76‑1 (MQ=25)
                                                                                                                                               cTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAg       <  1:470795/150‑1 (MQ=255)
                                                                                                                                               cTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACACAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAg       >  2:440818/1‑150 (MQ=255)
                                                                                                                                                 ccagccagAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTc     >  2:1880782/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCCGCCAGCTGAAGAAATCGCTAATTCTTGCAATGTTAGCCACTGGCTAATAGTATTGAGCTGTTAGATAAGAACTCTCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGACCTCCTTTCAAATGAATAGCCAACTCAAAATTCACACCTATTACCTTCCTCTGCACTTACACATTCGTTAAGTCATA  >  NC_000913/3648170‑3648467

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