Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,178,879 C→A 5.5% V633V (GTG→GTT)  PP_1033 ← sulfatase domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,178,8790CA 5.5% 227.8 / 2.5 73V633V (GTG→GTT) PP_1033sulfatase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (37/32);  new base A (2/2);  total (39/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 = 5.64e-01

ATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATGAAGAATGTGCG  >  NC_002947/1178732‑1179025
                                                                                                                                                   |                                                                                                                                                  
aagaacCAGAAAGCCCGCAATCATGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCa                                                                                                                                                   >  2:33075/5‑149 (MQ=39)
aagaacCAGAAAGCCCGCAATCATGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCa                                                                                                                                                   >  2:27132/5‑149 (MQ=39)
    acCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCAcgcg                                                                                                                                               >  7:303850/1‑149 (MQ=255)
       aGAAAGCCCGCATGATTGCGGGCTTTCTGGTTGTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCCCAACGCGtg                                                                                                                                             >  1:59696/1‑147 (MQ=255)
        gAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTgg                                                                                                                                           <  2:78304/149‑1 (MQ=255)
           aGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTggggg                                                                                                                                        >  6:129854/1‑149 (MQ=255)
            gCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGt                                                                                                                                       >  3:53871/1‑149 (MQ=255)
               cGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCt                                                                                                                                    <  7:234031/149‑1 (MQ=255)
               cGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCCCAACGCGTGGAGGGTgac                                                                                                                                    <  7:296828/149‑3 (MQ=255)
                gCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCCCAACGCGTGGAGGGTgacg                                                                                                                                   <  2:59696/149‑4 (MQ=255)
                        gCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGc                                                                                                                            <  4:304420/148‑1 (MQ=255)
                              tttCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTAACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCAt                                                                                                                     >  1:352147/1‑149 (MQ=255)
                                     ttCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAg                                                                                                              >  2:13982/1‑149 (MQ=255)
                                      tCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAgg                                                                                                             >  6:346946/1‑149 (MQ=255)
                                       cTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGt                                                                                                            >  3:381312/1‑149 (MQ=255)
                                           ttcgttcAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTa                                                                                                        <  4:161773/149‑1 (MQ=255)
                                              gttcAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTAtt                                                                                                      >  3:303529/1‑148 (MQ=255)
                                               ttcAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCTCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCAACCTCCACTCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTa                                                                                                    <  7:349711/149‑1 (MQ=255)
                                               ttcAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTa                                                                                                    <  3:123761/149‑1 (MQ=255)
                                                  aGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCNTGCGGGTCATCCTCAAGGTTGTATTTAGTc                                                                                                 <  1:33075/149‑1 (MQ=255)
                                                  aGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTc                                                                                                 <  1:27132/149‑1 (MQ=255)
                                                   gATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCCCGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGt                                                                                                  >  4:341517/1‑147 (MQ=255)
                                                     ttGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAgg                                                                                              >  7:30258/1‑149 (MQ=255)
                                                        cTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGgtg                                                                                           >  4:217823/1‑149 (MQ=255)
                                                         ttATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCgg                                                                                                                          <  3:321741/117‑1 (MQ=255)
                                                         ttATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCgg                                                                                                                          >  4:321741/1‑117 (MQ=255)
                                                               ttctGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCCCGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTAAGACAGGGTGTCAGAAc                                                                                    >  8:132176/1‑149 (MQ=255)
                                                                    gCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGATCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGa                                                                                 <  7:137031/147‑1 (MQ=255)
                                                                        tCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTc                                                                             <  4:53871/147‑1 (MQ=255)
                                                                            gCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGtat                                                                         <  5:346946/147‑1 (MQ=255)
                                                                            gCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTatat                                                                       <  7:221557/149‑1 (MQ=255)
                                                                            gCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTatat                                                                       >  5:259112/1‑149 (MQ=255)
                                                                             cAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATg                                                                      >  1:257907/1‑149 (MQ=255)
                                                                                ttAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGGCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACTCGTTGGGGGTGGTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTTAATTCGTATATGTAt                                                                   <  7:132176/148‑1 (MQ=255)
                                                                                aTAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTAt                                                                   <  4:303529/149‑1 (MQ=255)
                                                                                aTAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTAt                                                                   <  5:129854/149‑1 (MQ=255)
                                                                                aTAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTAt                                                                   <  2:352147/149‑1 (MQ=255)
                                                                                aTAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCACCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTGCGATCTGAATTCGTATATGTAt                                                                   <  2:336563/149‑1 (MQ=255)
                                                                                    cgcTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTAATCCTCAACGTTGTATTTAGACAGGGTGTCCGAACTGAATTCGAATATGTAATTg                                                                >  6:207785/1‑148 (MQ=255)
                                                                                       tttCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACg                                                            >  1:252831/1‑149 (MQ=255)
                                                                                           tGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTg                                                        >  1:351828/1‑149 (MQ=255)
                                                                                              tAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTTTTTGTACGGGTGCTc                                                     >  2:374562/1‑149 (MQ=255)
                                                                                                ggCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTcc                                                    <  2:257907/148‑1 (MQ=255)
                                                                                                      ataGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAg                                              <  4:381312/148‑1 (MQ=255)
                                                                                                       taGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGc                                            <  8:303850/149‑1 (MQ=255)
                                                                                                         ggAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCt                                           <  1:13982/148‑1 (MQ=255)
                                                                                                          gAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGa                                         >  3:32187/1‑149 (MQ=255)
                                                                                                          gAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGa                                         <  8:30258/149‑1 (MQ=255)
                                                                                                           aaaGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGAc                                        <  3:20668/149‑1 (MQ=255)
                                                                                                                 gcCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTTCTCCACCAAGGCTGACGCCAGt                                  <  7:244667/149‑1 (MQ=255)
                                                                                                                  ccTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGt                                                                     <  5:321715/113‑1 (MQ=255)
                                                                                                                  ccTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGt                                                                     >  6:321715/1‑113 (MQ=255)
                                                                                                                    tGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGATTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACAAAGGCTGACGCCAGTTgg                               >  5:380394/1‑149 (MQ=255)
                                                                                                                         ctcACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCa                            >  3:273739/1‑147 (MQ=255)
                                                                                                                          tcACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCa                            >  7:307716/1‑146 (MQ=255)
                                                                                                                                cGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTg                   <  7:355353/149‑1 (MQ=255)
                                                                                                                                 gATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTg                   >  1:299762/1‑148 (MQ=255)
                                                                                                                                   tGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGAc                <  4:273739/149‑1 (MQ=255)
                                                                                                                                    gAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACa               >  1:166155/1‑149 (MQ=255)
                                                                                                                                        tCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCaa                                               >  4:37422/1‑113 (MQ=255)
                                                                                                                                        tCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCaa                                               <  3:37422/113‑1 (MQ=255)
                                                                                                                                        tCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATgaa           >  6:334807/1‑149 (MQ=255)
                                                                                                                                        aaagcACCCCCAATGCGTGGAGGGTGATTTTCCTGCGGCTCATCCACAAGGTTGGATTTAGTCAGGGTGTCCGAACTGAATTCGTAAATGTAGTTGTAAGGATCCTACACCAAGGCTGACGCCAGTTGGTCAAAGAAAGTGGACATga            >  4:17824/5‑148 (MQ=255)
                                                                                                                                          gTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGt                    <  7:168185/138‑1 (MQ=255)
                                                                                                                                          gTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGt                    >  8:168185/1‑138 (MQ=255)
                                                                                                                                          gTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATgaag          >  4:98485/1‑148 (MQ=255)
                                                                                                                                          gTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATgaag          >  3:338378/1‑148 (MQ=255)
                                                                                                                                           tCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGc                                      <  7:76886/119‑1 (MQ=255)
                                                                                                                                           tCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGc                                      >  8:76886/1‑119 (MQ=255)
                                                                                                                                           tCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATgaagaa        <  4:32187/149‑1 (MQ=255)
                                                                                                                                           tCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATgaagaa        >  3:169759/1‑149 (MQ=255)
                                                                                                                                             aCCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATGAAGAAtg      >  6:381080/1‑149 (MQ=255)
                                                                                                                                                 tCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATGAAGAATGTGCg  >  5:151659/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGCGGGTCATCCTCAAGGTTGTATTTAGTCAGGGTGTCCGAACTGAATTCGTATATGTATTTGTACGGGTGCTCCACCAAGGCTGACGCCAGTTGGTCAAAGAACGTGGACATGAAGAATGTGCG  >  NC_002947/1178732‑1179025

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