Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,390,404 C→T 5.4% intergenic (‑44/+53) psd ← / ← rsgA phosphatidylserine decarboxylase/ribosome small subunit‑dependent GTPase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,390,4040CT 5.4% 247.3 / 3.3 74intergenic (‑44/+53)psd/rsgAphosphatidylserine decarboxylase/ribosome small subunit‑dependent GTPase A
Reads supporting (aligned to +/- strand):  ref base C (33/37);  new base T (2/2);  total (35/39)
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.30e-01

TCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCCGCGAT  >  NC_000913/4390261‑4390549
                                                                                                                                               |                                                                                                                                                  
tCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGcc                                                                                                                                               <  3:233469/149‑1 (MQ=255)
tCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGcc                                                                                                                                               >  3:369720/1‑149 (MQ=255)
  aGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCgg                                                                                                                                             >  4:98895/1‑149 (MQ=255)
  aGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCgg                                                                                                                                             >  3:48740/1‑149 (MQ=255)
    ccATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGAt                                                                                                                                           >  1:49476/1‑149 (MQ=255)
     cATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATc                                                                                                                                          <  4:48740/149‑1 (MQ=255)
     cATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGCAAGCCATAGTTTCGGCAGAACGTACTGTAGCGAAAGTTTAAGTGAATTTCACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCCGCCTACATGCCGGATc                                                                                                                                          <  3:430419/149‑1 (MQ=255)
          tGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCATGTAGGCCGGATc                                                                                                                                          <  2:403337/144‑1 (MQ=255)
          tGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCATGTAGGCCGGATc                                                                                                                                          >  1:403337/1‑144 (MQ=255)
             ccGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑aaaaa                                                                                                                                 >  2:169712/1‑149 (MQ=255)
               gcttgcttGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑aaaaaaa                                                                                                                               >  4:262973/1‑149 (MQ=255)
               gcttgcttGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑aaaaaaa                                                                                                                               <  2:310916/149‑1 (MQ=255)
                     ttgcGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGAc                                                                                                                         >  2:190116/1‑149 (MQ=255)
                     ttgcGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGAc                                                                                                                         >  2:362107/1‑149 (MQ=255)
                       gcgcCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCATAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGGCTCCAGGCCATTGTTTTGTCGTTCCTGGTCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGa                                                                                                                       <  3:580781/149‑1 (MQ=255)
                          ccccAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTc                                                                                                                    <  2:283585/149‑1 (MQ=255)
                          ccccAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTc                                                                                                                    >  1:73027/1‑149 (MQ=255)
                           cccAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGAtt                                                                                                                     <  4:292941/147‑1 (MQ=255)
                           cccAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGAtt                                                                                                                     >  3:292941/1‑147 (MQ=255)
                            ccAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCGGATCCGGCATGTAGGCCGGATcaggaa                                                                                                                                     <  2:377487/131‑6 (MQ=255)
                            ccAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACCGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCGGATCCGGCATGTAGGCCGGATcaggaa                                                                                                                                     >  1:377487/1‑126 (MQ=255)
                                    ccAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGCTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGt                                                                                                          <  3:674363/149‑1 (MQ=255)
                                      aGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGATTTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGCCGATTCTAACGACGGTTa                                                                                                        <  4:164536/149‑1 (MQ=255)
                                            gTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTATCCTCCAGGCCATTGTTTTGTCGTTGCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAac                                                                                                                >  3:36301/1‑135 (MQ=255)
                                            gTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAac                                                                                                                <  4:36301/135‑1 (MQ=255)
                                             tAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTAc                                                                                                                                                    >  1:153152/1‑99 (MQ=255)
                                             tAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTAc                                                                                                                                                    <  2:153152/99‑1 (MQ=255)
                                              aaGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAAtt                                                                                                <  3:411621/149‑1 (MQ=255)
                                              aaGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAAtt                                                                                                <  1:169712/149‑1 (MQ=255)
                                               aGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTg                                                                                               >  2:224539/1‑149 (MQ=255)
                                               aGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTg                                                                                               >  2:359823/1‑149 (MQ=255)
                                                  cATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTgtca                                                                                            <  3:262973/149‑1 (MQ=255)
                                                   aTAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTgtcag                                                                                           <  3:26070/149‑1 (MQ=255)
                                                        ttCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATc                                                                                      <  2:180595/149‑1 (MQ=255)
                                                        ttCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATc                                                                                      <  3:218200/149‑1 (MQ=255)
                                                         tCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATcc                                                                                     <  4:369720/149‑1 (MQ=255)
                                                             cAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTAt                                                                                 >  2:484213/1‑149 (MQ=255)
                                                             cAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTAt                                                                                 <  2:73027/149‑1 (MQ=255)
                                                                  tGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGaa                                                                            >  3:316207/1‑149 (MQ=255)
                                                                       tGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGtt                                                                       >  1:657306/1‑149 (MQ=255)
                                                                        gTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGttt                                                                      <  1:359823/149‑1 (MQ=255)
                                                                        gTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGttt                                                                      <  1:484213/149‑1 (MQ=255)
                                                                        gTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGttt                                                                      <  1:224539/149‑1 (MQ=255)
                                                                         tAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGttt                                                                      >  2:489500/1‑148 (MQ=255)
                                                                         tAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGttt                                                                      <  1:489500/148‑1 (MQ=255)
                                                                            cGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGttt                                                                      <  1:285293/145‑1 (MQ=255)
                                                                            cGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGttt                                                                      >  2:285293/1‑145 (MQ=255)
                                                                            cGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGtt                                                                       <  2:558379/144‑1 (MQ=255)
                                                                            cGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGtt                                                                       >  1:558379/1‑144 (MQ=255)
                                                                              aaaGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTAcgc                                                                >  2:645963/1‑149 (MQ=255)
                                                                                  tttAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAgg                                                                                                                             <  1:147325/84‑1 (MQ=255)
                                                                                  tttAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAgg                                                                                                                             >  2:147325/1‑84 (MQ=255)
                                                                                   ttAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGtttt                                                           <  1:5560/149‑1 (MQ=255)
                                                                                      aaTGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGaaaa                                                                          <  3:449856/131‑1 (MQ=255)
                                                                                      aaTGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAATCATCCGTATCAGaaaa                                                                          >  4:449856/1‑131 (MQ=255)
                                                                                       aTGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTAcc                                                       >  4:271617/1‑149 (MQ=255)
                                                                                          aaTTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTgc                                                    <  2:321154/149‑1 (MQ=255)
                                                                                            tttAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTgcgc                                                  >  1:68001/1‑149 (MQ=255)
                                                                                               aaCAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCAt                                               >  2:509231/1‑149 (MQ=255)
                                                                                                      tAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTcc                                        <  3:98895/149‑1 (MQ=255)
                                                                                                        gCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAg                                      >  3:112848/1‑149 (MQ=255)
                                                                                                             cAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATa                                  <  2:418307/148‑1 (MQ=255)
                                                                                                             cAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATa                                  >  1:418307/1‑148 (MQ=255)
                                                                                                                gCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGt                              <  4:466591/149‑1 (MQ=255)
                                                                                                                              gTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGaaa                >  3:72054/1‑149 (MQ=255)
                                                                                                                              gTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGaaa                <  3:271617/149‑1 (MQ=255)
                                                                                                                               ttCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATCGTTTTCGAAAc               <  1:362107/149‑1 (MQ=255)
                                                                                                                                 ccTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACgg             <  4:112848/149‑1 (MQ=255)
                                                                                                                                 ccTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACgg             <  1:509231/149‑1 (MQ=255)
                                                                                                                                   tGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGt           >  3:628070/1‑149 (MQ=255)
                                                                                                                                      tCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTc        >  1:509366/1‑149 (MQ=255)
                                                                                                                                       ccGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTcc       <  4:628070/149‑1 (MQ=255)
                                                                                                                                         ggCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCcgc     <  1:267070/149‑1 (MQ=255)
                                                                                                                                          gCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCcgcg    >  4:674688/1‑149 (MQ=255)
                                                                                                                                          gCCTACATGCCGGATCCCGAAAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCcgc     >  3:343893/1‑149 (MQ=255)
                                                                                                                                            cTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCCGCGAt  >  2:503151/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
TCAGCCATCCTGCCCGCTTGCTTGCGCCCCAACCCGCCAGGCGAGTAAGCCATAGTTTCGGCAGAATGTACTGTAGCGAAAGTTTAAATGAATTTAACAAGGTAGCCTCCAGGCCATTGTTTTGTCGTTCCTGATCCGGCCTACATGCCGGATCCTG‑AAAAAAAGGGGACGATTCTAACGACGGTTAGCTTAATTGTCAGTCATCCGTATCAGAAAAGTTTTTACGCGTTTTTACCTGCGCCATGCTTTCCAGAATACGGTGATAGTTTTCGAAACGGGTTTCCGCGAT  >  NC_000913/4390261‑4390549

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