Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 636,273 C→T 5.2% T99T (ACG→ACA ybdN ← PAPS reductase‑like domain protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913636,2730CT 5.2% 256.0 / 2.9 77T99T (ACG→ACAybdNPAPS reductase‑like domain protein
Reads supporting (aligned to +/- strand):  ref base C (31/42);  new base T (2/2);  total (33/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.42e-01

ACGTACAAATTGTTCGAAGGTCATGCCTGGCTGGTAAAAGCAGAAAAAGTCAGGGTCGGTTATCGCATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAGGGACTGAACATAGTTAATGGTGCAAGAGAATTGCGCCTCCCAGTCGATAAACAAAACGCAGATTTTTTTGCCCATTTGTC  >  NC_000913/636126‑636415
                                                                                                                                                   |                                                                                                                                              
aCGTACAAATTGTTCGAAGGTCATGCCTGGCTGGTAAAAGCAGAAAAAGTCAGGGTCGGTTATCGCATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGcg                                                                                                                                               >  3:474151/1‑149 (MQ=255)
  gTACAAATTGTTCGAAGGTCATGCCTGGCTGGTAAAAGCAGAAAAAGTCAGGGTCGGTTATCGCATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGcgtc                                                                                                                                             >  2:94895/1‑149 (MQ=255)
  gTACAAATTGTTCGAAGGTCATGCCTGGCTGGTAAAAGCAGAAAAAGTCAGGGTCGGTTATCGCATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGcgtc                                                                                                                                             >  3:168734/1‑149 (MQ=255)
   tACAAATTGTTCGAAGGTCATGCCTGGCTGGTAAAAGCAGAAAAAGTCAGGGTCGGTTATCGCATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGcgtcg                                                                                                                                            >  4:207181/1‑149 (MQ=255)
        aTTGTTCGAAGGTCATGCCTGGCTGGTAAAAGCAGAAAAAGTCAGGGTCGGTTATCGCATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGc                                                                                                                                       >  3:129092/1‑149 (MQ=255)
        aTTGTTCGAAGGTCATGCCTGGCTGGTAAAAGCAGAAAAAGTCAGGGTCGGTTATCGCATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGc                                                                                                                                       <  3:197922/149‑1 (MQ=255)
        aTTGTTCGAAGGTCATGCCTGGCTGGTAAAAGCAGAAAAAGTCAGGGTCGGTTATCGCATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGc                                                                                                                                       <  1:448988/149‑1 (MQ=255)
            ttCGAAGGTCATGCCTGGCTGGTAAAAGCAGAAAAAGTCAGGGTCGGTTATCGCATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGa                                                                                                                                   <  4:192799/149‑1 (MQ=255)
            ttCGAAGGTCATGCCTGGCTGGTAAAAGCAGAAAAAGTCAGGGTCGGTTATCGCATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGa                                                                                                                                   <  3:207181/149‑1 (MQ=255)
                   gTCATGCCTGGCTGGTAAAAGCAGAAAAAGTCAGGGTCGGTTATCGCATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCCTCGTAAGCGGTAGCGCAAc                                                                                                                            <  2:118094/149‑1 (MQ=255)
                           tggctggTAAAAGCAGAAAAAGTCAGGGTCGGTTATCGCATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTaaa                                                                                                                    >  3:96067/1‑149 (MQ=255)
                            ggctggTAAAAGCAGAAAAAGTCAGGGTCGGTTATCGCATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTaaaa                                                                                                                   <  4:470994/149‑1 (MQ=255)
                            ggctggTAAAAGCAGAAAAAGTCAGGGTCGGTTATCGCATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTaaaa                                                                                                                   <  2:456825/149‑1 (MQ=255)
                                ggTAAAAGCAGAAAAAGTCAGGGTCGGTTATCGCATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAActct                                                                                                               <  2:586302/149‑1 (MQ=255)
                                 gTAAAAGCAGAAAAAGTCAGGGTCGGTTAGCGCATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCtt                                                                                                              >  4:455708/1‑149 (MQ=255)
                                    aaaGCAGAAAAAGTCAGGGTCGGTTATCGCATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGa                                                                                                           <  1:540134/149‑1 (MQ=255)
                                     aaGCAGAAAAAGTCAGGGTCGGTTATCGCATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGAt                                                                                                          <  1:222819/149‑1 (MQ=255)
                                     aaGCAGAAAAAGTCAGGGTCGGTTATCGCATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGAt                                                                                                          <  2:100287/149‑1 (MQ=255)
                                     aaGCAGAAAAAGTCAGGGTCGGTTATCGCATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGAt                                                                                                          <  2:264328/149‑1 (MQ=255)
                                      aGCAGAAAAAGTCAGGGTCGGTTATCGCATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATg                                                                                                         >  4:440743/1‑149 (MQ=255)
                                         aGAAAAAGTCAGGGTCGGTTATCGCATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATg                                                                                                         <  1:183542/146‑1 (MQ=255)
                                         aGAAAAAGTCAGGGTCGGTTATCGCATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATg                                                                                                         >  2:183542/1‑146 (MQ=255)
                                          gAAAAAGTCAGGGTCGGTTATCGCATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACAt                                                                                                     >  2:183336/1‑149 (MQ=255)
                                           aaaaaGTCAGGGTCGGTTATCGCATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGTGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATc                                                                                                    >  4:164977/1‑149 (MQ=255)
                                                 tCAGGGTCGGTTATCGCATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTa                                                                                              >  4:98864/1‑149 (MQ=255)
                                                 tCAGGGTCGGTTATCGCATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTa                                                                                              >  1:586979/1‑149 (MQ=255)
                                                  cAGGGTCGGTTATCGCATCTTGGGGAGGTTGACGTACACATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTAc                                                                                             <  1:378571/149‑1 (MQ=255)
                                                   aGGGTCGGTTATCGCATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCTACCCAGTAAAACTCTTCGATGACATCGGTGTACa                                                                                            >  3:351317/1‑149 (MQ=255)
                                                    gggTCGGTTATCGCATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACaa                                                                                           <  2:465511/149‑1 (MQ=255)
                                                     ggTCGGTTATCGCATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAAc                                                                                          <  4:474151/149‑1 (MQ=255)
                                                      gTCGGTTATCGCATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACt                                                                                         >  2:17147/1‑149 (MQ=255)
                                                       tCGGTTATCGCATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTc                                                                                        >  3:576286/1‑149 (MQ=255)
                                                          gTTATCGCATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGTGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTcgcg                                                                                     >  3:151654/1‑149 (MQ=255)
                                                          gTTATCGCATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTcgcg                                                                                     >  1:141778/1‑149 (MQ=255)
                                                          gTTATCGCATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTcgcg                                                                                     >  1:260798/1‑149 (MQ=255)
                                                            tATCGCATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCa                                                                                   <  1:47053/149‑1 (MQ=255)
                                                              tCGCATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAgg                                                                                 >  2:477061/1‑149 (MQ=255)
                                                               cGCATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAggg                                                                                >  4:548392/1‑149 (MQ=255)
                                                                 cATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAGGGAc                                                                              <  2:159939/149‑1 (MQ=255)
                                                                 cATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAGGGAc                                                                              <  4:462500/149‑1 (MQ=255)
                                                                         ggAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAGGGACTGAACATa                                                                      <  1:477061/149‑1 (MQ=255)
                                                                           aGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAGGGACTGAACATAGt                                                                    <  2:586979/149‑1 (MQ=255)
                                                                           aGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAGGGACTGAACATAGt                                                                    <  1:145867/149‑1 (MQ=255)
                                                                             gTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGa                                                                                                           <  2:94605/108‑1 (MQ=255)
                                                                             gTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGa                                                                                                           >  1:94605/1‑108 (MQ=255)
                                                                              ttGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAGGGACTGAACATAGTTaa                                                                 <  1:230879/149‑1 (MQ=255)
                                                                              ttGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAGGGACTGAACATAGTTaa                                                                 <  4:168734/149‑1 (MQ=255)
                                                                              ttGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAGGGACTGAACATAGTTaa                                                                 <  4:192554/149‑1 (MQ=255)
                                                                                gACGTACCGATTCGACCTCAGGTTCCCAGCACTGCCATGCGGGTTTGTATTGTGATAGGGAATTTTGCGTCGTATGCGGGAGCGCAACCCAGTAAAACTAGTCGATGACATCGGTGTACAACGCGCGCACCGACTTAACATAGTTAATg                                                               <  2:49543/149‑1 (MQ=255)
                                                                                gACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAGGGACTGAACATAGTTAATg                                                               <  2:141778/149‑1 (MQ=255)
                                                                                gACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAGGGACTGAACATAGTTAATg                                                               <  4:576286/149‑1 (MQ=255)
                                                                                 aCGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAGGGACTGAACATAGTTAATgg                                                              >  4:58945/1‑149 (MQ=255)
                                                                                  cGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAGGGACTGAACATAGTTAATGGt                                                             <  2:209446/149‑1 (MQ=255)
                                                                                      cccATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAGGGACTGAACATAGTTAATGGTGc                                                           <  3:415810/147‑1 (MQ=255)
                                                                                      cccATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAGGGACTGAACATAGTTAATGGTGc                                                           >  4:415810/1‑147 (MQ=255)
                                                                                      cccATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAGGGACTGAACATAGTTAATGGTGCaa                                                         >  2:144721/1‑149 (MQ=255)
                                                                                      cccATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAGGGACTGAACATAGTTAATGGTGCaa                                                         >  3:380970/1‑149 (MQ=255)
                                                                                      cccATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAGGGACTGAACATAGTTAATGGTGCaa                                                         >  2:562963/1‑149 (MQ=255)
                                                                                            cGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAGGGACTGAACATAGTTAATGGTGCAAGAGAAt                                                   <  1:157186/149‑1 (MQ=255)
                                                                                             gACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAGGGACTGAACATAGTTAATGGTGCAAGAGAAtt                                                  <  1:461762/149‑1 (MQ=255)
                                                                                                     gTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAGGGACTGAACATAGTTAATGGTGCAAGAGAATTGCGCCTcc                                          <  2:212013/149‑1 (MQ=255)
                                                                                                     gTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAGGGACTGAACATAGTTAATGGTGCAAGAGAATTGCGCCTcc                                          <  4:397067/149‑1 (MQ=255)
                                                                                                     gTTCCCAGCACTGCCATTCGGGTTGGTATTGTCAAAGGGAACTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAGGGACTGAACATAGTTAATGGTGCAAGAGAATTGCGCCTcc                                          <  4:351317/149‑1 (MQ=255)
                                                                                                          cAGCACTGCCATTCGGGTTGGTATTGTGAAAGTGAATTTTGTGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCAGTGTACAACTCGCGCAGGGACTGAACATAGTTAATGGTGCAAGAGAATTGCGCCTCCCAGTc                                     <  4:151654/149‑1 (MQ=255)
                                                                                                          cAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAGGGACTGAACATAGTTAATGGTGCAAGAGAATTGCGCCTCCCAGTc                                     >  2:314882/1‑149 (MQ=255)
                                                                                                               cTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCATCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACATCTCGCGCAGGGACTGAACATAGTTAATGGTGCAAGAGAATTGCGCCTCCCAGTCGATaa                                >  4:345682/1‑149 (MQ=255)
                                                                                                                tGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAGGGACTGAACATAGTTAATGGTGCAAGAGAATTGCGCCTCCCAGTCGATaaa                               >  3:570779/1‑149 (MQ=255)
                                                                                                                tGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAGGGACTGAACATAGTTAATGGTGCAAGAGAATTGCGCCTCCCAGTCGATaaa                               <  1:17147/149‑1 (MQ=255)
                                                                                                                    aTTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAGGGACTGAACATAGTTAATGGTGCAAGAGAATTGCGCCTCCCAGTCGATAAACaaa                           >  3:138695/1‑149 (MQ=255)
                                                                                                                    aTTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAGGGACTGAACATAGTTAATGGTGCAAGAGAATTGCGCCTCCCAGTCGATAAACaaa                           >  3:221119/1‑149 (MQ=255)
                                                                                                                      tCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAGGGACTGAACATAGTTAATGGTGCAAGAGAATTGCGCCTCCCAGTCGATAAACAAAAc                         <  4:13344/149‑1 (MQ=255)
                                                                                                                        gggTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAGGGACTGAACATAGTTAATGGTGCAAGAGAATTGCGCCTCCCAGTCGATAAACAAAACGc                       <  3:393606/149‑1 (MQ=255)
                                                                                                                           tcggTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAGGGACTGAAc                                                                         >  4:286629/3‑96 (MQ=255)
                                                                                                                           tcggTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAGGGACTGAAc                                                                         <  3:286629/94‑1 (MQ=255)
                                                                                                                              gTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAGGGACTGAACATAGTTAATGGTGCAAGAGAATTGCGCCTCCCAGTCGATAAACAAAACGCAGAttt                 >  2:132762/1‑149 (MQ=255)
                                                                                                                               tATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAGGGACTGAACATAGTTAATGGTGCAAGAGAATTGCGCCTCCCAGTCGATAAACAAAACGCAGAtttt                >  2:595112/1‑149 (MQ=255)
                                                                                                                               gaTTGCGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGCTGACATCGGTGTACAACTCGCGCAGGGCCTGAACATAGTTAATGGTGCAAGAGAATTGCGCCTCCCAGTCGCTAACCAAAACGCAGAtttt                >  4:156195/2‑149 (MQ=255)
                                                                                                                                         gggAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAGGGACTGAACATAGTTAATGGTGCAAGAGAATTGCGCCTCCCAGTCGATAAACaaa                           >  2:340110/1‑128 (MQ=255)
                                                                                                                                         gggAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAGGGACTGAACATAGTTAATGGTGCAAGAGAATTGCGCCTCCCAGTCGATAAACaaa                           <  1:340110/128‑1 (MQ=255)
                                                                                                                                          ggAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAGGGACTGAACATAGTTAATGGTGCAAGAGAATTGCGCCTCCCAGTCGATAAACAAAACGCAGATTTTTTTGCCCAttt     <  3:98864/149‑1 (MQ=255)
                                                                                                                                             aTTTTGTGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAGGGACTGAACATAGTTAATGGTGCAAGAGAATTGCGCCTCCCAGTCGATAAACAAAACGCAGATTTTTTTGCCCATTTGTc  <  3:164977/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
ACGTACAAATTGTTCGAAGGTCATGCCTGGCTGGTAAAAGCAGAAAAAGTCAGGGTCGGTTATCGCATCTTGGGGAGGTTGACGTACCCATTCGACATCAGGTTCCCAGCACTGCCATTCGGGTTGGTATTGTGAAAGGGAATTTTGCGTCGTAAGCGGGAGCGCAACCCAGTAAAACTCTTCGATGACATCGGTGTACAACTCGCGCAGGGACTGAACATAGTTAATGGTGCAAGAGAATTGCGCCTCCCAGTCGATAAACAAAACGCAGATTTTTTTGCCCATTTGTC  >  NC_000913/636126‑636415

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