Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,238,702 T→C 5.1% A41A (GCT→GCC)  alx → putative membrane‑bound redox modulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,238,7020TC 5.1% 391.2 / 6.9 118A41A (GCT→GCC) alxputative membrane‑bound redox modulator
Reads supporting (aligned to +/- strand):  ref base T (48/64);  new base C (3/3);  total (51/67)
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.44e-01

Only 100 of 121 total aligned reads displayed.

CGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTCGCGCCGTTGCCGATCCACAGGCACTGGCCTTTCTCACAGGTTATCTGATTGAGAAATCGCTGGCG  >  NC_000913/3238554‑3238849
                                                                                                                                                    |                                                                                                                                                   
cGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCt                                                                                                                                                     >  4:474992/1‑149 (MQ=255)
cGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCt                                                                                                                                                     >  1:32065/1‑149 (MQ=255)
cGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGACGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCt                                                                                                                                                     >  3:568081/1‑149 (MQ=255)
 gTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTg                                                                                                                                                    >  2:802531/1‑149 (MQ=255)
   ttttttATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCg                                                                                                                                                  <  2:813945/149‑1 (MQ=255)
       ttATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCt                                                                                                                                              <  3:898958/149‑1 (MQ=255)
          tgtgTAAGGAACTTCTATGAATACTGTCGGCACGCCGTCGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCCGCAGCCTgtt                                                                                                                                           <  3:457733/149‑3 (MQ=255)
           gtgtAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTc                                                                                                                                          >  4:225561/1‑149 (MQ=255)
            tgtAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTTTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCCGCAGCCTgtttc                                                                                                                                         <  2:29330/149‑5 (MQ=255)
            tgtAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCg                                                                                                                                         >  1:6242/1‑149 (MQ=255)
             gtAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGc                                                                                                                                        <  2:304901/149‑1 (MQ=255)
                    aCTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTg                                                                                                                                 <  4:300379/149‑1 (MQ=255)
                        cTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTg                                                                                                                             >  1:806179/1‑149 (MQ=255)
                        cTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTg                                                                                                                             >  2:233170/1‑149 (MQ=255)
                          aTGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGAc                                                                                                                           <  1:362443/149‑1 (MQ=255)
                          aTGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCCGCAGCCTgtt                                                                                                                                           <  4:620785/133‑3 (MQ=255)
                          aTGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCCGCAGCCTgtt                                                                                                                                           >  3:620785/1‑131 (MQ=255)
                           tGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACg                                                                                                                          >  1:905389/1‑149 (MQ=255)
                              aTACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTg                                                                                                                       >  1:956646/1‑149 (MQ=255)
                                aCTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTc                                                                                                                     <  1:947062/149‑1 (MQ=255)
                                 cTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCg                                                                                                                    <  2:83298/149‑1 (MQ=255)
                                    tCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCATGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTgg                                                                                                                                >  2:222157/1‑134 (MQ=255)
                                    tCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTgg                                                                                                                                <  1:222157/134‑1 (MQ=255)
                                    tCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTa                                                                                                                 <  4:178993/149‑1 (MQ=255)
                                     cGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTAc                                                                                                                <  2:881845/149‑1 (MQ=255)
                                         aCGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGtt                                                                                                            <  1:216286/149‑1 (MQ=255)
                                             cGTTGCTAAGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTagc                                                                                                        <  4:505658/149‑3 (MQ=255)
                                                      gtggCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTc                                                                                               >  1:514278/3‑149 (MQ=255)
                                                      ggggCGGATTCGCTGTTGTTGTCGCCATTAGGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTc                                                                                               <  1:713757/149‑1 (MQ=255)
                                                        ggCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCtg                                                                                             <  1:448565/149‑1 (MQ=255)
                                                         gCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCtgg                                                                                            >  1:827533/1‑149 (MQ=255)
                                                           ggATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCtggtg                                                                                          <  1:506078/149‑1 (MQ=255)
                                                           ggATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCtggtg                                                                                          <  2:338656/149‑1 (MQ=255)
                                                              ttCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTa                                                                                       >  4:1002047/1‑149 (MQ=255)
                                                                  cTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTg                                                                                   <  2:956646/149‑1 (MQ=255)
                                                                   tgttgttgtCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTgg                                                                                  >  2:772576/1‑149 (MQ=255)
                                                                     ttgttgtCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTg                                                                                <  2:578892/149‑1 (MQ=255)
                                                                      tgttgtCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGc                                                                               <  3:225561/149‑1 (MQ=255)
                                                                         tgtCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCaaa                                                                            <  3:984685/149‑1 (MQ=255)
                                                                          gtCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAAc                                                                           <  1:1069727/149‑1 (MQ=255)
                                                                          gtCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAAc                                                                           <  2:1086880/149‑1 (MQ=255)
                                                                           tCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAAcc                                                                          >  3:478691/1‑149 (MQ=255)
                                                                           tCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAAcc                                                                          >  4:1021172/1‑149 (MQ=255)
                                                                             gCCATTATGCTGGCTGTCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGa                                                                        >  4:898380/1‑149 (MQ=255)
                                                                             gCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTc                                                                                                                                          >  3:666872/1‑83 (MQ=255)
                                                                             gCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTc                                                                                                                                          <  4:666872/83‑1 (MQ=255)
                                                                             gCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGa                                                                        >  2:705393/1‑149 (MQ=255)
                                                                             gCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGa                                                                        >  1:709703/1‑149 (MQ=255)
                                                                             gCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGa                                                                        >  1:755038/1‑149 (MQ=255)
                                                                              ccATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGaa                                                                       <  1:705393/149‑1 (MQ=255)
                                                                               cATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAg                                                                      <  3:896919/149‑1 (MQ=255)
                                                                                aTTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTgg                                                                                                                                <  2:486674/90‑1 (MQ=255)
                                                                                aTTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTgg                                                                                                                                >  1:486674/1‑90 (MQ=255)
                                                                                        ggCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTcg                                                                  >  4:329146/1‑144 (MQ=255)
                                                                                        ggCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTcg                                                                  <  3:727484/144‑1 (MQ=255)
                                                                                        ggCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTcg                                                                  <  3:329146/144‑1 (MQ=255)
                                                                                        ggCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTcg                                                                  >  4:727484/1‑144 (MQ=255)
                                                                                         gCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTCGCGCCGt                                                            >  3:338596/1‑149 (MQ=255)
                                                                                           tATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTAc                                                                                                                <  3:1069123/95‑1 (MQ=255)
                                                                                           tATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTAc                                                                                                                >  4:1069123/1‑95 (MQ=255)
                                                                                            aTCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTCGCGCCGTTGc                                                         <  2:6242/149‑1 (MQ=255)
                                                                                              cGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTCGCGCCGTTGCCg                                                       >  2:457550/1‑149 (MQ=255)
                                                                                               gACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTCGCGCCGTTGCCGa                                                      <  3:874556/149‑1 (MQ=255)
                                                                                               gACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTCGCGCCGTTGCCGa                                                      >  4:1051121/1‑149 (MQ=255)
                                                                                                 ccTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTCGCGCCGTTGCCGATc                                                    >  1:84784/1‑149 (MQ=255)
                                                                                                  cTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTCGCGCCGTTGCCGATcc                                                   >  4:895984/1‑149 (MQ=255)
                                                                                                   tgttgtTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTg                                                                                   <  1:944751/116‑1 (MQ=255)
                                                                                                   tgttgtTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTg                                                                                   >  2:944751/1‑116 (MQ=255)
                                                                                                      tgttgCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTCGCGCCGTTGCCGATCCACAg                                               <  2:642406/149‑1 (MQ=255)
                                                                                                      tgttgCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTCGCGCCGTTGCCGATCCACAg                                               <  4:332501/149‑1 (MQ=255)
                                                                                                        ttgCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTCGCGCCGTTGCCGATCCACAGGc                                             >  1:193570/1‑149 (MQ=255)
                                                                                                          gCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTCGGGCCGTTGCCGATCCACAGGCAc                                           >  3:663450/1‑149 (MQ=255)
                                                                                                             ggggCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTCGCGCCGTTGCCGATCCACAGGCACTgg                                        <  2:206365/149‑1 (MQ=255)
                                                                                                               ggCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTCGCGCCGTTGCCGATCCACAGGCACTGGcc                                      >  2:648019/1‑149 (MQ=255)
                                                                                                                 cgtcgtGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTCGCGCCGTTGCCGATCCACAGGCACTGGCCtt                                    <  3:98657/149‑1 (MQ=255)
                                                                                                                   tcgtGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTCGCGCCGTTGCCGATCCACAGGCACTGGCCTTtc                                  <  1:648019/149‑1 (MQ=255)
                                                                                                                    cgtGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTCGCGCCGTTGCCGATCCACAGGCACTGGCCTTtct                                 <  4:7195/149‑1 (MQ=255)
                                                                                                                            acaTGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTCGCGCCGTTGCCGATCCACAGGCACTGGCCTTTCTCACAGGtt                         <  3:343918/149‑1 (MQ=255)
                                                                                                                              aTGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTgg                                                                                  <  1:1076624/90‑1 (MQ=255)
                                                                                                                              aTGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTgg                                                                                  >  2:1076624/1‑90 (MQ=255)
                                                                                                                              aTGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTCGCGCCGTTGCCGATCCACAGGCACTGGCCTTTCTCACAGGTTAt                       >  3:559127/1‑149 (MQ=255)
                                                                                                                               tGCCATGAGCATGAGAGAGGCTGCGGCGTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTGAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTCGCGCCGTTGCCGATCCACAGGCACAGGCCTTTCTCACGGGTTATc                      <  2:970221/149‑1 (MQ=255)
                                                                                                                               tGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGGTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTCGCGCCGTTGCCGATCCACAGGCACTGGCCTTTCTCACAGGTTATc                      <  3:1021172/149‑1 (MQ=255)
                                                                                                                                   atgaccatgaAATAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTCGCGCCGTTGCCGATCCACAg                                               <  4:94100/120‑1 (MQ=255)
                                                                                                                                   atgaccatgaAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTCGCGCCGTTGCCGATCCACAg                                               >  3:94100/1‑120 (MQ=255)
                                                                                                                                    tgaccatgaAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTCGCGCCGTTGCCGATCCACAGGCACTGGCCTTTCTCACAGGTTATCTGAtt                 <  1:676262/149‑1 (MQ=255)
                                                                                                                                      gccatgaAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTCGCGCCGTTGCCGATCCACAGGCACTGGCCTTTCTCACAGGTTATCTGATTga               <  3:95762/148‑1 (MQ=255)
                                                                                                                                      accatgaAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGATGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTCGCGCCGTTGCCGATCCACAGGCACTGGCCTTTCTCACAGGTTATCTGATTga               <  2:499640/149‑1 (MQ=255)
                                                                                                                                        catgaAACAGTCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGCGCAAACCGAAGGTCGCGCCGTTGCCGATCCACAGGCACTGGCCTTTCTCACAGGTTATCTGATTgaga             >  3:69399/1‑149 (MQ=255)
                                                                                                                                         atgaAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTCGCGCCGTTGCCGATCCACAGGCACTGGCCTTTCTCACAGGTTATCTGa                   <  4:25506/142‑1 (MQ=255)
                                                                                                                                          tgaAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTCGCGCCGTTGCCGATCCACAGGCACTGGCCTTTCTCACAGGTTATCTGATTGAGaaa           >  3:108158/1‑149 (MQ=255)
                                                                                                                                            gaccAGGCCGCAGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTCGCGCCGTTGCCGATCCACAGGCACTGGCCTTTCTCACAGGTTATCTGATTGAGAAATc         >  4:894686/4‑149 (MQ=255)
                                                                                                                                            aaaCAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTCGCGCCGTTGCCGATCCACAGGCACTGGCCTTTCTCACAGGTTATCTGATTGAGAAATc         <  4:591123/149‑1 (MQ=255)
                                                                                                                                             accAGGCCGCAGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTCGCGCCGTTGCCGATCCACAGGCACTGGCCTTTCTCACAGGTTATCTGATTGAGAAATCg        >  1:315858/3‑149 (MQ=255)
                                                                                                                                             aaCAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTCGCGCCGTTGCCGATCCACAGGCACTGGCCTTTCTCACAGGTTATCTGATTGAGAAATCg        >  2:662303/1‑149 (MQ=255)
                                                                                                                                                aGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTCGCGCCGTTGCCGATCCACAGGCACTGGCCTTTCTCACAGGTTATCTGATTGAGAAATCGCTg     <  3:44482/149‑1 (MQ=255)
                                                                                                                                                aGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTCGCGCCGTTGCCGATCCACAGGCACTGGCCTTTCTCACAGGTTATCTGATTGAGAAATCGCTg     <  3:1059304/149‑1 (MQ=255)
                                                                                                                                                 ggCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTCGCGCCGTTGCCGATCCACAGGCACTGGCCTTTCTCACAGGTTATCTGATTGAGAAATCGCTgg    <  2:514278/149‑1 (MQ=255)
                                                                                                                                                   cTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTCGCGCCGTTGCCGATCCACAGGCACTGGCCTTTCTCACAGGTTATCTGATTGAGAAATCGCTGGCg  <  3:474992/149‑1 (MQ=255)
                                                                                                                                                   cTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTCGCGCCGTTGCCGATCCACAGGCACTGGCCTTTCTCACAGGTTATCTGATTGAGAAATCGCTGGCg  <  1:74282/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGCTATGGGGCGGATTCGCTGTTGTTGTCGCCATTATGCTGGCTATCGACCTGTTGTTGCAGGGGCGTCGTGGGGCACATGCCATGACCATGAAACAGGCTGCGGCCTGGTCGCTGGTCTGGGTGACGCTGTCGTTACTGTTTAACGCCGCTTTCTGGTGGTATCTGGTGCAAACCGAAGGTCGCGCCGTTGCCGATCCACAGGCACTGGCCTTTCTCACAGGTTATCTGATTGAGAAATCGCTGGCG  >  NC_000913/3238554‑3238849

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