Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,668,181 C→T 5.4% T336T (ACG→ACA)  ccp ← cytochrome c peroxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,668,1810CT 5.4% 294.4 / 3.3 93T336T (ACG→ACA) ccpcytochrome c peroxidase
Reads supporting (aligned to +/- strand):  ref base C (52/35);  new base T (3/2);  total (56/37)
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.20e-01

GACCAATATCCGCCGCCGTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCG  >  NC_000913/3668033‑3668329
                                                                                                                                                    |                                                                                                                                                    
gACCAATATCCGCCGCCGTAATTTCCCCAAAGTTAAAGTCTTTTTTCATCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATTGCCTTTTTTCTGTTGCGCCg                                                                                                                                                     <  2:548908/150‑1 (MQ=255)
gACCAATATCCGCCGCCGTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCg                                                                                                                                                     <  2:1090953/150‑1 (MQ=255)
  ccAATATCCGCCGCCGTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGGCAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTc                                                                                                                                                   >  1:1328037/1‑150 (MQ=255)
  ccAATATCCGCCGCCGTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTc                                                                                                                                                   >  2:338897/1‑150 (MQ=255)
    aaTATCCGCCGCCGTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAg                                                                                                                                                 >  2:904055/1‑150 (MQ=255)
    aaTATCCGCCGCCGTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAg                                                                                                                                                 >  1:933010/1‑150 (MQ=255)
    aaTATCCGCCGCCGTAATTTCCCCAAAGTTAAAGTCTTTCTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATGCCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAg                                                                                                                                                 <  2:1449923/150‑1 (MQ=255)
            ccgccgTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCTGACGGCGCaacag                                                                                                                                         <  2:200734/150‑5 (MQ=255)
            ccgccgTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCTGACGGCGCaacag                                                                                                                                         <  1:634473/150‑5 (MQ=255)
                cgTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCAtc                                                                                                                                     <  2:1391907/150‑1 (MQ=255)
                   aaTTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTcc                                                                                                                                  >  2:1570136/1‑150 (MQ=255)
                       tCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCgcgc                                                                                                                              <  1:1101998/150‑1 (MQ=255)
                        ccccAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCa                                                                                                                             >  2:596682/1‑150 (MQ=255)
                        ccccAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCa                                                                                                                             >  1:772059/1‑150 (MQ=255)
                        ccccAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCa                                                                                                                             <  1:1246513/150‑1 (MQ=255)
                        ccccAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCa                                                                                                                             >  1:1205315/1‑150 (MQ=255)
                        ccccAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCa                                                                                                                             <  2:11837/150‑1 (MQ=255)
                         cccAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTACGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATATTTGATAGCCTTTTTTCTGTTGCGCCTTCAGCGCATTTTCATCTCCGCGCaa                                                                                                                            >  2:1611019/1‑150 (MQ=255)
                                   aaGTCTTTTTTCAGCCCCAACGGTTCAAAGGACCGTCCGCCGAGAATAATACCACCATGAAAAGTTGCACATTTATTATCTTTAAATAATTGATAGCGTTTTTTCTGTTGCGCATTCAGCGAATTTTCATCTGCGCGCAACCATTTATCa                                                                                                                  >  2:656435/1‑150 (MQ=255)
                                    aGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTa                                                                                                                     >  1:497097/1‑146 (MQ=255)
                                    aGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTa                                                                                                                     <  2:497097/146‑1 (MQ=255)
                                      tCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAAt                                                                                                               >  1:1038560/1‑150 (MQ=255)
                                       cTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATg                                                                                                              >  2:786915/1‑150 (MQ=255)
                                       cTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATg                                                                                                              >  2:126233/1‑150 (MQ=255)
                                         tttttCAGCCCCAACGGTTCAAAGGAACGTCCTCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATggg                                                                                                            <  1:596682/150‑1 (MQ=255)
                                         tttttCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATggg                                                                                                            <  1:338897/150‑1 (MQ=255)
                                          ttttCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATgggg                                                                                                           <  2:1034398/150‑1 (MQ=255)
                                             tCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATggg                                                                                                            <  1:215051/146‑1 (MQ=255)
                                             tCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATggg                                                                                                            >  2:215051/1‑146 (MQ=255)
                                                gCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACACGTTGCACATTTTTTATCTTTAAATAATTGATAGCCTGTTTTCTGTTGCGCCGTCAGCGCATTTTCTGCTCCGCGCAACCATTTAGCTAATGGGGAATcgg                                                                                                     >  2:1593616/1‑148 (MQ=255)
                                                gCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCg                                                                                                     >  1:894875/1‑150 (MQ=255)
                                                  cccAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGc                                                                                                   <  2:534017/150‑1 (MQ=255)
                                                  cccAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGc                                                                                                   <  1:1337138/150‑1 (MQ=255)
                                                  cccAACGGTTCAAAGGAACGTCCGCAGAGAAAAATCCCACCAGGACAAGTTGGCCATTTATTATCTTTTAATAATTGATGGCCATTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGGGGACCGATTTATCAAAAGGGGAATCGGGc                                                                                                   >  2:1498787/1‑150 (MQ=255)
                                                     aaCGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGta                                                                                                <  1:823924/150‑1 (MQ=255)
                                                          ttCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGtaattaa                                                                                           >  1:726685/1‑150 (MQ=255)
                                                           tCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGtaattaat                                                                                          <  1:904055/150‑1 (MQ=255)
                                                             aaaGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGt                                                                                        >  1:1292176/1‑150 (MQ=255)
                                                               aGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGttt                                                                                      <  2:111285/150‑1 (MQ=255)
                                                                  aaCGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTtct                                                                                   <  2:44090/150‑1 (MQ=255)
                                                                    cGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCa                                                                                 >  1:187310/1‑150 (MQ=255)
                                                                         gccgAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAAtta                                                                            >  1:1074057/1‑149 (MQ=255)
                                                                         gccgAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTc                                                                            <  1:843979/150‑1 (MQ=255)
                                                                         gccgAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTc                                                                            >  2:771562/1‑150 (MQ=255)
                                                                          gcgAGAATAATACCACCATTACAAGTTGCACATTTATTATCTTTAAATTATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCa                                                                           <  2:727563/149‑1 (MQ=255)
                                                                                ataataCCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATg                                                                     <  1:153511/150‑1 (MQ=255)
                                                                                 taataCCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATgg                                                                    >  1:766881/1‑150 (MQ=255)
                                                                                   ataCCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTACGCGCAACCATTTATAAAATGGGGAATCCGGAGTAATTAATGTTTTCTCAAATACAGCAATGCCa                                                                  >  1:1524305/1‑150 (MQ=255)
                                                                                    taCCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCAt                                                                 >  2:258802/1‑150 (MQ=255)
                                                                                      ccaccaTGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTATTTTGCGCCGACAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCa                                                               >  2:706147/1‑150 (MQ=255)
                                                                                       caccaTGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAg                                                              >  2:791300/1‑150 (MQ=255)
                                                                                        accaTGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGt                                                             >  2:969413/1‑150 (MQ=255)
                                                                                         ccaTGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTa                                                            >  2:1041448/1‑150 (MQ=255)
                                                                                           aTGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAat                                                          >  2:370908/1‑150 (MQ=255)
                                                                                                  gTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCg                                                   >  2:539336/1‑150 (MQ=255)
                                                                                                   ttGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGc                                                  <  1:168268/150‑1 (MQ=255)
                                                                                                    tGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGcc                                                 <  1:1570136/150‑1 (MQ=255)
                                                                                                    tGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGcc                                                 >  2:552689/1‑150 (MQ=255)
                                                                                                     gCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCa                                                >  2:616288/1‑150 (MQ=255)
                                                                                                     gCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCa                                                >  2:1534134/1‑150 (MQ=255)
                                                                                                      cacaTTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCaa                                               >  1:974006/1‑149 (MQ=255)
                                                                                                        caTTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTg                                             <  1:258802/150‑1 (MQ=255)
                                                                                                        caTTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGAGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAAAAATGCGTAATTAATGTTTTGTCAAATTAAGCAATGGCTTCAGTAATTTTTTctctaaaa                                             >  1:716048/1‑143 (MQ=255)
                                                                                                              ttatCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCt                                       <  1:786915/150‑1 (MQ=255)
                                                                                                                atCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATATCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTg                                     >  1:1403399/1‑150 (MQ=255)
                                                                                                                  cTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCg                                   >  2:1575779/1‑150 (MQ=255)
                                                                                                                   tttAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCgg                                  >  1:764215/1‑150 (MQ=255)
                                                                                                                   tttAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCgg                                  >  2:806835/1‑150 (MQ=255)
                                                                                                                    ttAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGa                                 >  2:1411217/1‑150 (MQ=255)
                                                                                                                       aataatTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAg                              >  1:1356972/1‑150 (MQ=255)
                                                                                                                       aataatTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCAGCGCAACCATTTATCAAATTGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAg                              >  1:1507134/1‑150 (MQ=255)
                                                                                                                        ataatTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGa                             >  1:966421/1‑150 (MQ=255)
                                                                                                                         taatTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCa                                                >  1:1335404/1‑130 (MQ=255)
                                                                                                                         taatTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCa                                                <  2:1335404/130‑1 (MQ=255)
                                                                                                                         taatTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGAc                            <  2:1189425/150‑1 (MQ=255)
                                                                                                                          aatTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACt                           <  2:894875/150‑1 (MQ=255)
                                                                                                                          aatTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACt                           <  2:1038560/150‑1 (MQ=255)
                                                                                                                              gATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGa                       <  2:300218/150‑1 (MQ=255)
                                                                                                                              gATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACATTGCGGATAGACTTCGa                       >  2:1577381/1‑150 (MQ=255)
                                                                                                                                 aGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGa                    >  1:279526/1‑150 (MQ=255)
                                                                                                                                 aGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGa                    <  1:44304/150‑1 (MQ=255)
                                                                                                                                  gCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGaa                   >  2:701260/1‑150 (MQ=255)
                                                                                                                                   ccTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAAc                  <  2:1328037/150‑1 (MQ=255)
                                                                                                                                         tgaaaatGCGCTGACGGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGtt            >  2:1114919/7‑150 (MQ=255)
                                                                                                                                          gaaaatGCGCTGACGGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGttt           >  1:350139/6‑150 (MQ=255)
                                                                                                                                           cTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGtttt          >  2:1005520/1‑150 (MQ=255)
                                                                                                                                           aaaatGCGCTGACGGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGtttt          >  1:1129593/5‑150 (MQ=255)
                                                                                                                                            tGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTa         <  1:91429/150‑1 (MQ=255)
                                                                                                                                            tGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTa         >  1:729101/1‑150 (MQ=255)
                                                                                                                                            tGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTa         >  1:250230/1‑150 (MQ=255)
                                                                                                                                             gTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTaa        >  1:20632/1‑150 (MQ=255)
                                                                                                                                               tGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGc      <  1:1041448/150‑1 (MQ=255)
                                                                                                                                                  gcCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGc   <  1:891870/150‑1 (MQ=255)
                                                                                                                                                   ccGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCg  <  1:806835/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GACCAATATCCGCCGCCGTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCG  >  NC_000913/3668033‑3668329

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