Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,213,182 A→T 5.1% intergenic (+23/‑52) rrfE → / → yjaA 5S ribosomal RNA of rrnE operon/stress‑induced protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,213,1820AT 5.1% 380.7 / 4.3 127intergenic (+23/‑52)rrfE/yjaA5S ribosomal RNA of rrnE operon/stress‑induced protein
Reads supporting (aligned to +/- strand):  ref base A (63/57);  new base T (5/2);  total (68/59)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 127 total aligned reads displayed.

TTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACC  >  NC_000913/4213039‑4213330
                                                                                                                                                |                                                                                                                                                     
ttGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑c                                                                                                                                                >  2:337561/1‑150 (MQ=255)
 tGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑cc                                                                                                                                               <  2:233810/150‑1 (MQ=255)
   ccTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCgg                                                                                                                                             >  1:216663/1‑150 (MQ=255)
   ccTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCgg                                                                                                                                             <  2:60401/150‑1 (MQ=255)
     tGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCgggg                                                                                                                                           <  1:1647944/150‑1 (MQ=255)
         ggcAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGtttt                                                                                                                                       >  2:1552518/1‑150 (MQ=255)
         ggcAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGtttt                                                                                                                                       <  2:664876/150‑1 (MQ=255)
           cAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAAGGCCGTAGCGCCGATGGTAGTGGGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGG‑CCTTATGGACCGGGGtttttc                                                                                                                                     <  2:1377067/150‑2 (MQ=14)
           cAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAATTTA‑GAAAAACCGCGGTCCTTAAGG‑CCGGGGtatttt                                                                                                                                     >  1:90346/1‑150 (MQ=255)
            aGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTg                                                                                                                                    >  2:1692873/1‑150 (MQ=255)
             gTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACACAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGTGAACTGCCAGGCATCAAATTA‑GTAAATCCCCTGTCCATAAGG‑CCGGGGTTTTTTGc                                                                                                                                   <  1:253283/150‑1 (MQ=17)
               tgcgcGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCACGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCat                                                                                                                                 <  2:1603117/149‑1 (MQ=255)
                       ggtCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGTAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAatta                                                                                                                         >  2:943012/1‑150 (MQ=255)
                        gtCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCAAAAGG‑CCGGGGGTTTTTGCATATCAatttt                                                                                                                        >  1:72202/1‑148 (MQ=255)
                           ccACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGGTGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTg                                                                                                                     <  2:107572/150‑1 (MQ=255)
                            cACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGc                                                                                                                    >  1:463351/1‑150 (MQ=255)
                              ccTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCat                                                                                                                  <  2:527731/150‑1 (MQ=255)
                                tGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCatga                                                                                                                >  1:594845/1‑150 (MQ=255)
                                  acaccATCCCGAACTCAGAAGTGAAACGCCGTAGCTCCGATGGTAGTGTGGTGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCatgatg                                                                                                              <  2:262719/147‑1 (MQ=255)
                                  aCCCCATGCCGAACTCAGAAGTGAAACGCCTTAGCGCCGATGGTAGTGTGGTGTATCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCatgatg                                                                                                              <  2:912711/150‑1 (MQ=255)
                                   ccccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCAtgatga                                                                                                             >  1:417935/1‑150 (MQ=255)
                                             aaCTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAAtctc                                                                                                   >  1:756861/1‑150 (MQ=255)
                                                tCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATg                                                                                                <  2:916340/150‑1 (MQ=255)
                                                tCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATg                                                                                                >  2:658109/1‑150 (MQ=255)
                                                 cAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGt                                                                                               >  1:411541/1‑150 (MQ=255)
                                                    aaGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAg                                                                                            >  1:1611415/1‑150 (MQ=255)
                                                     aGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGt                                                                                           <  1:639491/150‑1 (MQ=255)
                                                        gAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCt                                                                                        >  2:1578125/1‑150 (MQ=255)
                                                         aaaCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTg                                                                                       >  1:1490452/1‑150 (MQ=255)
                                                         aaaCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTg                                                                                       <  2:370369/150‑1 (MQ=255)
                                                         aaaCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTg                                                                                       <  1:1385657/150‑1 (MQ=255)
                                                           aCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGta                                                                                     <  1:665358/150‑1 (MQ=255)
                                                           aCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGta                                                                                     >  1:379582/1‑150 (MQ=255)
                                                             gCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGtata                                                                                   <  2:41189/150‑1 (MQ=255)
                                                               cGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATc                                                                                 <  2:390688/150‑1 (MQ=255)
                                                                    cgcCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAAt                                                                            <  1:946674/150‑1 (MQ=255)
                                                                       cGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTTTTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATtcg                                                                         >  2:224431/1‑150 (MQ=255)
                                                                       cGATGGTAGTGTGGGGTCTCCCCATGAGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATtcg                                                                         >  2:1290992/1‑150 (MQ=255)
                                                                          tGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATtcgtcg                                                                      <  1:605235/150‑1 (MQ=255)
                                                                             tAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTa                                                                    >  1:1641088/1‑150 (MQ=255)
                                                                                ggtgGGTTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCAAGTCAGTTCTGTATATCCAAATTCGTCGTAATCa                                                                <  2:65970/149‑1 (MQ=255)
                                                                                   ggggTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATATCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAAt                                                             >  1:858428/1‑150 (MQ=255)
                                                                                    gggTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAAtt                                                            <  1:497226/150‑1 (MQ=255)
                                                                                      gTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTAc                                                          <  2:1360939/150‑1 (MQ=255)
                                                                                       tctcCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACt                                                         <  2:132634/150‑1 (MQ=255)
                                                                                         tcCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGt                                                       >  2:1536464/1‑150 (MQ=255)
                                                                                          ccccATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCTTGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGtt                                                      >  1:1078857/1‑150 (MQ=255)
                                                                                          ccccATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGtt                                                      >  2:927996/1‑150 (MQ=255)
                                                                                          ccccATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGtt                                                      >  1:1611041/1‑150 (MQ=255)
                                                                                            ccATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTcg                                                    <  2:463351/150‑1 (MQ=255)
                                                                                            ccATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTcg                                                    >  2:971833/1‑150 (MQ=255)
                                                                                            ccATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTcg                                                    >  1:106061/1‑150 (MQ=255)
                                                                                            ccATGCGAGAGTAGGAAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTcg                                                    >  2:461173/1‑150 (MQ=255)
                                                                                              aTGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTcgcg                                                  <  2:143473/150‑1 (MQ=255)
                                                                                               tGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCACGGTCATTAAGG‑ACGGGGTTTTTTGAATATCAATTATTTGCATGATGAAGGGACTCTCATGTCAGTTCTGTATATCAAAATTAGTCTTAATCAAATTACTGTTCGCGa                                                 >  1:845057/1‑150 (MQ=255)
                                                                                                gCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGAt                                                >  2:245536/1‑150 (MQ=255)
                                                                                                 cGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATc                                               <  1:268760/150‑1 (MQ=255)
                                                                                                  gagaGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCt                                              <  1:15134/150‑1 (MQ=255)
                                                                                                  gagaGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCt                                              <  1:205127/150‑1 (MQ=255)
                                                                                                    tagTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTg                                            <  2:1640428/149‑1 (MQ=255)
                                                                                                     agTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGa                                           <  1:1533315/150‑1 (MQ=255)
                                                                                                      gTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGaa                                          <  2:630002/150‑1 (MQ=255)
                                                                                                      gTAGGGAACTGCAAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATAATATGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTAGaa                                          >  1:48578/1‑150 (MQ=255)
                                                                                                       tAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGaaa                                         >  1:387376/1‑150 (MQ=255)
                                                                                                        aGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAg                                        >  2:581752/1‑150 (MQ=255)
                                                                                                        aGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAg                                        <  2:718558/150‑1 (MQ=255)
                                                                                                        aGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAg                                        <  2:605514/150‑1 (MQ=255)
                                                                                                         gggAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGc                                       >  2:856065/1‑150 (MQ=255)
                                                                                                         gggAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCAGCTCAGTTCTGTATATCCAATTTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGc                                       >  1:235172/1‑150 (MQ=255)
                                                                                                           gAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCaa                                     >  1:1617583/1‑150 (MQ=255)
                                                                                                             aCTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAAc                                   >  1:555542/1‑150 (MQ=255)
                                                                                                               tGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGt                                 <  2:498781/150‑1 (MQ=255)
                                                                                                                gCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTg                                <  1:1543676/150‑1 (MQ=255)
                                                                                                                    ggCATCAAATTAGGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAg                             <  2:1641088/150‑1 (MQ=255)
                                                                                                                    ggCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGGAATCAAATTACTGTTCGCGAGCTTGAAAGCAAACGTGAAgt                            >  1:1694696/1‑150 (MQ=255)
                                                                                                                       aTCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTc                         <  1:1578810/150‑1 (MQ=255)
                                                                                                                         cAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGATGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAg                       >  1:1444501/1‑150 (MQ=255)
                                                                                                                         cAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTAATGTTCGCTATCTTGAAAGCAAACGTGAAGTGTCAg                       >  2:120671/1‑150 (MQ=255)
                                                                                                                         aaaaTTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTTCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAg                       <  1:975271/149‑1 (MQ=255)
                                                                                                                         aaaaTTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAg                       <  2:1244515/149‑1 (MQ=255)
                                                                                                                          aaaTTA‑GAAAAACCCCGGTCCATAAGG‑CCTGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAgg                      <  1:943012/150‑1 (MQ=255)
                                                                                                                              tA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGAc                  >  2:1102595/1‑150 (MQ=255)
                                                                                                                               a‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACTTGAAGTGTCAGGAGACg                 >  2:85370/1‑150 (MQ=255)
                                                                                                                               a‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACg                 >  2:566032/1‑150 (MQ=255)
                                                                                                                                caaaaaaCCCCGG‑CCTTATGGACCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACg                 >  1:809221/3‑150 (MQ=255)
                                                                                                                                 gAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGc                >  1:794433/1‑150 (MQ=255)
                                                                                                                                 gAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGc                >  2:1359698/1‑150 (MQ=255)
                                                                                                                                 aaaaaaCCCCGG‑CCTTATGGACCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGc                >  2:1038266/2‑150 (MQ=255)
                                                                                                                                   aaaaCCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTg              >  1:1662824/1‑150 (MQ=255)
                                                                                                                                    aaaCCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGc             <  1:832766/150‑1 (MQ=255)
                                                                                                                                    aaaCCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGc             <  1:1003123/150‑1 (MQ=255)
                                                                                                                                    aaaCCCCGG‑CCTTATGGACCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGc             >  2:67706/1‑150 (MQ=255)
                                                                                                                                       ccccGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCtt          >  1:1071990/1‑150 (MQ=255)
                                                                                                                                        cccGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTc         >  2:684516/1‑150 (MQ=255)
                                                                                                                                        cccGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTc         <  1:1386265/150‑1 (MQ=255)
                                                                                                                                        cccGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTc         >  2:1479519/1‑150 (MQ=255)
                                                                                                                                          cGGTCCATCAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAg       <  2:53592/150‑1 (MQ=255)
                                                                                                                                          cGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAg       <  1:1692873/150‑1 (MQ=255)
                                                                                                                                           ggTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATAACCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGt      >  1:1184989/1‑150 (MQ=255)
                                                                                                                                               cATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAAcc  >  2:711318/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
TTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTA‑GAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAACC  >  NC_000913/4213039‑4213330

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