Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,980,924 G→T 6.1% P230Q (CCG→CAG)  otsA ← trehalose‑6‑phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,980,9240GT 6.1% 234.4 / 3.1 78P230Q (CCG→CAG) otsAtrehalose‑6‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base G (41/32);  new base T (3/2);  total (44/34)
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.80e-01

CGCTTCATAGGCGAGAAAACGCTCTGGCAAACCTTTGGAATAATCCAGCCGTTCGACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTATGGCTTTTTGCGCTACGTGTCGTGACGCGGGTCAGGTTAGAAAGACAATCCAGGAACGCCAGACGATCGTTTTCTGTCTGGAAACCCAGCAAA  >  NC_000913/1980776‑1981070
                                                                                                                                                    |                                                                                                                                                  
cGCTTCATAGGCGAGAAAACGCTCTGGCAAACCTTTGGAATAATCCAGCCGTTCGACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCg                                                                                                                                                    <  2:43095/149‑1 (MQ=255)
       tAGGCGAGAAAACGCTCTGGCAAACCTTTGGAATAATCCAGCCGTTCGACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCg                                                                                                                                                    >  6:368870/1‑142 (MQ=255)
       tAGGCGAGAAAACGCTCTGGCAAACCTTTGGAATAATCCAGCCGTTCGACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCg                                                                                                                                                    <  5:368870/142‑1 (MQ=255)
       tAGGCGAGAAAACGCTCTGGCAAACCTTTGGAATAATCCAGCCGTTCGACAGAAAAGATATTTAGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTACAAt                                                                                                                                             >  4:231002/1‑149 (MQ=255)
       tAGGCGAGAAAACGCTCTGGCAAACCTTTGGAATAATCCAGACGTTCAACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAATTTGCGCCAGTTTTGGCGGCAGTGGCACGTCAGCCTGTTTGGCTATTTCTTTCGGTTCAAt                                                                                                                                             >  2:401855/1‑149 (MQ=255)
           cGAGAAAACGCTCTGGCAAACCTTTGGAATAATCCAGCCGTTCGACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCaa                                                                                                                                              >  5:236972/1‑144 (MQ=255)
           cGAGAAAACGCTCTGGCAAACCTTTGGAATAATCCAGCCGTTCGACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCaa                                                                                                                                              <  6:236972/144‑1 (MQ=255)
            gagaAAACGCTCTGGCAAACCTTTGGAATAATCCAGCCGTTCGACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGa                                                                                                                                        >  3:270217/1‑149 (MQ=255)
             agaAAACGCTCTGGCAAACCTTTGGAATAATCCAGCCGTTCGACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGAt                                                                                                                                       >  6:13238/1‑149 (MQ=255)
               aaaaCGCTCTGGCAAACCTTTGGAATAATCCAGCCGTTCGACAGAAAAGATAGTTTGTACGTTTTTCAGTTCCGCTTTAAGGTGCGCCAGTTTTGGCGGCAGTGGCGCGGCAGCCTGGTTGGCTATTTCGTTCGGGTCATTGCCGAGCg                                                                                                                                     >  8:85820/1‑149 (MQ=255)
                 aaCGCTCTGGCAAACCTTTGGAATAATCCAGCCGTTCGACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCggg                                                                                                                                   >  1:357767/1‑149 (MQ=255)
                    gCTCTGGCAAACCTTTGGAATAATCCAGCCGTTCGACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCggg                                                                                                                                   >  1:301014/1‑146 (MQ=255)
                     ctctGGCAAACCTTTGGAATTATCCAGCCGTTCGACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGa                                                                                                                               <  3:240701/149‑1 (MQ=255)
                     ctctGGCAAACCTTTGGAATAATCCAGCCGTTCGACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCTGTTCAATGCCGATCGGGTAGa                                                                                                                               >  8:308395/1‑149 (MQ=255)
                     ctctGGCAAACCTTTGGAATAATCCAGCCGTTCGACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGa                                                                                                                               <  2:172676/149‑1 (MQ=255)
                      tctGGCAAACCTTTGGAATAATCCAGCCGTTCGACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTCTCGGTTCAATGCCGATCGGGTAGAc                                                                                                                              >  4:4490/1‑149 (MQ=255)
                          gCAAACCTTGGGAATAATCCAGCCGGTCGACAGAACAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTGTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTGCAATGCCGATCGGGTAGACtt                                                                                                                            <  3:340106/147‑1 (MQ=255)
                            aaaCCTTTGGAATAATCCAGCCGTTCGACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCggg                                                                                                                                   <  5:295346/138‑1 (MQ=255)
                            aaaCCTTTGGAATAATCCAGCCGTTCGACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCggg                                                                                                                                   >  6:295346/1‑138 (MQ=255)
                            aaaCCTTTGGAATAATCCAGCCGTTCGACAGAAAAGATATTTTGTACGTTTTTCAGTTACGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGt                                                                                                                        >  1:108686/1‑149 (MQ=255)
                             aaCCTTTGGAATAATCCAGCCGTTCGACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCTGTTCAATGCCGATCGGGTAGACTTCTGtt                                                                                                                       >  2:111463/1‑149 (MQ=255)
                             aaCCTTTGGAATAATCCAGCCGTTCGACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCTGTTCAATGCCGATCGGGTAGACTTCTGtt                                                                                                                       >  2:105626/1‑149 (MQ=255)
                              aCCGTTGGAATAATCCAGCCGTTCGACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTAGTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTc                                                                                                                      >  6:106050/1‑149 (MQ=255)
                                cTTTGGAATAATCGAGCCGTTCGACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGa                                                                                                                    <  5:238462/149‑1 (MQ=255)
                                cTTTGGAATAATCCAGCCGTTCGACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGGTATGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGAATTCTGTTCg                                                                                                                     >  1:184263/1‑148 (MQ=255)
                                 tttGGAATAATCCAGCCGTTCGACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCggg                                                                                                                                   >  5:342169/1‑133 (MQ=255)
                                 tttGGAATAATCCAGCCGTTCGACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCggg                                                                                                                                   <  6:342169/133‑1 (MQ=255)
                                    ggAATAATCCAGCCGTTCGACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACt                                                                                                                             >  5:179494/1‑136 (MQ=255)
                                    ggAATAATCCAGCCGTTCGACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACt                                                                                                                             <  6:179494/136‑1 (MQ=255)
                                      aataatCCAGCCGTTCGACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCt                                                                                                              >  6:236720/1‑149 (MQ=255)
                                       ataatCCAGCCGTTCGACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCtt                                                                                                             >  6:169205/1‑149 (MQ=255)
                                       ataatCCAGCCGTTCGACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCtt                                                                                                             <  1:330975/149‑1 (MQ=255)
                                       ataatCCAGCCGTTCGACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCtt                                                                                                             >  6:275552/1‑149 (MQ=255)
                                           tCCAGCCGTTCGACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCggg                                                                                                                                   <  2:357580/123‑1 (MQ=255)
                                           tCCAGCCGTTCGACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCggg                                                                                                                                   >  1:357580/1‑123 (MQ=255)
                                           tCCAGCCGTTCGACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGcc                                                                                                         >  1:44452/1‑149 (MQ=255)
                                                   ttCGACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTa                                                                                                                                 >  1:278760/1‑117 (MQ=255)
                                                   ttCGACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTa                                                                                                                                 <  2:278760/117‑1 (MQ=255)
                                                      gACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTAt                                                                                              >  7:26165/1‑149 (MQ=255)
                                                      gACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTAt                                                                                              <  5:130177/149‑1 (MQ=255)
                                                      gACAGAAAAGATATTTTGCACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTAt                                                                                              <  3:189571/149‑1 (MQ=255)
                                                       aCAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCATTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGACCCAGGCTGTATg                                                                                             >  5:56712/1‑149 (MQ=255)
                                                         aGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGc                                                                                                          >  1:379872/1‑134 (MQ=255)
                                                         aGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGc                                                                                                          <  2:379872/134‑1 (MQ=255)
                                                           aaaaGATATTTTGTACGTTTTTCAGTTCAGCTTTAAGTTTCGCCAGTTTTGGCGGAAGTGGCCAGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGAAGATCGGGTAGACTTCTGTTCGAAATGATTTGCCCCAGGCTGTATGGCtt                                                                                         >  5:129484/1‑149 (MQ=255)
                                                            aaaGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTATGGCttt                                                                                        <  4:148933/149‑1 (MQ=255)
                                                              aGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTATGGCttttt                                                                                      <  2:357767/149‑1 (MQ=255)
                                                               gATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTATGGCTTTTTg                                                                                     >  6:169217/1‑149 (MQ=255)
                                                                 tatTTTGTTCGTTTTTCAGTTGCGCTTTAAGTTGCGCCAGTTTTGGGGGCAGTGGGCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTATGGCTTTTTgc                                                                                    <  8:75419/148‑1 (MQ=255)
                                                                             ttttCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTATGGCTTTTTGCGCTACGTGTCGTg                                                                       >  1:19488/1‑149 (MQ=255)
                                                                                   gTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTATGGCTTTTTGCGCTACGTGGCGTGACGCgg                                                                 <  2:44452/149‑1 (MQ=255)
                                                                                                       aGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTATGGCTTTTTGCGCTACGTGTCGTGACGCGGGTCAGGTTAGAAAGACAATc                                             <  5:275552/149‑1 (MQ=255)
                                                                                                       aGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTATGGCTTTTTGCGCTACGTGTCGTGACGCGGGTCAGGTTAGAAAGACAATc                                             <  5:169205/149‑1 (MQ=255)
                                                                                                        gTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCTGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTATGGCTTTTTGCGCTACGTGTCGTGACGCGGGTCAGGTTAGAAAGACAATcc                                            <  1:111463/149‑1 (MQ=255)
                                                                                                        gTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCTGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTATGGCTTTTTGCGCTACGTGTCGTGACGCGGGTCAGGTTAGAAAGACAATcc                                            <  1:105626/149‑1 (MQ=255)
                                                                                                         ttttGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTATGGCTTTTTGCGCTACGTGTCGTGACGCGGGTCAGGTTAGAAAGACAATCCa                                           <  5:169217/149‑1 (MQ=255)
                                                                                                                ggcAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTATGGCTTTTTGCGCTACGTGTCGTGACGCGGGTCAGGTTAGAAAGACAATCCAGGAACGc                                    >  3:69387/1‑149 (MQ=255)
                                                                                                                       gCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTATGGCTTTTTGCGCTACGTGTCGTGACGCGGGTCTGGTTAGAAAGACAATCCAGGAACGCCAGACGa                             <  4:7201/149‑1 (MQ=255)
                                                                                                                         ccGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTATGGCTTTTTGCGCTACGTGTCGTGACGCGGGTCAgg                                                           >  5:54335/1‑117 (MQ=255)
                                                                                                                         ccGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTATGGCTTTTTGCGCTACGTGTCGTGACGCGGGTCAgg                                                           <  6:54335/117‑1 (MQ=255)
                                                                                                                         ccGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTAGGGCTTTTTGCGCTACGTGTCGTGACGCGGGTCAGGTTAGAAAGACAATCCAGGAACGCCAGACGATc                           <  2:283653/149‑1 (MQ=255)
                                                                                                                         ccGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGATGTATGGCTTTTTGCGCTACGTGTCGTGACGCGGGTCAGGTTAGAAAGACAATCCAGGAACGCCAGACGATc                           >  8:286300/1‑149 (MQ=255)
                                                                                                                          cGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTATGGCTTTTTGCGCTACGTGTCGTGACGCGGGTCAGGTTAGAAAGACAATCCAGGAACGCCAGACGATCg                          <  4:69387/149‑1 (MQ=255)
                                                                                                                          cGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTATGGCTTTTTGCGCTACGTGTCGTGACGCGGGTCAGGTTAGAAAGACAATCCAGGAACGCCAGACGATCg                          >  1:185201/1‑149 (MQ=255)
                                                                                                                               gCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTATGGCTTTTTGCGCTACGTGTCGTGACGCGGGTCAGGTTAGAAAGACAATCCAGGAACGCCAGACGATCGTTTtc                     <  3:231002/149‑1 (MQ=255)
                                                                                                                                    tttGGCTATTTCTTTCGGTTCAATGCTGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTATGGCTTTTTGCGCTACGTGTCGTGACGCGGGTCAGGTTAGAAAGACAATCCAGGAACGCCAGACGATCGTTTtctgtc                 <  1:401855/148‑1 (MQ=255)
                                                                                                                                        gCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTATGGCTTTTTGCGCTACGTGTCGTGACGCGGGTCAGGTTAGAAAGACAATCCAGGAACGCCAGACGATCGTTTTCTGTCTGGa             >  6:359423/1‑148 (MQ=255)
                                                                                                                                         cTATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTATGGCTTTTTGCGCTACGTGTCGTGACGCGGGTCAGGTTAGAAAGACAATCCAGGAACGCCAGACGATCGTTTTCTGTCTGGaaa           <  7:211996/149‑1 (MQ=255)
                                                                                                                                          tATTTCTTTCGGTTCAATGCCGATCTGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTATGGCTTTTTGCGCTACGTGTCGTGACGCGGGTCAGGTTAGAAAGACAATCCAGGAACGCCAGACGATCGTTTTCTGTCTGGAAAc          >  1:307324/1‑149 (MQ=255)
                                                                                                                                          tATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTATGGCTTTTTGCGCTACGTGTCGTGACGCGGGTCAGGTTAGAAAGACAATCCAGGAACGCCAGACGATCGTTTTCTGTCTGGAAAc          >  7:309982/1‑149 (MQ=255)
                                                                                                                                          tATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTATGGCTTTTTGCGCTACGTGTCGTGACGCGGGTCAGGTTAGAAAGACAATCCAGGAACGCCAGACGATCGTTTTCTGTCTGGAAAc          >  7:315520/1‑149 (MQ=255)
                                                                                                                                          tATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATCCTTTGCCCCAGGCTGTATGGCTTTTTGCGCTACGTGTCGTGACGCGGGACAGGTTAGAAAGACAATCCAGGAACGCCAGACGATCGTTTTCTGTCTGGAAAc          >  7:128259/1‑149 (MQ=255)
                                                                                                                                            tttctttcGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTATGGCTTTTTGCGCTACGTGTCGTGACGCGGGTCAGGt                                                          >  2:53715/1‑99 (MQ=255)
                                                                                                                                            tttctttcGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTATGGCTTTTTGCGCTACGTGTCGTGACGCGGGTCAGGt                                                          <  1:53715/99‑1 (MQ=255)
                                                                                                                                            tttctttcGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTATGGCTTTTTGCGCTACGTGTCGTGACGCGGGTCAGGTTAGa                                                      <  5:382588/103‑1 (MQ=255)
                                                                                                                                            tttctttcGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTATGGCTTTTTGCGCTACGTGTCGTGACGCGGGTCAGGTTAGa                                                      >  6:382588/1‑103 (MQ=255)
                                                                                                                                            tttctttcGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTATGGCTTTTTGCGCTACGTGTCGTGACGCGGGTCAGGTTAGAAAGACAATCCAGGAACGCCAGACGATCGTTTTCTGTCTGGAAAcc         >  7:64590/1‑148 (MQ=255)
                                                                                                                                             ttctttcGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTATGGCTTTTTGCGCTACGTGTCGTGACGCGGGTCAGGTTAGAAAGACAATCCAGGAACGCCAGACGATCGTTTTCTGTCTGGAAAcc         >  7:114078/1‑147 (MQ=255)
                                                                                                                                              tctttcGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTATGGCTTTTTGCGCTACGTGTCGTGACGCGGGTCAGGTTAGAAAGACAATCCAGGAACGCCAGACGATCGTTTTCTGTCTGGAAACCCa       <  8:44585/148‑1 (MQ=255)
                                                                                                                                                tttcGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTATGGCTTTTTGCGCTACGTGTCGTGACGCGGGTCAGGTTAGAAAGACAATCCAGGAACGCCAGACGATCGTTTTCTGTCTGGAAACCCAGc     >  4:1533/1‑148 (MQ=255)
                                                                                                                                                   cGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGATGTATGGCTTTTTGCGCTACGTGTCGTGACGCGGGTCAGGTTAGAAAGACAATCCAGGAACGCCAGACGATCGTTTTCTGTCTGGAAACCCAGCaaa  <  7:286300/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CGCTTCATAGGCGAGAAAACGCTCTGGCAAACCTTTGGAATAATCCAGCCGTTCGACAGAAAAGATATTTTGTACGTTTTTCAGTTCCGCTTTAAGTTGCGCCAGTTTTGGCGGCAGTGGCCCGGCAGCCTGTTTGGCTATTTCTTTCGGTTCAATGCCGATCGGGTAGACTTCTGTTCGAAATGCTTTGCCCCAGGCTGTATGGCTTTTTGCGCTACGTGTCGTGACGCGGGTCAGGTTAGAAAGACAATCCAGGAACGCCAGACGATCGTTTTCTGTCTGGAAACCCAGCAAA  >  NC_000913/1980776‑1981070

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