Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,396,962 T→G 6.2% R278R (CGA→CGC rng ← ribonuclease G

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,396,9620TG 6.2% 396.8 / 12.9 112R278R (CGA→CGCrngribonuclease G
Reads supporting (aligned to +/- strand):  ref base T (52/53);  new base G (3/4);  total (55/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 121 total aligned reads displayed.

CGCTTCAATATTGGTATTGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGAGTCAACGCGAATGCGG  >  NC_000913/3396815‑3397109
                                                                                                                                                   |                                                                                                                                                   
   ttCAATATTGGTATTGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGcgctcgct                                                                                                                                                 >  3:844052/1‑149 (MQ=255)
   ttCAATATTGGTATTGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGcgctcgct                                                                                                                                                 >  3:606147/1‑149 (MQ=255)
    tCAATATTGGTATTGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGcgctcgct                                                                                                                                                 <  2:271914/148‑1 (MQ=255)
    tCAATATTGGTATTGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGcgctcgct                                                                                                                                                 >  1:271914/1‑148 (MQ=255)
       atatTGGTATTGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGa                                                                                                                                              <  1:556023/148‑1 (MQ=255)
       atatTGGTATTGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGAGCGCTGGaa                                                                                                                                             <  3:291122/149‑2 (MQ=255)
        tatTGGTATTGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGAtt                                                                                                                                            <  4:723978/149‑1 (MQ=255)
        tatTGGTATTGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGAGCGCTGGaac                                                                                                                                            <  3:452226/149‑3 (MQ=255)
        tatTGGTATTGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGAGCGCTGGaa                                                                                                                                             <  1:203266/148‑2 (MQ=255)
         atTGGTATTGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGAttt                                                                                                                                           >  4:500336/1‑149 (MQ=255)
         atTGGTATTGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGAtt                                                                                                                                            <  2:169383/148‑1 (MQ=255)
              tATTGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGtt                                                                                                                                       <  2:439650/148‑1 (MQ=255)
               aTTGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGttt                                                                                                                                      >  1:141809/1‑148 (MQ=255)
               aTTGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGttt                                                                                                                                      >  1:83951/1‑148 (MQ=255)
                ttGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTc                                                                                                                                    <  3:442345/149‑1 (MQ=255)
                ttGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTc                                                                                                                                    >  3:661140/1‑149 (MQ=255)
                   aaaaTGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGAc                                                                                                                                 >  4:454156/1‑149 (MQ=255)
                     aaTGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACa                                                                                                                                >  1:234950/1‑148 (MQ=255)
                         gtgtCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCaa                                                                                                                            <  1:142195/148‑1 (MQ=255)
                           gtcgtcCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAaga                                                                                                                         <  4:661140/149‑1 (MQ=255)
                              gtcCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGAt                                                                                                                      >  4:863778/1‑149 (MQ=255)
                               tcCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTTGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGAt                                                                                                                      <  1:430455/148‑1 (MQ=255)
                               tcCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATc                                                                                                                     <  4:720366/149‑1 (MQ=255)
                                   gATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGaaa                                                                                                                 <  3:842954/149‑1 (MQ=255)
                                             aCCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTgg                                                                                                        <  1:893252/148‑1 (MQ=255)
                                                gACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGt                                                                                                     >  2:531489/1‑148 (MQ=255)
                                                 aCAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTcc                                                                                                   >  3:586377/1‑149 (MQ=255)
                                                    aaCGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCtgt                                                                                                <  3:139051/149‑1 (MQ=255)
                                                       gCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTa                                                                                             >  4:768489/1‑149 (MQ=255)
                                                           cGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGc                                                                                         <  4:25336/149‑1 (MQ=255)
                                                            ggTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGc                                                                                         >  2:184122/1‑148 (MQ=255)
                                                              tATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTc                                                                                       >  1:148985/1‑148 (MQ=255)
                                                               aTTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCa                                                                                     <  3:506320/149‑1 (MQ=255)
                                                                ttGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAg                                                                                    >  3:588089/1‑149 (MQ=255)
                                                                 tGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCAAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAgc                                                                                   >  4:197269/1‑149 (MQ=255)
                                                                    tGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAgctt                                                                                 >  2:910867/1‑148 (MQ=255)
                                                                      tCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAgcttgc                                                                               >  1:565431/1‑148 (MQ=255)
                                                                       ccACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAgcttgctt                                                                             <  3:320562/149‑1 (MQ=255)
                                                                        cACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGcttgcttg                                                                            >  3:683520/1‑149 (MQ=255)
                                                                            gtggtCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCa                                                                         >  1:752484/1‑148 (MQ=255)
                                                                                 cATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCgg                                                                   <  3:863778/149‑1 (MQ=255)
                                                                                  aTCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCgg                                                                   <  2:83951/148‑1 (MQ=255)
                                                                                   tCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCggg                                                                  >  1:755165/1‑148 (MQ=255)
                                                                                     gCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAAt                                                               <  4:532291/149‑1 (MQ=255)
                                                                                       ttCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAag                                                                                                                          >  3:738764/1‑88 (MQ=255)
                                                                                       ttCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGAt                                                                                                                      <  4:336437/92‑1 (MQ=255)
                                                                                       ttCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGAt                                                                                                                      >  3:336437/1‑92 (MQ=255)
                                                                                            tCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGcgctcgc                                                                                                                                                  >  4:871438/1‑59 (MQ=255)
                                                                                            tCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGcgctcgc                                                                                                                                                  <  3:871438/59‑1 (MQ=255)
                                                                                            tCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTcc                                                        >  3:703347/1‑149 (MQ=255)
                                                                                              tGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGa                                                      <  3:820252/149‑1 (MQ=255)
                                                                                                gTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAgg                                                    >  4:713067/1‑149 (MQ=255)
                                                                                                    aTAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGaa                                                <  4:586377/149‑1 (MQ=255)
                                                                                                    aTAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGa                                                 <  2:755165/148‑1 (MQ=255)
                                                                                                     tAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGaa                                                >  1:342887/1‑148 (MQ=255)
                                                                                                        tGAGATAACCACCGGATTTCAGTTCCACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAAGTGAACTCa                                            <  4:844052/149‑1 (MQ=255)
                                                                                                         gagaTAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCaa                                           <  3:258387/149‑1 (MQ=255)
                                                                                                         gagaTAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCa                                            <  2:141809/148‑1 (MQ=255)
                                                                                                          agaTAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCCCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCat                                           <  1:877797/148‑2 (MQ=255)
                                                                                                            aTAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCa                                            <  2:370270/145‑1 (MQ=255)
                                                                                                            aTAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCCAGGTGAACTCAAGt                                         >  1:607536/1‑148 (MQ=255)
                                                                                                              aaCCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTaa                                       <  1:910867/148‑1 (MQ=255)
                                                                                                                 caccGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCt                                   >  4:205041/1‑149 (MQ=255)
                                                                                                                 caccGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCt                                   >  4:132468/1‑149 (MQ=255)
                                                                                                                  accGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCtt                                  >  3:10116/1‑149 (MQ=255)
                                                                                                                  accGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCtt                                  >  3:353528/1‑149 (MQ=255)
                                                                                                                   ccGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACAGCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCtt                                  >  1:796883/1‑148 (MQ=255)
                                                                                                                    cGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCa                                            >  2:546595/1‑137 (MQ=255)
                                                                                                                    cGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCa                                            <  1:546595/137‑1 (MQ=255)
                                                                                                                       aTTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTaa                             <  4:112184/149‑1 (MQ=255)
                                                                                                                       aTTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTaa                             <  3:132468/149‑1 (MQ=255)
                                                                                                                       aTTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTa                              <  2:234950/148‑1 (MQ=255)
                                                                                                                          tCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGt                           <  2:148985/148‑1 (MQ=255)
                                                                                                                          tCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTc                          <  3:713067/149‑1 (MQ=255)
                                                                                                                            aGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGt                                                   >  3:237029/1‑122 (MQ=255)
                                                                                                                            aGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGt                                                   <  4:237029/122‑1 (MQ=255)
                                                                                                                             gTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAgg                       <  4:10116/149‑1 (MQ=255)
                                                                                                                              ttCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTaa                             >  4:46540/1‑142 (MQ=255)
                                                                                                                              ttCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTaa                             <  3:46540/142‑1 (MQ=255)
                                                                                                                              ttCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGc                      >  3:233951/1‑149 (MQ=255)
                                                                                                                              ttCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGc                      >  3:2479/1‑149 (MQ=255)
                                                                                                                               tCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTa                                                                                             <  1:230522/77‑1 (MQ=255)
                                                                                                                               tCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTa                                                                                             >  2:230522/1‑77 (MQ=255)
                                                                                                                               tCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCg                     >  3:210123/1‑149 (MQ=255)
                                                                                                                                  aCTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTg                   <  2:752484/148‑1 (MQ=255)
                                                                                                                                  aCTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGa                  <  3:226022/149‑1 (MQ=255)
                                                                                                                                  aCTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGa                  <  3:454156/149‑1 (MQ=255)
                                                                                                                                   cTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGa                  <  1:184122/148‑1 (MQ=255)
                                                                                                                                   cTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGAg                 >  4:424197/1‑149 (MQ=255)
                                                                                                                                    tttGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGAGt                >  4:624817/1‑149 (MQ=255)
                                                                                                                                    tttGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGAGt                <  4:703347/149‑1 (MQ=255)
                                                                                                                                       gaaatCCAGCGAGCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGa                  <  4:395954/140‑1 (MQ=255)
                                                                                                                                        cGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGAGTCAAc            <  3:424197/149‑1 (MQ=255)
                                                                                                                                          ttCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGAGTCAAcg           >  2:754347/1‑148 (MQ=255)
                                                                                                                                          atCCAGCGAGCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCCTGCTTGTCCTCTTGGGAATGTAGTCCGCGGTGAACTCAAGTAACGCTTCGTAAGACAGGCGTGAGTCAAccc          >  4:496718/2‑147 (MQ=255)
                                                                                                                                            ccAGCGAGCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGAGTCAAcgcg         >  2:410851/1‑148 (MQ=255)
                                                                                                                                             cAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGAGTCAACGCGa        >  1:358784/1‑148 (MQ=255)
                                                                                                                                             cAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGAGTCAACGCGa        <  2:876139/148‑1 (MQ=255)
                                                                                                                                                cgctcgctGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGAGTCAACGCGAATGc    <  4:2479/149‑1 (MQ=255)
                                                                                                                                                  ctcgctGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGAGTCAACGCGAATGCgg  <  4:233951/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CGCTTCAATATTGGTATTGAAAATGGTGTCGTCCAGATTGCGATGACCGACAAACGCTCCGGTATTGATGTCCACGGTGGTCATCGCTTCGGTCTGGTCGATAATGAGATAACCACCGGATTTCAGTTCTACTTTGCGTTCCAGCGCTCGCTGGATTTCGTTTTCGACATCAAAGAGATCGAAAATCGGCTGGCGTCCTGTGTAATGCTCCAGCTTGCTTGTCATCTCGGGAATGTACTCCGAGGTGAACTCAAGTAACGCTTCGTAAGTCAGGCGTGAGTCAACGCGAATGCGG  >  NC_000913/3396815‑3397109

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: