Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,004,529 A→T 13.9% Q184L (CAG→CTG)  PP_0866 → transmembrane sensor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,004,5290AT13.9% 190.2 / 22.2 72Q184L (CAG→CTG) PP_0866transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base A (27/35);  new base T (5/5);  total (32/40)
Fisher's exact test for biased strand distribution p-value = 7.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.93e-01

CTCGCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCATATAGAT  >  NC_002947/1004389‑1004677
                                                                                                                                            |                                                                                                                                                    
cTCGCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTTAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCTGGGCAgtgt                                                                                                                                             <  2:921607/150‑1 (MQ=255)
cTCGCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTCATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAgtgt                                                                                                                                             >  1:216155/1‑150 (MQ=255)
   gCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGAt                                                                                                                                          <  2:651804/150‑1 (MQ=255)
        tGGAGCTCAATGGCAACACGCAGTTGTTGCACGTGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTca                                                                                                                                     >  1:749001/1‑150 (MQ=255)
           aGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCCGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAg                                                                                                                                  <  1:1017083/150‑1 (MQ=255)
           aGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAg                                                                                                                                  <  2:153013/150‑1 (MQ=255)
              tCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGt                                                                                                                                       >  1:83435/1‑142 (MQ=255)
              tCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGt                                                                                                                                       <  2:83435/142‑1 (MQ=255)
              tCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTa                                                                                                                               >  1:1269697/1‑150 (MQ=255)
               cAATGGCAACACGCTGATGCGGGGCGGGTTCTGCACCAGCCGGTGGCCCGTGTGCCTGAACGCGGGTGTAGTGATGTTCAGCGGCGCCCATGACAGTAACCGGCCGTCTGTGGTCGACACTGCCCAGGGCAGTGTTATCGTTACAGGTAc                                                                                                                              <  2:665732/150‑1 (MQ=255)
               cAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTAc                                                                                                                              <  1:896186/150‑1 (MQ=255)
                      aaCACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCttt                                                                                                                       <  2:855142/150‑1 (MQ=255)
                      aaCACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCttt                                                                                                                       >  2:427754/1‑150 (MQ=255)
                       acacGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTg                                                                                                                      >  2:775529/1‑150 (MQ=255)
                        cacGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGACCGATTGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGa                                                                                                                     <  1:374473/150‑1 (MQ=255)
                         acGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGAc                                                                                                                    >  2:834629/1‑150 (MQ=255)
                            cAGTTGCAGGTCGAGTTACGTGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTc                                                                                                                 <  1:783026/150‑1 (MQ=255)
                              gTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCTGGGCAgtg                                                                                                                                              >  1:1172399/1‑119 (MQ=255)
                              gTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCTGGGCAgtg                                                                                                                                              <  2:1172399/119‑1 (MQ=255)
                              gTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTcc                                                                                                                <  2:454539/149‑1 (MQ=255)
                              gTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTcc                                                                                                                >  1:454539/1‑149 (MQ=255)
                                      tCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGAcc                                                                                                       >  1:141845/1‑150 (MQ=255)
                                          gTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTCATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACccagc                                                                                                   <  2:216155/150‑1 (MQ=255)
                                              aGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCcagccagc                                                                                               <  1:431224/150‑1 (MQ=255)
                                              aGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGGGAAGTGATGTTCGGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCcagccagc                                                                                               >  2:561377/1‑150 (MQ=255)
                                                  ccGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCAccc                                                                                           <  1:1154995/150‑1 (MQ=255)
                                                       cAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGtgtg                                                                                      >  1:1192831/1‑150 (MQ=255)
                                                        aGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTgg                                                                                     >  1:119077/1‑150 (MQ=255)
                                                            gCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGt                                                                                 <  1:882775/150‑1 (MQ=255)
                                                            gCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGt                                                                                 <  2:1242649/150‑1 (MQ=255)
                                                                    gcCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACacg                                                                         <  2:1349401/150‑3 (MQ=255)
                                                                     ccTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAggg                                                                        <  2:1192831/150‑1 (MQ=255)
                                                                       tGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAggggc                                                                      <  1:1011025/150‑2 (MQ=255)
                                                                       tGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGt                                                                      >  1:406510/1‑150 (MQ=255)
                                                                           cGCTGGTGAAGGGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGt                                                                  <  2:141845/150‑1 (MQ=255)
                                                                             cTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCTCAAATCATTGGACTGTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGc                                                                <  2:749001/150‑1 (MQ=255)
                                                                               ggTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGGTCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCgg                                                              <  1:262669/150‑1 (MQ=255)
                                                                                      tGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAgg                                                       <  1:723748/150‑1 (MQ=255)
                                                                                        aTGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGc                                                     <  2:1269697/150‑1 (MQ=255)
                                                                                         tGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCtt                                                                                                                        >  1:755339/1‑82 (MQ=255)
                                                                                         tGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCtt                                                                                                                        <  2:755339/82‑1 (MQ=255)
                                                                                           ttCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAg                                                  <  2:406510/150‑1 (MQ=255)
                                                                                            tCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGATGCGGGTGCAGGGCAAgg                                                 <  2:119077/150‑1 (MQ=255)
                                                                                              aGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGt                                                                                 <  2:1058723/116‑1 (MQ=255)
                                                                                              aGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGt                                                                                 >  1:1058723/1‑116 (MQ=255)
                                                                                               gCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAgg                                                       >  1:292889/1‑141 (MQ=255)
                                                                                               gCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAgg                                                       <  2:292889/141‑1 (MQ=255)
                                                                                               gCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTATGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGGCGGTGCGGGTGCAGGGCAAGGGCg                                              >  2:702427/1‑150 (MQ=255)
                                                                                                      ccATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGcc                                       >  1:242193/1‑150 (MQ=255)
                                                                                                             aGTAACCGGCCGTTTGTGGTCGACACTGCCCATGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGAc                                            >  2:983564/1‑138 (MQ=255)
                                                                                                             aGTAACCGGCCGTTTGTGGTCGACACTGCCCATGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGAc                                            <  1:983564/138‑1 (MQ=255)
                                                                                                             aGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCgg                                                              <  2:346728/120‑1 (MQ=255)
                                                                                                             aGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCgg                                                              >  1:346728/1‑120 (MQ=255)
                                                                                                             aGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCaa                                <  1:427754/150‑1 (MQ=255)
                                                                                                                 aCCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGa                            >  1:723583/1‑150 (MQ=255)
                                                                                                                 aCCGGCCGTTTGTGATCGACACAGCACAAGGCAGTGTGATCTTCACAGGTACCCGCTGTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGa                            >  2:1378227/1‑150 (MQ=255)
                                                                                                                   cGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCgg                                                              >  1:1144497/1‑114 (MQ=255)
                                                                                                                   cGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCgg                                                              <  2:1144497/114‑1 (MQ=255)
                                                                                                                   cGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACg                          >  2:197880/1‑150 (MQ=255)
                                                                                                                    ggCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGAGCGGGTGCAGGGCAAGGGCGACTCGCTGGAGCACCTGACgg                         >  1:440099/1‑150 (MQ=255)
                                                                                                                        gTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGCGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCgg                     >  1:786043/1‑150 (MQ=255)
                                                                                                                            gtgGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAg                 <  1:775529/150‑1 (MQ=255)
                                                                                                                                ataacacTGCCCTGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGt             >  1:123580/4‑150 (MQ=255)
                                                                                                                                 tcacacTGCCCTGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCgg                                                              <  2:110657/98‑1 (MQ=255)
                                                                                                                                 tcacacTGCCCTGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCgg                                                              >  1:110657/3‑100 (MQ=255)
                                                                                                                                 tcacacTGCCCTGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAggg                                                >  1:552621/3‑114 (MQ=255)
                                                                                                                                 tcacacTGCCCTGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAggg                                                <  2:552621/112‑1 (MQ=255)
                                                                                                                                 tcacacTGCCCTGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTc            >  2:729897/3‑150 (MQ=255)
                                                                                                                                 cGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTc            <  1:896259/150‑1 (MQ=255)
                                                                                                                                      cTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCata       <  1:902489/150‑1 (MQ=255)
                                                                                                                                        gCCCTGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCAtata     <  2:123580/150‑1 (MQ=255)
                                                                                                                                           cAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCATATAGAt  >  2:51252/1‑150 (MQ=255)
                                                                                                                                            |                                                                                                                                                    
CTCGCATCTGGAGCTCAATGGCAACACGCAGTTGCAGGTCGAGTTCAGCGCCGGCCAGCGGCACGTGCGCCTGAACGCTGGTGAAGTGATGTTCAGCGTCGCCCATGACAGTAACCGGCCGTTTGTGGTCGACACTGCCCAGGGCAGTGTGATCGTCACAGGTACCCGCTTTGACGTCCGCCTGGACCCAGCCAGCACCCGTGTGGCGGTGGAACAGGGGTCGGTGCGGGTGCAGGGCAAGGGCGACTCGCCGGCGCAACTGACGGCCGGGCAGGGGTCGCATATAGAT  >  NC_002947/1004389‑1004677

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: ‑