Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,210,435 Δ1 bp 5.3% coding (109/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,4350C. 5.3% 323.4 / 10.1 95coding (109/1014 nt)tsaDN(6)‑L‑threonylcarbamoyladenine synthase, TsaD subunit
Reads supporting (aligned to +/- strand):  ref base C (41/49);  new base . (3/2);  total (44/51)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.97e-01

CTAATCCAGGGCCTGCGGTATAGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTTTATTACCGCGCAGTGTAGGACC  >  NC_000913/3210291‑3210584
                                                                                                                                                |                                                                                                                                                      
cTAATCCAGGGCCTGCGGTATAGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCa                                                                                                                                                   >  1:735979/1‑150 (MQ=255)
  aaTCCAGGGCCTGCGGTATAGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAAt                                                                                                                                                 <  1:981552/150‑1 (MQ=255)
            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                                                                                                                                            >  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                                                                                                                     <  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                                                                                                >  1:1126944/1‑150 (MQ=255)
                                                  aaaCCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCAtcgtcg                                                                                                <  2:694376/150‑1 (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                                                                                        <  1:1050030/150‑1 (MQ=255)
                                                          cTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGc                                                                                        <  2:169300/150‑1 (MQ=255)
                                                            ccTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGa                                                                                      >  1:960309/1‑150 (MQ=255)
                                                            ccTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGa                                                                                      <  2:1014854/150‑1 (MQ=255)
                                                             cTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGAt                                                                                     >  1:1200017/1‑150 (MQ=255)
                                                             cTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGAt                                                                                     <  2:564341/150‑1 (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:20250/150‑1 (MQ=255)
                                                                             aTCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAg                                                                     <  2:1505221/150‑1 (MQ=255)
                                                                             aTCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAg                                                                     >  1:1333491/1‑150 (MQ=255)
                                                                               cAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGa                                                                   >  1:1083491/1‑150 (MQ=255)
                                                                               cAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGa                                                                   >  1:43459/1‑150 (MQ=255)
                                                                                aaCGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCa                                                                                                                                             <  2:208673/76‑1 (MQ=255)
                                                                                aaCGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCa                                                                                                                                             >  1:208673/1‑76 (MQ=255)
                                                                                 acggtacggtTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAg                                                                 >  2:976480/1‑150 (MQ=255)
                                                                                    gtacggtTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATg                                                                                          >  2:618986/1‑122 (MQ=255)
                                                                                    gtacggtTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATg                                                                                          <  1:618986/122‑1 (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:1233723/1‑150 (MQ=255)
                                                                                                                ggCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGttt                                  >  1:130127/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: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)
                                                                                                                                                |                                                                                                                                                      
CTAATCCAGGGCCTGCGGTATAGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTTTATTACCGCGCAGTGTAGGACC  >  NC_000913/3210291‑3210584

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: