Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,947,051 A→C 6.2% intergenic (+19/+77) trpT → / ← hdfR tRNA‑Trp/flhDC operon transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,947,0510AC 6.2% 413.6 / 11.5 141intergenic (+19/+77)trpT/hdfRtRNA‑Trp/flhDC operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (61/71);  new base C (3/6);  total (64/77)
Fisher's exact test for biased strand distribution p-value = 5.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

Only 100 of 143 total aligned reads displayed.

GCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTATTTTGC  >  NC_000913/3946902‑3947197
                                                                                                                                                     |                                                                                                                                                   
gcctgTCACGCTGGGGGTCGGGGGTGCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGaaa                                                                                                                                                     <  1:158456/150‑1 (MQ=255)
 cctgTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAg                                                                                                                                                    <  1:1056990/150‑1 (MQ=255)
  ctgTCACGCAGGGGGTCGCGTGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTTGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGc                                                                                                                                                   <  1:1531404/150‑1 (MQ=255)
  ctgTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGc                                                                                                                                                   <  1:1351303/150‑1 (MQ=255)
   tgTCACGCTGGGGGTCGCGGGTTCGAGGCCCGTCCGTTGCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGATCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCTAGTCTCTCCGCCCCTGCCAGAAATCATCCTTTTCGAAAGCt                                                                                                                                                  <  2:662318/150‑1 (MQ=255)
          cAGGGGGTCGCGGGTGCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCTTTCGCTAAGGA‑tg                                                                                                                                          <  2:45697/150‑2 (MQ=255)
            gggggTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGGCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑tttt                                                                                                                                        <  1:998502/150‑1 (MQ=255)
             ggggGCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑ttttt                                                                                                                                       <  2:1251863/150‑1 (MQ=255)
              gggTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCATAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑tttttt                                                                                                                                      <  2:856341/150‑1 (MQ=255)
               ggTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑ttttttt                                                                                                                                     <  1:690357/150‑1 (MQ=255)
                gTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑tttttttt                                                                                                                                    <  1:1086893/150‑1 (MQ=255)
                 tCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTa                                                                                                                                   <  1:753340/150‑1 (MQ=255)
                 tCGCGGGTTCGAGGCCCGTCCGTTCCGCCCCCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCTAAACCGGGTGGTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTa                                                                                                                                   <  2:1724281/150‑1 (MQ=255)
                   gcgGGTTCGCGTCCCGTCCGTTCCGCCCCCCTAATTAGGGGCGTAGTTCAATTGGTCGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAACGCTAAGGA‑TTTTTTTTATc                                                                                                                                 <  1:1710150/150‑1 (MQ=255)
                   gcgGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATc                                                                                                                                 <  2:261385/150‑1 (MQ=255)
                      ggTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCACCAAAACCGGGAGTTGGGAGTTAGAGTCTCTCCGCCACACCCAGAAATCATCCTTTGCGAAAGCTAACGA‑GTTTTTTTATCTGa                                                                                                                              >  2:414739/1‑150 (MQ=255)
                       gTTCGAGTCCCGTCCGGTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGaa                                                                                                                             <  2:1176289/150‑1 (MQ=255)
                        ttCGAGTCACGTCCGTTCCGCCACCCTAATTAGGTGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTAAGCGAAAGCTAAGGA‑TTTTTTTTATCTGaaa                                                                                                                            <  1:561325/150‑1 (MQ=255)
                        gtCGAGTCGCGGCCGTTCCGCCACCCTAATTCGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGTTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGaaa                                                                                                                            <  1:251197/149‑1 (MQ=255)
                           gAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCAACGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATaa                                                                                                                         >  2:1009594/1‑150 (MQ=255)
                             gacccGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAAcc                                                                                                                       <  1:1228118/148‑1 (MQ=255)
                                  gtccgtTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCTTTCGCTAAGGa                                                                                                                                             <  2:551485/124‑1 (MQ=255)
                                  gtccgtTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCtctc                                                                                                                  <  1:1694970/150‑1 (MQ=255)
                                     cgtTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGa                                                                                                               <  1:1385701/150‑1 (MQ=255)
                                       ttCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCTAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑‑TTTTTTTATCTGAAATAACCCTCTCCGAAGt                                                                                                            <  1:197451/150‑1 (MQ=255)
                                       ttCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAg                                                                                                             >  2:1346922/1‑150 (MQ=255)
                                        tCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGt                                                                                                            <  1:1343330/150‑1 (MQ=255)
                                         ccGCCACCCTAATTAGGTGCGTAGTTCAATTGGTAGAGCACCGGTCTCGAATTCCGGGTGTTGGGAGTTGGAGTCTCTACGCCCAGGCCTGAAATCATCCTTAGCGAAAGCTAAGGA‑TGGTTGTTATCTGACATACCCCTCTCCTCAGTa                                                                                                           >  1:132028/1‑150 (MQ=255)
                                         ccGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAACCCGGGTGTTGGGAGGTCGAGTCTCTCCGCCCCTGCCAGAAATCCTCCTTAGCGAAAGCTAAGGA‑GTTTTTGTATCTGAAATAACCCTCTCAGAAGTa                                                                                                           >  1:126569/1‑150 (MQ=255)
                                         ccGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTa                                                                                                           >  1:625525/1‑150 (MQ=255)
                                         ccGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATACTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTa                                                                                                           >  2:1024913/1‑150 (MQ=255)
                                           gCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTaaa                                                                                                         >  1:628207/1‑150 (MQ=255)
                                           gCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTaaa                                                                                                         >  1:41881/1‑150 (MQ=255)
                                            ccACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAAt                                                                                                        <  2:926477/150‑1 (MQ=255)
                                              aCCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATcc                                                                                                      <  1:549806/150‑1 (MQ=255)
                                                    aTTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTAc                                                                                                <  1:951875/150‑1 (MQ=255)
                                                       aGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCgg                                                                                             <  1:32164/150‑1 (MQ=255)
                                                              aaGTACAATTGGTAGAGCACCGGTCTCCAAAACCGCGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAAACATCCTTAGCGAAAGCTAACGA‑TTTTATTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCtt                                                                                      >  1:1778759/2‑150 (MQ=255)
                                                                 ttCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATACCCCTCTCCGAAGTAAATCCTTCTCCCGGCATCCTTgcc                                                                                   >  1:319255/1‑150 (MQ=255)
                                                                  tCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTgcca                                                                                  >  2:1324130/1‑150 (MQ=255)
                                                                  tCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTgcca                                                                                  <  1:823506/150‑1 (MQ=255)
                                                                   cAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTgccag                                                                                 >  1:119508/1‑150 (MQ=255)
                                                                       tGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCAt                                                                             >  1:1450875/1‑150 (MQ=255)
                                                                         gTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTc                                                                           >  1:1047108/1‑150 (MQ=255)
                                                                         gTAGAGCACCGGGATCCAAAACCGGGTGTTGGGAGTTCGTGGCTCTCCGACCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTc                                                                           <  2:256335/150‑1 (MQ=255)
                                                                          tAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑GTTTTTTTATCTGAAATAACCCTCGCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCa                                                                          >  1:116158/1‑150 (MQ=255)
                                                                          tAGAGAACCGGTCTCCAAAACCGGGTGGTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTACGGA‑TTTTTTTTATCGGAAATAACCCTCTCCGAAGTAAATCATTCTACCGTCATCCTTGCAAGCCATTcg                                                                          >  1:116162/1‑149 (MQ=255)
                                                                             atcACCGTTCTCCAAATCAGTGTGTTGGGAGTTCGATTCTCTCCGACCCTGCCAGAAATCATCCTTAGCGAAAGCTTATGA‑TTTTTTTTATCTGAAATAACCCTCTGCGAAGTTAATCCTTCTACCGGCATCCTTGCCAGCCATTCatat                                                                       <  1:229892/148‑1 (MQ=255)
                                                                              tcACCGGTATCCAAAAACGGGTGGTGGGATTTCGAGTCTCTCCGCCCCTTCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATAtt                                                                      <  2:1387256/149‑1 (MQ=255)
                                                                               cACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTATCTCAGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTa                                                                     <  1:234417/150‑1 (MQ=255)
                                                                                 cgaactcccAACACCCGGTTTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGATAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAAt                                                                   >  2:706950/8‑150 (MQ=255)
                                                                                  cGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATa                                                                  <  2:16661/150‑1 (MQ=255)
                                                                                  cGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATa                                                                  <  2:435271/150‑1 (MQ=255)
                                                                                  cGGTCTCCAAAACCGGGTGTTGGGAGCTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTACCAGCCATTCATATTAATa                                                                  <  2:1597199/150‑1 (MQ=255)
                                                                                        ccAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACtta                                                            >  1:1201152/1‑149 (MQ=255)
                                                                                        ccAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTATTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTc                                                            >  2:1227965/1‑150 (MQ=255)
                                                                                        ccAAAACCGGGTGTTGGGAGTTCGAGTCTATCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATATGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTc                                                            >  1:469836/1‑150 (MQ=255)
                                                                                        ccAAAACCGGGTGTTGGGAGTTCGAGTCGCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATACCCCTCTCCGAAGTAAATCCTTCTACCGGAAGCCGTGCCAGCCATTCATATTAAGACACGTc                                                            >  2:1161986/1‑150 (MQ=255)
                                                                                         cAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCa                                                           >  1:193419/1‑150 (MQ=255)
                                                                                          aaatCCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAATCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCAt                                                          <  2:1160477/150‑1 (MQ=255)
                                                                                          aaaaCCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCAt                                                          >  1:1768441/1‑150 (MQ=255)
                                                                                            aaCCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATcc                                                        <  1:673562/150‑1 (MQ=255)
                                                                                               cGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCGTATTAATACACTTCATCCAGc                                                     >  2:1776660/1‑150 (MQ=255)
                                                                                                  gtgtTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACg                                                  <  1:1665480/150‑1 (MQ=255)
                                                                                                   ggtTGGGAGTTCGAGTCTCTCCGCCCCTTCCAGAAATCATCCTTAGCGAAAGCTAAGTA‑TTTTTTTTATCTGAAATAACCCTCTCAGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGt                                                 <  1:73667/149‑1 (MQ=255)
                                                                                                      tGGGAGTTCGCGTCTCTCCGCCCCTGACAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTAACTGAAATAACCCTCTCCGAACTAACTCCTTCTACCTGCATCCTTACCAGCAAGTCATATTAATACACATCATCCAGCACGTTaa                                              >  1:894779/1‑150 (MQ=255)
                                                                                                       gggAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAAt                                             >  1:1453938/1‑150 (MQ=255)
                                                                                                           gTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCTTTAATATTAATACACTTCATCCAGAACGTTAATTTTc                                         <  2:350633/150‑1 (MQ=255)
                                                                                                           gTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTc                                         <  2:779135/150‑1 (MQ=255)
                                                                                                            ttCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTATTTTAACTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTAGCCAGCCATACATATTAATACACTTCATCCAGCACGTTAATTGTCa                                        >  1:872079/1‑150 (MQ=255)
                                                                                                              cGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCaaa                                      >  1:11731/1‑150 (MQ=255)
                                                                                                              cGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCaaa                                      >  1:982176/1‑150 (MQ=255)
                                                                                                              cGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑‑TTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCaaaa                                     >  2:424429/1‑150 (MQ=255)
                                                                                                               gAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCaaaa                                     >  1:1654630/1‑150 (MQ=255)
                                                                                                               gAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCaaaa                                     >  2:1807732/1‑150 (MQ=255)
                                                                                                                aGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAg                                    >  2:1048839/1‑150 (MQ=255)
                                                                                                                  tctctcCGCCCCTGCCAGATATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTTCCGGCTTCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGAt                                  <  2:1549510/150‑1 (MQ=255)
                                                                                                                  tctctcCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTTCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGAt                                  <  2:1542446/150‑1 (MQ=255)
                                                                                                                   ctctcCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATc                                 >  2:379628/1‑150 (MQ=255)
                                                                                                                     ctcCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATcgc                               >  2:1772255/1‑150 (MQ=255)
                                                                                                                      tcCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATcgcg                              <  1:339235/150‑1 (MQ=255)
                                                                                                                      tcCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATcgcg                              <  1:1703119/150‑1 (MQ=255)
                                                                                                                       ccGCCCCTGCCAGAAATCATCCTTAGCGAACGCTAATGA‑TTTTTTTTATCTGAAATAACACTCTCAGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGa                             <  2:1253735/150‑1 (MQ=255)
                                                                                                                       ccGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGa                             >  2:1569232/1‑150 (MQ=255)
                                                                                                                              tGCAAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATACATTTTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACg                      <  2:321655/150‑1 (MQ=255)
                                                                                                                               gCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTATATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCAAGAACGTTAATTTTCAAAAGATCGCGAATCAACGc                     >  1:437906/1‑150 (MQ=255)
                                                                                                                                ccAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCa                    <  1:549601/150‑1 (MQ=255)
                                                                                                                                acAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTATTTTTCAAAAGATCGCGAATCAACGCa                    <  1:60475/149‑1 (MQ=255)
                                                                                                                                   aaaaaaatcCTTAGCTTTCGCTAAGGATTTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGATTCTACGCAtt                  >  1:379666/7‑150 (MQ=255)
                                                                                                                                   aaaaaaatcCTTAGCTTTCGCTAAGGA‑‑TTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCAtttt                >  2:1549670/7‑150 (MQ=255)
                                                                                                                                    aaaTCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCAtttt                <  2:587410/150‑1 (MQ=255)
                                                                                                                                    aaaTCATCCTTAGCGAAAGCGAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTAATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATAAACGCAtttt                >  2:557971/1‑150 (MQ=255)
                                                                                                                                        catcCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATc            <  2:292356/150‑1 (MQ=255)
                                                                                                                                         atcCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCg           <  2:1047108/150‑1 (MQ=255)
                                                                                                                                            cTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTa        <  2:502379/150‑1 (MQ=255)
                                                                                                                                             ttAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCt         <  2:1200965/148‑1 (MQ=255)
                                                                                                                                                 cGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTTCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTATTTTg   <  1:1324130/150‑1 (MQ=255)
                                                                                                                                                 cGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTATTTTg   <  2:1450875/150‑1 (MQ=255)
                                                                                                                                                  gAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTATTTTGc  <  2:628207/150‑1 (MQ=255)
                                                                                                                                                  gAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTATTTTGc  >  1:579394/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCCTGTCACGCAGGGGGTCGCGGGTTCGAGTCCCGTCCGTTCCGCCACCCTAATTAGGGGCGTAGTTCAATTGGTAGAGCACCGGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGA‑TTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCATTCATATTAATACACTTCATCCAGCACGTTAATTTTCAAAAGATCGCGAATCAACGCATTTTTATCGCTATTTTGC  >  NC_000913/3946902‑3947197

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