Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,398,552 G→A 5.4% A24V (GCT→GTT)  uspE ← stress‑induced protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,398,5520GA 5.4% 235.1 / 2.2 74A24V (GCT→GTT) uspEstress‑induced protein
Reads supporting (aligned to +/- strand):  ref base G (37/33);  new base A (2/2);  total (39/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.49e-01

GATCCAGGCTGTACGCTGGCTGATGACGCCCTGACGCATAGCGGTACGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCTGGTGAATGAAACAGTGA  >  NC_000913/1398407‑1398699
                                                                                                                                                 |                                                                                                                                                   
gATCCAGGCTGTACGCTGGCTGATGACGCCCTGACGCATAGCGGTACGTTCGTCCGGGGAGAGCAGGGTTGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTc                                                                                                                                                  >  6:64715/1‑149 (MQ=255)
 aTCCAGGCTGTACGCTGGCTGATGACGCCCTGACGCATAGCGGTACGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTcg                                                                                                                                                 >  3:85364/1‑149 (MQ=255)
    cAGGCTGTACGCTGGCTGATGACGCCCTGACGCATAGCGGTACGTTCGTCCGGGGAGAGCAGTGTGGTCATTTCTTAAGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTcgccg                                                                                                                                              <  1:337908/149‑1 (MQ=255)
     aGGCTGTACGCTGGCTGATGACGCCCTGACGCATAGCGGTACGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTcgccg                                                                                                                                              >  3:31853/1‑148 (MQ=255)
     aGGCTGTACGCTGGCTGATGACGCCCTGACGCATAGCGGTACGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTcgccg                                                                                                                                              >  7:143325/1‑148 (MQ=255)
      ggCTGTACGCTGGCTGATGCCGCCCTGACGCATAGCGGTACGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTcgccgc                                                                                                                                             >  3:262090/1‑148 (MQ=255)
      ggCTGTACGCTGGCTGATGACGCCCTGACGCATAGCGGTACGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTcgccgc                                                                                                                                             >  2:314956/1‑148 (MQ=255)
             cGCTGGCTGATGACGCCCTGACGCATAGCGGTACGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCtgg                                                                                                                                     <  2:219018/149‑1 (MQ=255)
                tggctgATGACGCCCTGACGCATAGCGGTACGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCtggttg                                                                                                                                  >  5:339189/1‑149 (MQ=255)
                  gctgATGACGCCCTGACGCATAGCGGTACGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCtggttg                                                                                                                                  <  6:43579/147‑1 (MQ=255)
                         aCGCCCTGACGCATAGCGGTACGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTg                                                                                                                         <  1:282643/149‑1 (MQ=255)
                            cccTGACGCATAGCGGTACGTTCGTCCTGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTTCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGtt                                                                                                                      <  6:161527/149‑1 (MQ=255)
                              cTGACGCATAGCGGTACGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAg                                                                                                                    <  4:31853/149‑1 (MQ=255)
                                     aTAGCGGTACGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGa                                                                                                              >  7:145222/1‑148 (MQ=255)
                                     aTAGCGGTACGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGAt                                                                                                             <  1:314956/149‑1 (MQ=255)
                                          ggTACGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATaacaa                                                                                                        >  8:175920/1‑149 (MQ=255)
                                           gTACGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATaacaac                                                                                                       <  6:72787/149‑1 (MQ=255)
                                               gTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACg                                                                                                      <  6:235321/146‑1 (MQ=255)
                                               gTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACg                                                                                                      >  5:235321/1‑146 (MQ=255)
                                               gTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAg                                                                                                    >  2:337183/1‑148 (MQ=255)
                                                   gTCCGGGTAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGAAAATAAACAGCTCGCCGCAATGCTCGATGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGt                                                                                               >  2:222895/1‑149 (MQ=255)
                                                     ccGGGGAGAGCAGGGTGGTCATTTCGTATTAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATTTAAATAAACAGCTCGCCGCAATGCTGGTGGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTc                                                                                             <  1:240339/149‑1 (MQ=255)
                                                     ccGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAACTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTc                                                                                             >  5:41265/1‑149 (MQ=255)
                                                      cGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTc                                                                                             <  4:262090/148‑1 (MQ=255)
                                                      cGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCt                                                                                            <  8:143325/149‑1 (MQ=255)
                                                        gggAGAGCAGGGTGGTCATTTCGTATTATAAGTCATAGATAGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATTATGGTTGGTCTTCCTGGTTAGGATAGATAACAACGAGAATGTTCTGa                                                                                          >  3:261020/1‑149 (MQ=255)
                                                          gagaGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATa                                                                                        <  4:85364/149‑1 (MQ=255)
                                                            gagCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACa                                                                                      >  8:348683/1‑149 (MQ=255)
                                                              gCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATa                                                                                    <  6:301744/149‑1 (MQ=255)
                                                               cAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGc                                                                                                   >  5:627/1‑133 (MQ=255)
                                                               cAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGc                                                                                                   <  6:627/133‑1 (MQ=255)
                                                               cAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAg                                                                                   <  3:266568/149‑1 (MQ=255)
                                                               cAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAACTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAg                                                                                   <  6:41265/149‑1 (MQ=255)
                                                                    tggtCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATa                                                                              <  7:300558/149‑1 (MQ=255)
                                                                       tCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAg                                                                           >  1:162024/1‑149 (MQ=255)
                                                                        cATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAgg                                                                          <  6:339189/149‑1 (MQ=255)
                                                                         aTTTCGTATGAGAAGTCATAGATCGGCAAAAATGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAgg                                                                          <  4:261020/148‑1 (MQ=255)
                                                                         aTTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAggg                                                                         >  8:361491/1‑149 (MQ=255)
                                                                          tttCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGt                                                                        >  7:131376/1‑149 (MQ=255)
                                                                          tttCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTAGACGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGt                                                                        >  5:176661/1‑149 (MQ=255)
                                                                               tATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCt                                                                   <  7:175920/149‑1 (MQ=255)
                                                                               tATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCt                                                                   >  4:180619/1‑149 (MQ=255)
                                                                                 tGAGAAGTCATAGATCGGCAAAAAGTCCTTAATTTTGCCACCAATCCGTTTATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGCTAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTa                                                                 <  1:222895/149‑1 (MQ=255)
                                                                                   agaAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTac                                                                >  1:370280/1‑148 (MQ=255)
                                                                                    gaAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGtt                                                                                              <  5:239922/117‑1 (MQ=255)
                                                                                    gaAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGtt                                                                                              >  6:239922/1‑117 (MQ=255)
                                                                                        tCATAGATCGGCAAAAAGGCTTTAATTTTGCCCCCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACaacaac                                                          >  1:106653/1‑149 (MQ=255)
                                                                                         cATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACaacaac                                                          <  1:337183/148‑1 (MQ=255)
                                                                                          aTAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTg                                                        >  4:321525/1‑149 (MQ=255)
                                                                                              aTCGGCAAAAAGGCTTTAATTTTGCCCCCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCaa                                                    >  6:247213/1‑149 (MQ=255)
                                                                                                     aaaaGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAACTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGtt                                             >  4:255869/1‑149 (MQ=255)
                                                                                                     aaaaGGCTTTAATTTTGCCACCAAACCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATACCAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGtt                                             >  6:102416/1‑149 (MQ=255)
                                                                                                      aaaGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGttt                                            >  8:300593/1‑149 (MQ=255)
                                                                                                      aaaGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAACTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGttt                                            <  3:255869/149‑1 (MQ=255)
                                                                                                         ggCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTa                                         <  3:105392/149‑1 (MQ=255)
                                                                                                            tttAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAAt                                       >  2:295197/1‑148 (MQ=255)
                                                                                                              tAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGATGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACACCTGTCAACGCAGTATGTAATTa                                     >  4:47907/1‑148 (MQ=255)
                                                                                                               aaTTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTaa                                    >  8:348264/1‑148 (MQ=255)
                                                                                                                aTTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTaaa                                   <  7:161173/148‑1 (MQ=255)
                                                                                                                        accaATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAAcc                          >  3:34189/1‑149 (MQ=255)
                                                                                                                         ccgaTCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTa                                         <  3:9868/130‑1 (MQ=255)
                                                                                                                         ccgaTCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTa                                         >  4:9868/4‑133 (MQ=255)
                                                                                                                         ccaATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAAcc                          >  8:85014/1‑148 (MQ=255)
                                                                                                                         ccaATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAAcc                          >  1:154404/1‑148 (MQ=255)
                                                                                                                                gTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCTg                  <  4:34189/149‑1 (MQ=255)
                                                                                                                                gTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCTg                  >  3:281892/1‑149 (MQ=255)
                                                                                                                                gTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCTg                  <  2:154404/149‑1 (MQ=255)
                                                                                                                                gTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTACTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCTg                  <  8:145222/149‑1 (MQ=255)
                                                                                                                                 ttGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCTgg                 <  1:295197/149‑1 (MQ=255)
                                                                                                                                        aaataaaCAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCTGGtgaatg           >  8:283718/1‑148 (MQ=255)
                                                                                                                                        aaataaaCAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCAGGtgaatg           >  1:134280/1‑148 (MQ=255)
                                                                                                                                          ataaaCAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCTGGTGAATGaaa        <  2:162024/149‑1 (MQ=255)
                                                                                                                                               cAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCTGGTGAATGAAACAGtg   <  7:300593/149‑1 (MQ=255)
                                                                                                                                                aGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCTGGTGAATGAAACAGtga  <  8:340145/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GATCCAGGCTGTACGCTGGCTGATGACGCCCTGACGCATAGCGGTACGTTCGTCCGGGGAGAGCAGGGTGGTCATTTCGTATGAGAAGTCATAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCCTTACAACAACTGTCAACGCAGTTTGTAATTAAAAGATTAACCCATATCTGGTGAATGAAACAGTGA  >  NC_000913/1398407‑1398699

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