Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,307,070 G→T 9.3% A264D (GCT→GAT)  deaD ← ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,307,0700GT 9.3% 326.3 / 18.9 97A264D (GCT→GAT) deaDATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base G (42/46);  new base T (6/3);  total (48/49)
Fisher's exact test for biased strand distribution p-value = 3.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.70e-01

ACGCTCAACGTCCAGGCCACGGGCTGCAACGTCGGTCGCAATCAGGATGTCCAGACGACCATCTTTCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTCGTTTTTGCGCATACCCCAGACAGTCCAGTAGCTCTGGCTGATGTCAGGACG  >  NC_000913/3306922‑3307212
                                                                                                                                                    |                                                                                                                                              
aCGCTCAACGTCCAGGCCACGGGCTGCAACGTCGGTCGCAATCAGGATGTCCAGACGACCATCTTTCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAagag                                                                                                                                                <  2:847829/149‑1 (MQ=255)
 cGCTCAACGTCCAGGCCACGGGCTGCAACGTCGGTCGCAATCAGGATGTCCAGACGACCATCTTTCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGc                                                                                                                                               <  4:123045/149‑1 (MQ=255)
  gCTCAACGTCCAGGCCACGGGCTGCAACGTCGGTCGCAATCAGGATGTCCAGACGACCATCTTTCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCt                                                                                                                                              >  1:303493/1‑149 (MQ=255)
  gCTCAACGTCCAGGCCACGGGCTGCAACGTCGGTCGCAATCAGGATGTCCAGACGACCATCTTTCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCt                                                                                                                                              >  3:332758/1‑149 (MQ=255)
   cTCAACGTCCAGGCCACGGGCTGCAACGTCGGTCGCAATCAGGATGTCCAGACGACCATCTTTCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGATCtt                                                                                                                                             >  4:831019/1‑149 (MQ=255)
   cTCAACGTCCAGGCCACGGGCTGCAACGTCGGTCGCAATCAGGATGTCCAGACGACCATCTTTCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCtt                                                                                                                                             <  2:128728/149‑1 (MQ=255)
    tCAACGTCCAGGCCACGGGCTGCAACGTCGGTCGCAATCAGGATGTCCAGACGACCATCTTTCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTc                                                                                                                                            <  2:876441/149‑1 (MQ=255)
     cAACGTCCAGGCCACGGGCTGCAACGTCGGTCGCAATCAGGATGTCCAGACGACCATCTTTCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCCTCa                                                                                                                                           <  4:539723/149‑1 (MQ=255)
      aaCGTCCAGGCCACGGGCTGCAACGTCGGTCGCAATCAGGATGTCCAGACGACCATCTTTCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAg                                                                                                                                          >  1:465227/1‑149 (MQ=255)
             aGGCCACGGGCTGCAACGTCGGTCGCAATCAGGATGTCCAGACGACCATCTTTCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTc                                                                                                                                   <  3:127912/149‑1 (MQ=255)
              ggCCACGGGCTGCAACGTCGGTCGCAATCAGGATGTCCAGACGACCATCTTTCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGATCTTCAGCCACTTcc                                                                                                                                  >  3:375239/1‑149 (MQ=255)
               gCCACGGGCTGCAACTTTGGTCGCAATCAGGATGTGCAGACCACCATCTATCAGGCGTTCCAGTGTCTGGTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGACGTTACGCTCAAGAGCTTGAGCCACTTCCa                                                                                                                                 >  4:571635/1‑149 (MQ=255)
                 cACGGGCTGCAACGTCGGTCGCAATCAGGATGTCCAGACGACCATCTTTCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCaga                                                                                                                               >  1:176962/1‑149 (MQ=255)
                 cACGGGCTGCAACGTCGGTCGCAATCAGGATGTCCAGACGACCATCTTTCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCaga                                                                                                                               >  3:713349/1‑149 (MQ=255)
                 cACGGGCTGCAACGTCGGTCGCAATCAGGATGTCCAGACGACCATCTTTCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCaga                                                                                                                               <  3:383241/149‑1 (MQ=255)
                        tGCAACGTCGGTCGCAATCAGGATGTCCAGACGACCATCTTTCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGt                                                                                                                        <  4:832519/149‑1 (MQ=255)
                         gCAACGTCGGTCGCAATCAGGATGTCCAGACGACCATCTTTCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGtt                                                                                                                       <  3:555150/149‑1 (MQ=255)
                            aCGTCGGTCGCAATCAGGATGTCCAGACGACCATCTTTCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGttttt                                                                                                                    <  4:329064/149‑1 (MQ=255)
                            aCGTCGGTCGCAATCAGGATGTCCAGACGACCATCTTTCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGttttt                                                                                                                    <  2:176962/149‑1 (MQ=255)
                            aCGTCGGTCGCAATCAGGATGTCCAGACGACCATCTTTCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGttttt                                                                                                                    <  4:603807/149‑1 (MQ=255)
                                      cAATCAGGATGTCCAGACGACCATCTTTCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTacgaacg                                                                                                          >  4:787573/1‑149 (MQ=255)
                                        aTCAGGATGTCCAGACGACCATCTTTCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTAcgaacgaa                                                                                                        <  4:601468/149‑1 (MQ=255)
                                        aTCAGGATGTCCAGACGACCATCTTTCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTAcgaacgaa                                                                                                        <  2:360797/149‑1 (MQ=255)
                                         tCAGGATGTCCAGACGACCATCTTTCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAg                                                                                                       >  4:492749/1‑149 (MQ=255)
                                         tCAGGATGTCCAGACGACCATCTTTCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAg                                                                                                       >  4:727201/1‑149 (MQ=255)
                                            ggATGTCCAGACCACCATCTTTCAGGCGTTCCAGTGTCTGTTAACGCAGCGCCTGGGTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATa                                                                                                    >  4:83198/1‑149 (MQ=255)
                                             gATGTCCAGACGACCATCTTTCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATaa                                                                                                   >  3:73663/1‑149 (MQ=255)
                                                   cAGACGACCATCTTTCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATcgccg                                                                                             >  2:844586/1‑149 (MQ=255)
                                                    aGACGACCATCTTTCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATcgccgc                                                                                            >  2:906847/1‑149 (MQ=255)
                                                        gacCATCTTTCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCa                                                                                        >  3:488530/1‑149 (MQ=255)
                                                           cATCTTTCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCaaaa                                                                                     <  1:561773/149‑1 (MQ=255)
                                                           cATCTTTCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCaaaa                                                                                     >  1:705518/1‑149 (MQ=255)
                                                            aTCTTTCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTATGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAAt                                                                                    >  1:122309/1‑149 (MQ=255)
                                                            aTCTTTCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAAt                                                                                    >  4:608221/1‑149 (MQ=255)
                                                             tCTTTCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATc                                                                                   >  1:59927/1‑149 (MQ=255)
                                                               tttCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGATCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCtt                                                                                 >  3:88856/1‑149 (MQ=255)
                                                               tttCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCtt                                                                                 >  1:794163/1‑149 (MQ=255)
                                                                ttCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATc                                                                                                 >  3:667064/1‑132 (MQ=255)
                                                                ttCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATc                                                                                                 <  4:667064/132‑1 (MQ=255)
                                                                ttCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTc                                                                                <  4:332758/149‑1 (MQ=255)
                                                                    ggCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCt                                                                            >  2:296089/1‑149 (MQ=255)
                                                                        ttCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCa                                                                        <  4:437939/149‑1 (MQ=255)
                                                                             gtgtCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTATGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGaaa                                                                   <  2:122309/149‑1 (MQ=255)
                                                                             gtgtCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGaaa                                                                   <  2:303493/149‑1 (MQ=255)
                                                                             gtgtCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGaaa                                                                   <  1:756408/149‑1 (MQ=255)
                                                                              tgtCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGATCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAAc                                                                  >  2:371687/1‑149 (MQ=255)
                                                                              tgtCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTAcgaacgaa                                                                                                        >  1:742051/1‑111 (MQ=255)
                                                                              tgtCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTAcgaacgaa                                                                                                        <  2:742051/111‑1 (MQ=255)
                                                                              tgtCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAAc                                                                  <  3:97432/149‑1 (MQ=255)
                                                                              tgtCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAAc                                                                  <  3:727201/149‑1 (MQ=255)
                                                                                 cTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTa                                                               <  4:756882/149‑1 (MQ=255)
                                                                                  tGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTAc                                                              <  2:504242/149‑1 (MQ=255)
                                                                                  tGTTCACGCAGCGCCTGGTTAATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTAc                                                              <  4:888467/149‑1 (MQ=255)
                                                                                      cACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGATCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGt                                                          <  4:88591/149‑1 (MQ=255)
                                                                                      cACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGt                                                          >  4:569593/1‑149 (MQ=255)
                                                                                      cACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGt                                                          >  1:786203/1‑149 (MQ=255)
                                                                                       aCGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTg                                                         >  3:217645/1‑149 (MQ=255)
                                                                                         gCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCt                                                       <  4:217645/149‑1 (MQ=255)
                                                                                           aGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCATTTCCAGAGTCGCGTTTTTTCTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGTTTc                                                     <  4:586694/149‑1 (MQ=255)
                                                                                           aGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTc                                                     <  2:59927/149‑1 (MQ=255)
                                                                                           aGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGCAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTc                                                     <  3:608221/149‑1 (MQ=255)
                                                                                                cTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCg                                                                              >  4:225376/1‑119 (MQ=255)
                                                                                                cTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCg                                                                              <  3:225376/119‑1 (MQ=255)
                                                                                                 tGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTCGttttt                                               <  1:309180/149‑1 (MQ=255)
                                                                                                 tGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTCGttttt                                               <  3:83198/149‑1 (MQ=255)
                                                                                                  ggTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTCGTTTTTg                                              <  2:465227/149‑1 (MQ=255)
                                                                                                      cATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTCGTTTTTGCGCa                                          >  3:427570/1‑149 (MQ=255)
                                                                                                         gTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTCGTTTTTGCGCATAc                                       >  4:403278/1‑149 (MQ=255)
                                                                                                           cACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTCGTTTTTGCGCATAccc                                     >  2:49669/1‑149 (MQ=255)
                                                                                                           cACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTCGTTTTTGCGCATAccc                                     <  2:786203/149‑1 (MQ=255)
                                                                                                            aCCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTCGTTTTTGCGCATAcccc                                    >  1:491709/1‑149 (MQ=255)
                                                                                                             ccGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTCGTTTTTGCGCATACCCCa                                   >  3:172563/1‑149 (MQ=255)
                                                                                                               gTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTCGTTTTTGCGCATACCCCAGa                                 <  3:787573/149‑1 (MQ=255)
                                                                                                                    gcgcGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTCGTTTTTGCGCATACCCCAGACAGTc                            <  3:841961/149‑1 (MQ=255)
                                                                                                                    gcgcGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTCGTTTTTGCGCATACCCCAGACAGTc                            >  3:71561/1‑149 (MQ=255)
                                                                                                                      gcggcgCTGTTGTAGCCGTTACGCTCAAGATCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTCGTTTTTGCGCATACCCCAGACAGTCCa                          <  1:371687/149‑1 (MQ=255)
                                                                                                                        ggcgCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCAATTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTCGTTTTTGCGCATACCCCAGACAGTCCTGt                        <  4:71561/149‑1 (MQ=255)
                                                                                                                           gcTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTCGTTTTTGCGCATACCCCAGACAGTCCAGTAGc                     >  1:794748/1‑149 (MQ=255)
                                                                                                                             tgttgtAGCCGTTACGCTCAAGATCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTCGtt                                                  <  4:285528/118‑1 (MQ=255)
                                                                                                                             tgttgtAGCCGTTACGCTCAAGATCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTCGtt                                                  >  3:285528/1‑118 (MQ=255)
                                                                                                                             tgttgtAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTCGTTTTTGCGCATACCCCAGCCATTCCAGTAGctc                   >  1:932055/1‑149 (MQ=255)
                                                                                                                                tgtAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTCGTTTTTGCGCATACCCCAGACAGTCCAGTAGCTCTgg                <  3:238713/149‑1 (MQ=255)
                                                                                                                                 gtAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGTGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTCGTTTTTGCGCATACCCCAGACAGTCCAGTAGCTCTGGc               >  2:681474/1‑149 (MQ=255)
                                                                                                                                  tAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTCGTTTTTGCGCATACCCCAGACAGTCCAGTAGCTCTGGCt              <  1:906847/149‑1 (MQ=255)
                                                                                                                                  tAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTCGTTTTTGCGCATACCCCAGACAGTCCAGTAGCTCTGGCt              >  1:643366/1‑149 (MQ=255)
                                                                                                                                    gCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTCGTTTTTGCGCATACCCCAGACAGTCCAGTAGCTCTGGCTGa            <  4:427570/149‑1 (MQ=255)
                                                                                                                                    gCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTCGTTTTTGCGCATACCCCAGACAGTCCAGTAGCTCTGGCTGa            >  3:656412/1‑149 (MQ=255)
                                                                                                                                    gCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTCGTTTTTGCGCATACCCCAGACAGTCCAGTAGCTCTGGCTGa            >  2:567231/1‑149 (MQ=255)
                                                                                                                                     ccGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTCGTTTTTGCGCATACCCCAGACAGTCCAGTAGCTCTGGCTGAt           <  2:794748/149‑1 (MQ=255)
                                                                                                                                     ccGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTCGTTTTTGCGCATACCCCAGACAGTCCAGTAGCTCTGGCTGAt           <  2:309737/149‑1 (MQ=255)
                                                                                                                                       gTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTCGTTTTTGCGCATACCCCAGACAGTCCAGTAGCTCTGGCTGATGt         <  3:388950/149‑1 (MQ=255)
                                                                                                                                       gTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTCGTTTTTGCGCATACCCCAGACAGTCCAGTAGCTCTGGCTGATGt         >  4:328651/1‑149 (MQ=255)
                                                                                                                                            gCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTCGTTTTTGCGCATACCCCAGACAGTCCAGTAGCTCTGGCTGATGTCAGGa    >  1:773496/1‑149 (MQ=255)
                                                                                                                                            gCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTCGTTTTTGCGCATACCCCAGACAGTCCAGTAGCTCTGGCTGATGTCAGGa    <  4:172563/149‑1 (MQ=255)
                                                                                                                                            gCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTCGTTTTTGCGCATACCCCAGACAGTCCAGTAGCTCTGGCTGATGTCAGGa    >  3:731065/1‑149 (MQ=255)
                                                                                                                                              tCAAGATCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTCGTTTTTGCGCATACCCCAGACAGTCCAGTAGCTCTGGCTGATGTCAGGACg  >  4:643466/1‑149 (MQ=255)
                                                                                                                                              tCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTCGTTTTTGCGCATACCCCAGACAGTCCAg                         >  2:72988/1‑126 (MQ=255)
                                                                                                                                              tCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTCGTTTTTGCGCATACCCCAGACAGTCCAg                         <  1:72988/126‑1 (MQ=255)
                                                                                                                                                 agagCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTCGTTTTTgc                                             >  2:332176/1‑103 (MQ=255)
                                                                                                                                                 agagCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTCGTTTTTgc                                             <  1:332176/103‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
ACGCTCAACGTCCAGGCCACGGGCTGCAACGTCGGTCGCAATCAGGATGTCCAGACGACCATCTTTCAGGCGTTCCAGTGTCTGTTCACGCAGCGCCTGGTTCATGTCACCGTTCAGCGCGGCGCTGTTGTAGCCGTTACGCTCAAGAGCTTCAGCCACTTCCAGAGTCGCGTTTTTGGTACGAACGAAGATAATCGCCGCATCAAAATCTTCCGCTTCCAGGAAACGTACCAGTGCTTCGTTTTTGCGCATACCCCAGACAGTCCAGTAGCTCTGGCTGATGTCAGGACG  >  NC_000913/3306922‑3307212

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