Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,210,444 T→G 5.2% K34Q (AAA→CAA)  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,4440TG 5.2% 311.7 / 7.1 97K34Q (AAA→CAA) tsaDN(6)‑L‑threonylcarbamoyladenine synthase, TsaD subunit
Reads supporting (aligned to +/- strand):  ref base T (41/50);  new base G (3/2);  total (45/52)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

CTGCGGTATAGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTTTATTACCGCGCAGTGTAGGACCAATGCGGGT  >  NC_000913/3210303‑3210593
                                                                                                                                             |                                                                                                                                                      
cTGCGGTATAGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGAc                                                                                                                                               <  1:904081/150‑1 (MQ=255)
 tGCGGTATAGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACt                                                                                                                                              <  1:525477/150‑1 (MQ=255)
 tGCGGTATAGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACt                                                                                                                                              <  1:1114035/150‑1 (MQ=255)
 tGCGGTATAGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGGTGAAATTGCACGCTGACt                                                                                                                                              <  2:401166/150‑1 (MQ=255)
   cGGTATAGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACtat                                                                                                                                            >  1:187260/1‑150 (MQ=255)
    ggTATAGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCAc                                                                                                                                                     >  1:490523/1‑141 (MQ=255)
    ggTATAGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCAc                                                                                                                                                     <  2:490523/141‑1 (MQ=255)
      tataGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACa                                                                                                                                         >  1:1330004/1‑150 (MQ=255)
       ataGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGGTGAAATTGCACGCTGACtacggcg                                                                                                                                        <  2:1276480/150‑6 (MQ=255)
           gCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAAttgg                                                                                                                                    >  2:105542/1‑150 (MQ=255)
                 gCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCt                                                                                                                              >  1:916347/1‑150 (MQ=255)
                 gCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCt                                                                                                                              <  2:916347/150‑1 (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                                                                                                  <  1:470101/132‑1 (MQ=255)
                                                               ggATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATgg                                                                                                  >  2:470101/1‑132 (MQ=255)
                                                                 aTCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAg                                                                              >  1:1333491/1‑150 (MQ=255)
                                                                 aTCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAg                                                                              <  2:1505221/150‑1 (MQ=255)
                                                                 aTCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAg                                                                              <  1:20250/150‑1 (MQ=255)
                                                                   cAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGa                                                                            >  1:1083491/1‑150 (MQ=255)
                                                                   cAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGa                                                                            >  1:43459/1‑150 (MQ=255)
                                                                    aaCGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCa                                                                                                                                                      >  1:208673/1‑76 (MQ=255)
                                                                    aaCGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCa                                                                                                                                                      <  2:208673/76‑1 (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:273802/150‑1 (MQ=255)
                                                                                                     gCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGtttt                                          <  2:1333491/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:1248526/1‑150 (MQ=255)
                                                                                                                cACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTtta                               <  2:237470/150‑1 (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                                >  1:614710/1‑134 (MQ=255)
                                                                                                                               gTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCttt                                <  2:614710/134‑1 (MQ=255)
                                                                                                                               gTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTTTATTACCGCGCAGTGTa                <  1:768722/150‑1 (MQ=255)
                                                                                                                               gTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTTTATTACCGCGCAGTGTa                <  1:229053/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)
                                                                                                                                             |                                                                                                                                                      
CTGCGGTATAGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTTTATTACCGCGCAGTGTAGGACCAATGCGGGT  >  NC_000913/3210303‑3210593

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: