Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,668,185 A→G 5.5% L335P (CTG→CCG)  ccp ← cytochrome c peroxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,668,1850AG 5.5% 283.1 / 4.0 91L335P (CTG→CCG) ccpcytochrome c peroxidase
Reads supporting (aligned to +/- strand):  ref base A (51/35);  new base G (3/2);  total (54/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 = 1.00e+00

AATATCCGCCGCCGTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGAT  >  NC_000913/3668037‑3668332
                                                                                                                                                    |                                                                                                                                                   
aaTATCCGCCGCCGTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAg                                                                                                                                                    >  2:904055/1‑150 (MQ=255)
aaTATCCGCCGCCGTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAg                                                                                                                                                    >  1:933010/1‑150 (MQ=255)
aaTATCCGCCGCCGTAATTTCCCCAAAGTTAAAGTCTTTCTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATGCCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAg                                                                                                                                                    <  2:1449923/150‑1 (MQ=255)
        ccgccgTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCTGACGGCGCaacag                                                                                                                                            <  1:634473/150‑5 (MQ=255)
        ccgccgTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCTGACGGCGCaacag                                                                                                                                            <  2:200734/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:1205315/1‑150 (MQ=255)
                    ccccAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCa                                                                                                                                <  1:1246513/150‑1 (MQ=255)
                    ccccAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCa                                                                                                                                <  2:11837/150‑1 (MQ=255)
                    ccccAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCa                                                                                                                                >  1:772059/1‑150 (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:126233/1‑150 (MQ=255)
                                   cTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATg                                                                                                                 >  2:786915/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                                                                                                               >  2:215051/1‑146 (MQ=255)
                                         tCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATggg                                                                                                               <  1:215051/146‑1 (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                                                                               >  2:771562/1‑150 (MQ=255)
                                                                     gccgAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTc                                                                               <  1:843979/150‑1 (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                                                   <  2:1335404/130‑1 (MQ=255)
                                                                                                                     taatTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCa                                                   >  1:1335404/1‑130 (MQ=255)
                                                                                                                     taatTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGAc                               <  2:1189425/150‑1 (MQ=255)
                                                                                                                      aatTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACt                              <  2:1038560/150‑1 (MQ=255)
                                                                                                                      aatTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACt                              <  2:894875/150‑1 (MQ=255)
                                                                                                                          gATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGa                          <  2:300218/150‑1 (MQ=255)
                                                                                                                          gATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACATTGCGGATAGACTTCGa                          >  2:1577381/1‑150 (MQ=255)
                                                                                                                             aGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGa                       <  1:44304/150‑1 (MQ=255)
                                                                                                                             aGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGa                       >  1:279526/1‑150 (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:729101/1‑150 (MQ=255)
                                                                                                                                        tGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTa            >  1:250230/1‑150 (MQ=255)
                                                                                                                                        tGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTa            <  1:91429/150‑1 (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)
                                                                                                                                                  tCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGAt  <  1:1534134/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AATATCCGCCGCCGTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGAT  >  NC_000913/3668037‑3668332

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