Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,210,447:1 +G 5.3% coding (97/1014 nt) tsaD ← N(6)‑L‑threonylcarbamoyladenine synthase, TsaD subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,210,4471.G 5.3% 398.4 / 15.5 95coding (97/1014 nt)tsaDN(6)‑L‑threonylcarbamoyladenine synthase, TsaD subunit
Reads supporting (aligned to +/- strand):  ref base . (42/48);  new base G (3/2);  total (45/50)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCGGTATAGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTTTATTACCGCGCAGTGTAGGACCAATGCGGGTTGA  >  NC_000913/3210303‑3210596
                                                                                                                                                 |                                                                                                                                                     
cTGCGGTATAGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGAc                                                                                                                                                  <  1:904081/150‑1 (MQ=255)
 tGCGGTATAGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACt                                                                                                                                                 <  1:1114035/150‑1 (MQ=255)
 tGCGGTATAGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACt                                                                                                                                                 <  1:525477/150‑1 (MQ=255)
 tGCGGTATAGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGGTGAAATTGCACGCTGACt                                                                                                                                                 <  2:401166/150‑1 (MQ=255)
   cGGTATAGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACtat                                                                                                                                               >  1:187260/1‑150 (MQ=255)
    ggTATAGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCAc                                                                                                                                                        <  2:490523/141‑1 (MQ=255)
    ggTATAGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCAc                                                                                                                                                        >  1:490523/1‑141 (MQ=255)
      tataGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACa                                                                                                                                            >  1:1330004/1‑150 (MQ=255)
       ataGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGGTGAAATTGCACGCTGACtacggcg                                                                                                                                           <  2:1276480/150‑6 (MQ=255)
           gCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAAttgg                                                                                                                                       >  2:105542/1‑150 (MQ=255)
                 gCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCt                                                                                                                                 <  2:916347/150‑1 (MQ=255)
                 gCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCt                                                                                                                                 >  1:916347/1‑150 (MQ=255)
                  cATCAATGTCTTTTGGCGTTAAACCAGACTACTTTAGCGCCGCCTGGATCAACGGTAGTGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTATTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTa                                                                                                                                <  2:1364076/150‑1 (MQ=255)
                    tCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAAc                                                                                                                              <  1:237406/150‑1 (MQ=255)
                     cAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACa                                                                                                                             <  2:60357/150‑1 (MQ=255)
                     cAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACa                                                                                                                             >  1:177566/1‑150 (MQ=255)
                       aTGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACaaa                                                                                                                           >  2:732701/1‑150 (MQ=255)
                       aTGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACaaa                                                                                                                           <  1:333876/150‑1 (MQ=255)
                        tGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAAc                                                                                                                          <  1:183820/150‑1 (MQ=255)
                           cTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGCGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCtt                                                                                                                       >  2:1536076/1‑150 (MQ=255)
                            ttttGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCttt                                                                                                                      >  1:273802/1‑150 (MQ=255)
                                 ccGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCAt                                                                                                                 <  2:763899/150‑1 (MQ=255)
                                  cGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCAtc                                                                                                                >  2:353227/1‑150 (MQ=255)
                                      aaaCCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCAtcgtcg                                                                                                            >  2:1393246/1‑150 (MQ=255)
                                      aaaCCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCAtcgtcg                                                                                                            <  1:1393246/150‑1 (MQ=255)
                                      aaaCCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCAtcgtcg                                                                                                            <  2:694376/150‑1 (MQ=255)
                                      aaaCCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCAtcgtcg                                                                                                            >  1:1126944/1‑150 (MQ=255)
                                       aaCCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATcgtcgt                                                                                                           >  1:1014854/1‑150 (MQ=255)
                                        aCCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTa                                                                                                          <  2:886699/150‑1 (MQ=255)
                                         ccAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTaa                                                                                                         >  1:894459/1‑150 (MQ=255)
                                           aGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAAt                                                                                                       <  1:353227/150‑1 (MQ=255)
                                             aCTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATgg                                                                                                     <  2:444707/150‑1 (MQ=255)
                                              cTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGc                                                                                                    <  2:169300/150‑1 (MQ=255)
                                              cTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGc                                                                                                    <  1:1050030/150‑1 (MQ=255)
                                                ccTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGa                                                                                                  <  2:1014854/150‑1 (MQ=255)
                                                ccTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGa                                                                                                  >  1:960309/1‑150 (MQ=255)
                                                 cTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGAt                                                                                                 <  2:564341/150‑1 (MQ=255)
                                                 cTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGAt                                                                                                 >  1:1200017/1‑150 (MQ=255)
                                                 cTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGACCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGAt                                                                                                 >  2:393276/1‑150 (MQ=255)
                                                       cgccgcCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTCACTATACAATTGGTTGGCTAACAACCGTTTTTCATCGTCGTAAATGGCGATGCCGGt                                                                                           >  1:975359/1‑150 (MQ=255)
                                                         ccgccTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGttt                                                                                         <  2:1516183/150‑1 (MQ=255)
                                                              tGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCg                                                                                    <  1:98234/150‑1 (MQ=255)
                                                               ggATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATgg                                                                                                     >  2:470101/1‑132 (MQ=255)
                                                               ggATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATgg                                                                                                     <  1:470101/132‑1 (MQ=255)
                                                                 aTCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAg                                                                                 <  2:1505221/150‑1 (MQ=255)
                                                                 aTCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAg                                                                                 >  1:1333491/1‑150 (MQ=255)
                                                                 aTCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAg                                                                                 <  1:20250/150‑1 (MQ=255)
                                                                   cAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGa                                                                               >  1:43459/1‑150 (MQ=255)
                                                                   cAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGa                                                                               >  1:1083491/1‑150 (MQ=255)
                                                                     acggtacggtTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAg                                                                             >  2:976480/1‑150 (MQ=255)
                                                                        gtacggtTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATg                                                                                                      <  1:618986/122‑1 (MQ=255)
                                                                        gtacggtTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATg                                                                                                      >  2:618986/1‑122 (MQ=255)
                                                                            ggtTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAAt                                                                      <  2:1126944/150‑1 (MQ=255)
                                                                               tttACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATAcc                                                                   <  2:187904/150‑1 (MQ=255)
                                                                                ttACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATAccc                                                                  <  2:380632/150‑1 (MQ=255)
                                                                                ttACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATAccc                                                                  <  2:374987/150‑1 (MQ=255)
                                                                                  aCGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAg                                                                <  2:894459/150‑1 (MQ=255)
                                                                                   cGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGt                                                               >  2:1397057/1‑150 (MQ=255)
                                                                                   cGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGt                                                               >  2:1411053/1‑150 (MQ=255)
                                                                                     cacaTGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCTCAGGAAGTTTCAATACCCAGTac                                                             <  2:930866/150‑1 (MQ=255)
                                                                                     cacaTGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTac                                                             >  1:633553/1‑150 (MQ=255)
                                                                                      acaTGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTaca                                                            >  1:375756/1‑150 (MQ=255)
                                                                                          gATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCa                                                        <  2:1200017/150‑1 (MQ=255)
                                                                                          gATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTACACCTGACTATACAATTGGTTGGCTAACACACCTTTTTCAGCGTCGTAAATGGCGATGCCGGTGTCATCGCAGGAAGGTTCAATACCCAGTACACGCa                                                        >  1:1434957/1‑150 (MQ=255)
                                                                                            tCGCGGGTGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATg                                                      <  2:74797/150‑1 (MQ=255)
                                                                                             cgcgGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGa                                                     >  2:226947/1‑150 (MQ=255)
                                                                                              gcgGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGAc                                                    <  1:1195585/150‑1 (MQ=255)
                                                                                                   aGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGtt                                               <  1:307807/150‑1 (MQ=255)
                                                                                                    ggCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGttt                                              >  1:130127/1‑150 (MQ=255)
                                                                                                    ggCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGttt                                              >  1:1233723/1‑150 (MQ=255)
                                                                                                     gCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGtttt                                             <  2:1333491/150‑1 (MQ=255)
                                                                                                     gCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGtttt                                             <  2:273802/150‑1 (MQ=255)
                                                                                                       cAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTAc                                           >  1:428755/1‑150 (MQ=255)
                                                                                                        aGTTCATGGACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTAcc                                          <  1:1327807/150‑1 (MQ=255)
                                                                                                               gCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCttt                                   <  2:735979/150‑1 (MQ=255)
                                                                                                                cACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTtta                                  <  2:237470/150‑1 (MQ=255)
                                                                                                                cACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTtta                                  >  2:1248526/1‑150 (MQ=255)
                                                                                                                  cgacgCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTttatt                                >  2:102181/1‑150 (MQ=255)
                                                                                                                  cgacgCCGCCGTAGTCAGCGTGCAATTTCAACTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTttatt                                >  2:237502/1‑150 (MQ=255)
                                                                                                                   gacgCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTttatta                               <  1:649876/150‑1 (MQ=255)
                                                                                                                      aattgtatAGTCAGGTGAAATTGCACGCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTTTATTACcg                            >  1:189311/8‑150 (MQ=255)
                                                                                                                       attgtatAGTCAGGTGAAATTGCACGCTGCCTATACAATTGGTTGGCTAACACCCCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCTCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTTTATTACcgc                           >  1:1520255/7‑150 (MQ=255)
                                                                                                                        cgccgTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGNCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTTTATTACcgcg                          <  1:1233957/150‑1 (MQ=255)
                                                                                                                         gccgTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACATCGCTTTATTACCgcgc                         >  1:339268/1‑150 (MQ=255)
                                                                                                                            atAGTCAGGTGAAATTGCACGCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTAGCTTTATTACCGCGCAgt                      >  1:647400/2‑150 (MQ=255)
                                                                                                                             tAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTTTATTACCGCGCAgtg                     <  1:1411053/150‑1 (MQ=255)
                                                                                                                               gTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCttt                                   <  2:614710/134‑1 (MQ=255)
                                                                                                                               gTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCttt                                   >  1:614710/1‑134 (MQ=255)
                                                                                                                               gTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTTTATTACCGCGCAGTGTa                   <  1:229053/150‑1 (MQ=255)
                                                                                                                               gTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTTTATTACCGCGCAGTGTa                   <  1:768722/150‑1 (MQ=255)
                                                                                                                                 cAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTTTATTACCGCGCAGTGTAgg                 >  2:758780/1‑150 (MQ=255)
                                                                                                                                    cGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTTTATTACCGCGCAGTGTAGGAcc              <  1:324284/150‑1 (MQ=255)
                                                                                                                                      tGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTTTATTACCGCGCAGTGTa                   >  2:893355/1‑143 (MQ=255)
                                                                                                                                      tGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTTTATTACCGCGCAGTGTa                   <  1:893355/143‑1 (MQ=255)
                                                                                                                                      tGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTTTATTACCGCGCAGTGTAGGACCaa            >  2:641227/1‑150 (MQ=255)
                                                                                                                                             tCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATAGTCGTAAATGGCGATCCCGCTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCAAGACTTGTTTTACCTCGCTTTATTACCGAGCAGTGTAGGACCAATGAGGGt     >  1:1003021/1‑150 (MQ=255)
                                                                                                                                              cACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTTTATTACCGCGCAGTGTAGGACCAATGCGGGtt    >  2:459004/1‑150 (MQ=255)
                                                                                                                                                ccTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTTTATTACCGCGCAGTGTAGGACCAATGCGGGTTGa  >  1:593036/1‑150 (MQ=255)
                                                                                                                                                ccTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATAGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTTTATTACCGCGCAGTGTAGGACAAATGCGGGTTGa  >  2:647815/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CTGCGGTATAGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTTTATTACCGCGCAGTGTAGGACCAATGCGGGTTGA  >  NC_000913/3210303‑3210596

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