Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 404,189 T→A 6.4% Q624L (CAG→CTG)  PP_0337 ← sensory box protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947404,1890TA 6.4% 239.7 / 7.1 78Q624L (CAG→CTG) PP_0337sensory box protein
Reads supporting (aligned to +/- strand):  ref base T (38/35);  new base A (2/3);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.95e-01

GCGCTTCCAGGCCCAGCAGGCGACCGCTGCGCAGGCACAGCTTGGGCTGGTAGAACACGTCCAGCTCGTTCTGGGTCAGGGCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTTCATCAGGCTGGCCGGGTCGACGCCGTGCAACGGTGCACAGGCCAGGCCTACCGAGGCGGTGACGTTGATCAGCTGGTTGTCGACGAACATCGGC  >  NC_002947/404040‑404338
                                                                                                                                                     |                                                                                                                                                      
tctcTTCCAGGCCCAGCAGGCGACCGCTGCGCAGGCACAGCTTGGGCTGGTAGAACACGTCCAGCTCGTTCTGGGTCAGGGCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCt                                                                                                                                                        <  2:678795/147‑1 (MQ=255)
    ttCCAGGCCCAGCAGGCGACCGCTGCGNAGGCACAGCTTGNGCTGGTNGAACACGTCCAGCTCGTTCTGGGTCAGGGCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTgaga                                                                                                                                                    >  1:678795/1‑147 (MQ=255)
      ccAGGCCCAGCAGGCGACCGCTGCGCAGGCACAGCTTGGGCTGGTAGAACACGTCCAGCTCGTTCTGGGTCAGGGCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCt                                                                                                                                                 >  1:168086/1‑150 (MQ=255)
        aGGCCCAGCAGGCGACCGCTGCGCAGGCACAGCTTGGGCTGGTAGAACACGTCCAGCTCGTTCTGGGTCAGGGCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTgg                                                                                                                                               >  2:778281/1‑150 (MQ=255)
        aGGCCCAGCAGGCGACCGCTGCGCAGGCACAGCTTGGGCTGGTAGAACACGTCCAGCTCGTTCTGGGTCAGGGCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTgg                                                                                                                                               >  1:573842/1‑150 (MQ=255)
         ggCCCAGCCGGAGACCGCTGCGCAGGCACAGCTTGGGCTGGTAGAACACGTCCAGCACGTTATTTCAGACGTGACGCAGCACCTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGgt                                                                                                                                              <  1:732003/150‑1 (MQ=255)
                  ggCGACCGCTGCGCAGGCACAGCTTGGGCTGGTAGAACACGTCCAGCTCGTTCTGGGTCAGGGCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGt                                                                                                                                     >  1:495129/1‑150 (MQ=255)
                   gCGACCGCTGCGCAGGCACAGCTTGGGCTGGTAGAACACGTCCAGCTCGTTCTGGGTCAGGGCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTAcc                                                                                                                                       >  2:649396/1‑147 (MQ=255)
                   gCGACCGCTGCGCAGGCACAGCTTGGGCTGGTAGAACACGTCCAGCTCGTTCTGGGTCAGGGCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTAcc                                                                                                                                       <  1:649396/147‑1 (MQ=255)
                   gCGACCGCTGCGCAGGCACAGCTTGGGCTGGTAGAACACGTCCAGCTCGTTCTGGGTCAGGGCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGtt                                                                                                                                    <  2:219683/150‑1 (MQ=255)
                   gCGACCGCTGCGCAGGCACAGCTTGGGCTGGTAGAACACGTCCAGCTCGTTCTGGGTCAGGGCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGtt                                                                                                                                    <  1:341963/150‑1 (MQ=255)
                   gCGACCGCTGCGCAGGCACAGCTTGGGCTGGTAGAACACGTCCAGCTCGTTCTGGGTCAGGGCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGtt                                                                                                                                    <  2:621642/150‑1 (MQ=255)
                      aCCGCTGCGCAGGCACAGCTTGGGCTGGTAGAACACGTCCAGCTCGTTCTGGGTCAGGGTACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGttggc                                                                                                                                 <  1:339123/150‑1 (MQ=255)
                        cGCTGCGCAGGCACAGCTTGGGCTGGTAGAACACGTCCAGCTCGTTCTGGGTCAGGGCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGttggcct                                                                                                                               <  2:495129/150‑1 (MQ=255)
                          cTGCGCAGGCACAGCTTGGGCTGGTAGAACACGTCCAGCTCGTTCTGGGTCAGGGCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGttggccttg                                                                                                                             >  2:1209444/1‑150 (MQ=255)
                            gcgcAGGCACAGCTTGGGCTGGTAGAACACGTCCAGCTCGTTCTGGGTCAGGGCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGttggccttggc                                                                                                                           <  1:778281/150‑1 (MQ=255)
                                  gCACAGCTTGGGCTGGTAGAACACGTCCAGCTCGTTCTGGGTCAGGGCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTtgtg                                                                                                                     >  1:116359/1‑150 (MQ=255)
                                  gCACAGCTTGGGCTGGTAGAACACGTCCAGCTCGTTCTGGGTCAGGGCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTtgtg                                                                                                                     <  1:536573/150‑1 (MQ=255)
                                    acaGCTTGGGCTGGTAGAACACGTCCAGCTCGTTCTGGGTCAGGGCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCATtatgaa                                                                                                                   >  2:461146/1‑148 (MQ=255)
                                      aGCTTGGGCTGGTAGAACACGTCCAGCTCGTTCTGGGTCAGGGCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGtgcaat                                                                                                                 >  1:1105195/1‑150 (MQ=255)
                                       gCTTGGGCTGGTAGAACACGTCCAGCTCGTTCTGGGTCAGGGCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGtgcaatg                                                                                                                >  1:1015496/1‑150 (MQ=255)
                                        cTTGGGCTGGTAGAACACGTCCAGCTCGTTCTGGGTCAGGGCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGtgcaatgc                                                                                                               >  2:193200/1‑150 (MQ=255)
                                            ggCTGGTAGAACACGTCCAGCTCGTTCTGGGTCAGGGCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAgg                                                                                                           >  1:773909/1‑150 (MQ=255)
                                               tGGTAGAACACGTCCAGCTCGTTCTGGGTCAGGGCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCg                                                                                                        <  2:1015496/150‑1 (MQ=255)
                                                    gAACACGTCCAGCTCGTTCTGGGTCAGGGCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCAtt                                                                                                   <  2:168086/150‑1 (MQ=255)
                                                     aaCACGTCCAGCTCGTTCTGGGTCAGGGCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCAttc                                                                                                  <  1:193200/150‑1 (MQ=255)
                                                      acacGTCCAGCTCGTTCTGGGTCAGGGCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCt                                                                                                                                                 <  2:572694/102‑1 (MQ=255)
                                                      acacGTCCAGCTCGTTCTGGGTCAGGGCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCt                                                                                                                                                 >  1:572694/1‑102 (MQ=255)
                                                        acGTCCAGCTCGTTCTGGGTCAGGGCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCAttcttc                                                                                               <  1:206530/150‑1 (MQ=255)
                                                          gTCCAGCTCGTTCTGGGTCAGGGCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTtcat                                                                                             <  2:773909/150‑1 (MQ=255)
                                                          gTCCAGCTCGTTCTGGGTCAGGGCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTtcat                                                                                             >  1:131817/1‑150 (MQ=255)
                                                           tCCAGCTCGTTCTGGGTCAGGGCACGGCGCAGGTTGGTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTtcatc                                                                                            >  1:425980/1‑150 (MQ=255)
                                                             cAGCTCGTTCTGGGTCAGGGCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTTCATCAg                                                                                          >  2:1047661/1‑150 (MQ=255)
                                                               gCTCGTTCTGGGTCAGGGCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTTCATCAGGc                                                                                        <  1:326126/150‑1 (MQ=255)
                                                                    ttCTGGGTCAGGGCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGttgg                                                                                                                                  >  1:196295/1‑103 (MQ=255)
                                                                    ttCTGGGTCAGGGCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGttgg                                                                                                                                  <  2:196295/103‑1 (MQ=255)
                                                                    ttCTGGGTCAGGGCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTTCATCAGGCTGGcc                                                                                   >  2:1112271/1‑150 (MQ=255)
                                                                        gggTCAGGGCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTTCATCAGGCTGGCCGGGt                                                                               <  2:1105195/150‑1 (MQ=255)
                                                                        gggTCAGGGCACGGCGCAGGTTGTTCTCGACGAACAGATTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTTCATCAGGCTGGCCGGGt                                                                               <  1:461146/150‑1 (MQ=255)
                                                                          gTCAGGGCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTTCATCAGGCTGGCCGGGTCg                                                                             >  1:787268/1‑150 (MQ=255)
                                                                              gggCACGGCGCAGGTTGTTCTCGACGAACAGCTGGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTTCATCAGGCTGGCCGGGTCGACGc                                                                         >  2:59349/1‑150 (MQ=255)
                                                                               ggCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTTCATCAGGCTGGCCGGGTc                                                                              <  2:882885/144‑1 (MQ=255)
                                                                               ggCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTTCATCAGGCTGGCCGGGTc                                                                              >  1:882885/1‑144 (MQ=255)
                                                                                  aCGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCg                                                                                                                                      <  2:886145/85‑1 (MQ=255)
                                                                                  aCGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCg                                                                                                                                      >  1:886145/1‑85 (MQ=255)
                                                                                  aCGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCACGGCCGGCCTTCTTCATCCGGCTCGCCGGCTCGACGCCGTg                                                                     >  2:737962/1‑150 (MQ=255)
                                                                                  aCGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTTCATCAGGCTGGCCGGGTCGACGCCGTg                                                                     <  1:523243/150‑1 (MQ=255)
                                                                                    ggCGCAGGTTGTTCGCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTTCATCAGGCTGGCCGGGTCGACGCCGTGCa                                                                   >  2:373803/1‑150 (MQ=255)
                                                                                       gcAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGGGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTTCATCAGGCTGGCCGGGTCGACGCCGTGCAACg                                                                <  2:425980/150‑1 (MQ=255)
                                                                                          ggTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTTCATCAGGCTGGCCGGGTCGACGCCGTGCAACGGTg                                                             <  1:373803/150‑1 (MQ=255)
                                                                                            ttgttCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTTCATCAGGCTGGCCGGGTCGACGCCGTGCAACAGCGca                                                           >  2:879496/1‑150 (MQ=255)
                                                                                                  tcGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGGCAGGGGCAATGCAAGGCCGGCATTCTTCATCAGGCTGGCCGGGTCGACGCCGTGCAACGGTGCAcaggcc                                                     >  1:720926/1‑150 (MQ=255)
                                                                                                    gacgaACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTTCATCAGGCTGGCCGGGTCGACGCCGTGCAACGGTGCAcaggccag                                                   >  1:1763/1‑150 (MQ=255)
                                                                                                    gacgaACAGCTTGTAGCTGGCTTCGGCGTTGAGCACCTCGGTAACCACCTGTACCTGGTGTTTACCGTTGGCCTTGGGCTTGTGCAATGCAAGGCCGGCATTCTTCATCAGGCTGGCCGGGTCGACGCCGTGCAACGGTGCAcaggccag                                                   <  1:1048946/150‑1 (MQ=255)
                                                                                                     acgaACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTTCATCAGGCTGGCCGGGTCGACGCCGTGCAACGGTGCAcaggccagg                                                  >  2:1069993/1‑150 (MQ=255)
                                                                                                       gaACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCATGTCTACTGCAAGGCCGGCATTCTTCATCAGTCTGGCCGGGTCCACGCCCGGCAACGGTTCACaggccaggcc                                                >  1:347650/1‑150 (MQ=255)
                                                                                                          cAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACCCCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTGTTCATCAGGCTGGCCGGGTCGACGCCGTGCAACGGTGCACTGGCCAGGCCTAc                                             >  1:1225912/1‑150 (MQ=255)
                                                                                                             cTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTtcatc                                                                                            >  1:691421/1‑100 (MQ=255)
                                                                                                             cTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTtcatc                                                                                            <  2:691421/100‑1 (MQ=255)
                                                                                                                 ttgcattgcacaaggccaaggccaaCGGTAAACACCAGGTACAGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTTCATCAGGCTGGCCGGGTCGACGCCGTGCAACGGTGCACAGGCCAGGCCTACCGAGGCg                                      <  2:228809‑M1/125‑1 (MQ=255)
                                                                                                                      ggCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTTCATCAGGCTGGCCGGGTCGACGCCGTGCAACGGTGCACAGGCCAGGCCTACCGAGGCGGTGAc                                 <  1:1112271/150‑1 (MQ=255)
                                                                                                                      ggCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTTCATCAGGCTGGCCGGGTCGACGCCGTGCAACGGTGCACAGGCCAGGCCTACCGAGGCGGTGAc                                 <  1:1069993/150‑1 (MQ=255)
                                                                                                                         ttCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTTCATCAGGCTGGCCGGGt                                                                               >  1:1084000/1‑101 (MQ=255)
                                                                                                                         ttCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTTCATCAGGCTGGCCGGGt                                                                               <  2:1084000/101‑1 (MQ=255)
                                                                                                                              cGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTTCATCAGGCTGGCCGGGTCGACGCCGTGCAACGGTGCACAGGCCAGGCCTACCGAGGCGGTGACGTTGATCa                         <  2:1002568/150‑1 (MQ=255)
                                                                                                                                 caaggccaaCGGTAAACACCAGGTACAGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCAttcttc                                                                                               >  2:415772‑M1/10‑77 (MQ=255)
                                                                                                                                 caaggccaaCGGTAAACACCAGGTACAGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCAttcttc                                                                                               <  1:415772‑M1/68‑1 (MQ=255)
                                                                                                                                  gAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTTCATCAGGCTGGCCGGGTCGACGCCGTGCAACGGTGCACAGGCCAGGCCTACCGAGGCGGTGACGTTGATCAGCTg                     >  2:1182422/1‑150 (MQ=255)
                                                                                                                                       caaCGGTAAACACCAGGTACAGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTTCATCAGGCTGGCCGGGTCGACGCCGTGCAACGGTGCACAGGCCAGGCCTACCGAGGCGGTGACGTTGATCAGCTGGTTGt                >  1:228809‑M1/4‑150 (MQ=255)
                                                                                                                                        ttCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTTCATCAGGCTGGCCGGGTCGACGCCGTGGAACGGTGCACAGGCCAGGCCTACCGAGGCGGTGACGTTGATCAGCTGGTTGTc               >  2:334064/1‑150 (MQ=255)
                                                                                                                                        aaCGGTAAACACCAGGTACAGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTTCATCAGGCTGGCCGGGTCGACGCCGTGCAACGCTGCACAGGCCAGGCCTAGCGAGGAGGTGGTGTTGATCGGCTGGTTgtt               <  2:522987‑M1/148‑2 (MQ=255)
                                                                                                                                          CGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTTCATCAGGCTGGCCGGGTCGACGCCGTGCAACGGTGCACAGGCCAGGCCTACCGAGGCGGTGACGTTGATCAGCTGGTTGTcga             >  2:720926‑M1/1‑150 (MQ=255)
                                                                                                                                           GGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTTCATCAGGCTGGCCGGGTCGACGCCGTGCAACGGTGCACAGGCCAGGCCTACCGAGGCGGTGACGTTGATCAGCTGGTTGTcgac            >  1:1047661‑M1/1‑150 (MQ=255)
                                                                                                                                             TAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTTCATCAGGCTGGCCGGGTCGACGCCGTGCAACGGTGCACAGGCCAGGCCTACCGAGGCGGTGACGTTGATCAgctggttgagatcgg          <  1:914391‑M1/150‑8 (MQ=255)
                                                                                                                                             TAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTTCATCAGGCTGGCCGGGTCGACGCCGTGCAACGGTGCACAGGCCAGGCCTACCGAGGCGGTGACGTTGATCAGCTGGTTGTcgacga          <  1:569204‑M1/150‑1 (MQ=255)
                                                                                                                                             TAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTTCATCAGGCTGGCCGGGTCGACGCCGTGCAACGGTGCACAGGCCAGGCCTAACGAGGCGGTGACGTTGATCAGCTGGTTGTcgacga          <  1:572278‑M1/150‑1 (MQ=255)
                                                                                                                                                  ACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTTCATCAGGCTg                                                                                      <  2:570588‑M1/69‑1 (MQ=255)
                                                                                                                                                  ACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTTCATCAGGCTg                                                                                      >  1:570588‑M1/1‑69 (MQ=255)
                                                                                                                                                  ACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTTCATCAGGCTGGCCGGGTCGACGCCGTGCAACGGTGCACAGGCCAGGCCTACCGAGGCGGTGACGTTGATCAGCTGGTTGTCGACGAACATc     <  2:702805‑M1/150‑1 (MQ=255)
                                                                                                                                                     tGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTTCATCAGGCTGGCCGGGTCGACGCCGTGCAACGGTGCACAGGCCAGGCCTACCGAGGCGGTGACGTTGATCAGCTGGTTGTCGACGAACATCGGc  <  2:1763/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GCGCTTCCAGGCCCAGCAGGCGACCGCTGCGCAGGCACAGCTTGGGCTGGTAGAACACGTCCAGCTCGTTCTGGGTCAGGGCACGGCGCAGGTTGTTCTCGACGAACAGCTTGTAGCTGGCTTCGGCGTTGAGCACTTCGGTAAACACCTGTACCTGGTGTTTACCGTTGGCCTTGGCCTTGTGCAATGCAAGGCCGGCATTCTTCATCAGGCTGGCCGGGTCGACGCCGTGCAACGGTGCACAGGCCAGGCCTACCGAGGCGGTGACGTTGATCAGCTGGTTGTCGACGAACATCGGC  >  NC_002947/404040‑404338

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