Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,518,197 C→G 8.1% noncoding (35/76 nt) alaW ← tRNA‑Ala

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,518,1970CG 8.1% 426.9 / 23.1 123noncoding (35/76 nt)alaWtRNA‑Ala
Reads supporting (aligned to +/- strand):  ref base C (58/55);  new base G (5/5);  total (63/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

Only 100 of 129 total aligned reads displayed.

GATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAAC  >  NC_000913/2518057‑2518343
                                                                                                                                            |                                                                                                                                                   
gATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAg                                                                                                                                             >  3:527629/1‑149 (MQ=255)
    gAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCTCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATAGAACCGCTGACCTCTTGCATGCCATGCAAGCGct                                                                                                                                         >  1:1193350/1‑149 (MQ=255)
    gAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTGGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGct                                                                                                                                         <  2:1439958/149‑1 (MQ=255)
     aaCCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatgcca                                                                                                                                                   <  1:581795/138‑1 (MQ=255)
     aaCCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatgcca                                                                                                                                                   >  2:581795/1‑138 (MQ=255)
     aaCCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAgc                                                                                                                                            <  1:1325531/145‑1 (MQ=255)
     aaCCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAgc                                                                                                                                            >  2:1325531/1‑145 (MQ=255)
               ctctTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGa                                                                                                                              <  1:733185/149‑1 (MQ=255)
               ctctTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGa                                                                                                                              >  4:275431/1‑149 (MQ=255)
               ctctTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGa                                                                                                                              >  2:733185/1‑149 (MQ=255)
               ctctTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCAATTACGTACCAAGTTTGATGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACTGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGa                                                                                                                              >  3:780091/1‑149 (MQ=255)
                  ttGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCt                                                                                                                           >  3:712184/1‑149 (MQ=255)
                   tGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCta                                                                                                                          >  2:1145902/1‑149 (MQ=255)
                          catgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGc                                                                                                                            >  3:945619/1‑140 (MQ=255)
                           atgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCg                                                                                                                  >  3:226035/1‑149 (MQ=255)
                             gcAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTc                                                                                                                >  3:617618/1‑149 (MQ=255)
                                aGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCCGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACg                                                                                                             >  2:122673/1‑149 (MQ=255)
                                aGCGCTCTCCAAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCGCTTGCATGGCATGCAAga                                                                                                                                            >  2:321780/1‑117 (MQ=255)
                                aGCGCTCTCCAAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCGCTTGCATGGCATGCAAga                                                                                                                                            <  1:321780/118‑2 (MQ=37)
                                 gcgcTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGt                                                                                                            <  4:712184/149‑1 (MQ=255)
                                         ccAGCTGAGCTATAGCCGCACGATGCGTTTACGTACCAAGGTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTg                                                                                                    >  4:927654/1‑149 (MQ=255)
                                            gCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGGCATGCAAga                                                                                                                                            >  4:1158300/1‑105 (MQ=255)
                                            gCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCg                                                                                                 >  2:1452610/1‑149 (MQ=255)
                                            gCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCg                                                                                                 >  4:239686/1‑149 (MQ=255)
                                             cTGATCTATAACGCCACGATGCGTTTACGTCCCAAGTTTGCAGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCCATCCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGa                                                                                                <  4:780091/149‑1 (MQ=21)
                                              tGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTg                                                                                                    <  3:1204191/144‑1 (MQ=255)
                                              tGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTg                                                                                                    >  4:1204191/1‑144 (MQ=255)
                                              tGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAg                                                                                               <  4:576300/149‑1 (MQ=255)
                                              tGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAg                                                                                               >  4:12277/1‑149 (MQ=255)
                                               gAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGt                                                                                              <  1:1452610/149‑1 (MQ=255)
                                                   tataGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCtata                                                                                                                        <  4:1453264/119‑1 (MQ=255)
                                                   tataGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGAc                                                                                          >  2:1272002/1‑149 (MQ=255)
                                                    ataGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACg                                                                                         >  2:1207674/1‑149 (MQ=255)
                                                    ataGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACg                                                                                         <  1:122673/149‑1 (MQ=255)
                                                     taGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGt                                                                                              <  2:951830/143‑1 (MQ=255)
                                                      aGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGGCATGCAAga                                                                                                                                            <  1:1295914/96‑2 (MQ=255)
                                                      aGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGGCATGCAAga                                                                                                                                            >  2:1295914/1‑95 (MQ=255)
                                                        ccccACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGggcg                                                                                     <  4:617618/149‑1 (MQ=255)
                                                         cccACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCaa                                                                                                                                              <  1:724203/91‑1 (MQ=255)
                                                           cACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTccc                                                                                                                                    <  3:444669/99‑1 (MQ=255)
                                                           cACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTccc                                                                                                                                    >  4:444669/1‑99 (MQ=255)
                                                           cACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGca                                                                                  >  3:233281/1‑149 (MQ=255)
                                                                tGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCata                                                                            <  3:275431/149‑1 (MQ=255)
                                                                   gTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGa                                                                         >  3:611497/1‑149 (MQ=255)
                                                                   gTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGa                                                                         >  3:1106949/1‑149 (MQ=255)
                                                                          tACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGc                                                                  >  1:1094626/1‑149 (MQ=255)
                                                                           aCCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGcc                                                                 <  1:1272002/149‑1 (MQ=255)
                                                                            ccAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCt                                                                                                                           <  4:1076564/91‑1 (MQ=255)
                                                                            ccAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCt                                                                                                                           >  3:1076564/1‑91 (MQ=255)
                                                                            ccAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGca‑tc                                                                               <  4:683124/134‑1 (MQ=255)
                                                                            ccAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGccc                                                                <  2:1094626/149‑1 (MQ=255)
                                                                             cAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatgccatgc                                                                                                                                                >  4:554295/1‑69 (MQ=255)
                                                                             cAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatgccatgc                                                                                                                                                <  3:554295/69‑1 (MQ=255)
                                                                             cAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGa                                                                                                                              >  3:130005/1‑87 (MQ=255)
                                                                               aGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGaa                                                             <  1:1145902/149‑1 (MQ=255)
                                                                               aGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGaa                                                             <  3:12277/149‑1 (MQ=255)
                                                                                  ttGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAAtgt                                                          <  1:1207674/149‑1 (MQ=255)
                                                                                   tGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAAtgtg                                                         >  1:1414860/1‑149 (MQ=255)
                                                                                     cTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTc                                                       <  4:611497/149‑1 (MQ=255)
                                                                                      tGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCa                                                      >  1:84068/1‑149 (MQ=255)
                                                                                      tGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCa                                                      >  1:558765/1‑149 (MQ=255)
                                                                                       gggTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCaa                                                     <  2:558765/149‑1 (MQ=255)
                                                                                            cAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGGCATGCAAga                                                                                                                                            >  4:694355/1‑57 (MQ=255)
                                                                                            cAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGGCATGCAAga                                                                                                                                            <  3:694355/58‑2 (MQ=255)
                                                                                                    ggtggAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCaa                                        <  4:226035/149‑1 (MQ=255)
                                                                                                    ggtggAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCaa                                        <  3:239686/149‑1 (MQ=255)
                                                                                                     gtggAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAAc                                       <  4:88871/149‑1 (MQ=255)
                                                                                                     gtggAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAAc                                       <  2:732137/149‑1 (MQ=255)
                                                                                                     gtggAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAAc                                       <  3:927654/149‑1 (MQ=255)
                                                                                                      tggAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGc                                                                  >  3:1226496/1‑121 (MQ=255)
                                                                                                      tggAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGc                                                                  <  4:1226496/121‑1 (MQ=255)
                                                                                                        gAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACgtg                                    <  4:1106949/149‑1 (MQ=255)
                                                                                                          gCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTa                                  <  2:1414860/149‑1 (MQ=255)
                                                                                                           cTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAg                                 <  3:681126/149‑1 (MQ=255)
                                                                                                            tAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGAATTCATATGatgc                                                                      >  1:271191/1‑109 (MQ=18)
                                                                                                                cGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGt                                                                                                   <  4:1355624/79‑1 (MQ=255)
                                                                                                                cGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGt                                                                                                   >  3:1355624/1‑79 (MQ=255)
                                                                                                                  ggATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAAt                          <  2:84068/149‑1 (MQ=255)
                                                                                                                   gATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATc                         >  3:1318361/1‑149 (MQ=255)
                                                                                                                    aTCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGaa                                                             <  1:775451/112‑1 (MQ=255)
                                                                                                                    aTCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGaa                                                             >  2:775451/1‑112 (MQ=255)
                                                                                                                      cGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCt                      <  4:1318361/149‑1 (MQ=255)
                                                                                                                        aaCCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGa                    >  1:772981/1‑149 (MQ=255)
                                                                                                                        aaCCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGa                    >  1:1418656/1‑149 (MQ=255)
                                                                                                                           cGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCTCGTCACGTAAAGCTTGTCGAGTTGACGGGCAGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCaa                           <  2:1191112/139‑1 (MQ=255)
                                                                                                                           cGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCCCGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCaa                           >  1:1191112/1‑139 (MQ=255)
                                                                                                                             cTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAAt               >  2:1194191/1‑149 (MQ=255)
                                                                                                                              tGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATc              >  2:847245/1‑149 (MQ=255)
                                                                                                                               gACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCa             >  4:117467/1‑149 (MQ=255)
                                                                                                                               gACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCa             <  3:769489/149‑1 (MQ=255)
                                                                                                                                aCCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAg            <  3:117467/149‑1 (MQ=255)
                                                                                                                                  ctctTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCg                        >  4:105457/1‑135 (MQ=255)
                                                                                                                                  ctctTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGc          >  1:695755/1‑149 (MQ=255)
                                                                                                                                   tctTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACgtgt                                   >  2:1197915/1‑123 (MQ=255)
                                                                                                                                   gctTGCATGGCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGa                                                                                                >  3:1071324/2‑63 (MQ=255)
                                                                                                                                     ttGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAg                                                                                                        >  1:526594/1‑53 (MQ=255)
                                                                                                                                     ttGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAg                                                                                                        <  2:526594/53‑1 (MQ=255)
                                                                                                                                     ttGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCaaa       >  3:1038806/1‑149 (MQ=255)
                                                                                                                                          tgccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAAc  >  1:1256797/1‑149 (MQ=255)
                                                                                                                                          tgccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAAc  >  1:1275312/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
GATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAAC  >  NC_000913/2518057‑2518343

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: ‑