Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,376,114 T→A 6.2% L142Q (CTG→CAG)  PP_2083 → lipase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,376,1140TA 6.2% 243.8 / 6.0 80L142Q (CTG→CAG) PP_2083lipase
Reads supporting (aligned to +/- strand):  ref base T (39/36);  new base A (2/3);  total (41/39)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

CCGGAAATGCGAGGCCACGGCCTGTCGCCGCGTAACCATGACTGGAAACACAACAGCGTTGCCGCCTACGCCCGTGATGACCTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGTTCGGTACCCAGATCAGCCGTGTCTATTGGCCGTTGAAAGTGCCGCCGCTGGCTTGGGGCGCGAAGCTGGTGCT  >  NC_002947/2375966‑2376261
                                                                                                                                                    |                                                                                                                                                   
ccGGAAATGCGAGGCCACGGCCTGTCGCCGCGTAACCATGACTGGAAACACAACAGCGTTGCCGCCTACGCCCGTGATGACCTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTg                                                                                                                                                    >  2:715614/1‑150 (MQ=255)
 cGGAAATGCGAGGCCACGGCCTGTCGCCGCGTAACCATGACTGGAAACACAACAGCGTTGCCGCCTACGCCCGTGATGACCTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCAgg                                                                                                                                                   >  1:370166/1‑150 (MQ=255)
  ggAAATGCGAGGCCACGGCCTGTCGCCGCGTAACCATGACTGGAAACACAACAGCGTTGCCGCCTACGCCCGTGATGACCTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTggg                                                                                                                                                  >  1:1321543/1‑150 (MQ=255)
   gAAATGCGAGGCCACGGCCTGTCGCCGCGTAACCATGACTGGAAACACAACAGCGTTGCCGCCTACGCCCGTGATGACCTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGt                                                                                                                                                 <  1:715614/150‑1 (MQ=255)
          gAGGCCACGGCCTGTCGCCGCGTAACCATGACTGGAAACACAACAGCGTTGCCGCCTACGCCCGTGATGACCTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGcacca                                                                                                                                          >  2:178836/1‑150 (MQ=255)
             gCCACGGCCTGTCGCCGCGTAACCATGACTGGAAACACAACAGCGTTGCCGCCTACGCCCGTGATGACCTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGTGCCAccc                                                                                                                                       <  2:66122/150‑1 (MQ=255)
             gCCACGGCCTGTCGCCGCGTAACCATGACTGGAAACACAACAGCGTTGCCGCCTACGCCCGTGATGACCTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCAccc                                                                                                                                       >  2:209336/1‑150 (MQ=255)
                 cGGCCTGTCGCCGCGTAACCATGACTGGAAACACAACAGCGTTGCCGCCTACGCCCGTGATGACCTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCAGGGTGGTGCCACCCAggg                                                                                                                                   <  2:964058/150‑2 (MQ=255)
                    ccTGTCGCCGCGTAACCATGACTGGAAACACAACAGCGTTGCCGCCTACGCCCGTGATGACCTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGc                                                                                                                                >  2:505352/1‑150 (MQ=255)
                    ccTGTCGCCGCGTAACCATGACTGGAAACACAACAGCGTTGCCGCCTACGCCCGTGATGACCTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCAGGGTGGTGCCAcccagggaat                                                                                                                                <  2:370166/150‑8 (MQ=255)
                          gccgcGTAACCATGACTGGAAACACAACAGCGTTGCCGCCTACGCCCGTGATGACCTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTgg                                                                                                                                    <  2:323888/140‑1 (MQ=255)
                          gccgcGTAACCATGACTGGAAACACAACAGCGTTGCCGCCTACGCCCGTGATGACCTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCAGGGTGGTGCCAcccaggga                                                                                                                                  >  2:512224/1‑137 (MQ=255)
                          gccgcGTAACCATGACTGGAAACACAACAGCGTTGCCGCCTACGCCCGTGATGACCTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCAGGGTGGTGCCAcccaggga                                                                                                                                  <  1:512224/142‑6 (MQ=255)
                          gccgcGTAACCATGACTGGAAACACAACAGCGTTGCCGCCTACGCCCGTGATGACCTGCCGCTGATCAACGCCTTCGTGCGCGAGCAAACGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTgg                                                                                                                                    >  1:323888/1‑140 (MQ=255)
                           ccgcGTAACCATGACTGGAAACACAACAGCGTTGCCGCCTACGCCCGTGATGACCTGCCGCTGATCAACGCATTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTTGTCATTCCCTGGGTGGCACCACCCTGCTAGCTGCGCTg                                                                                                                         <  2:170406/150‑1 (MQ=255)
                             gcgTAACCATGACTGGAAACACAACAGCGTTGCCGCCTACGCCCGTGATGACCTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTggg                                                                                                                       <  2:1321543/150‑1 (MQ=255)
                               gTAACCATGACTGTAAACACAACAGCGTTGCCGCCTACGCCCGTGATGACCTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCg                                                                                                                     >  2:1307344/1‑150 (MQ=255)
                                     aTGACTGGAAACACAACAGCGTTGCCGCCTACGCCCGTGATGACCTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCAGCACCCTGGTAGCTGCGCTGGGCGGTGGtt                                                                                                               >  2:1247922/1‑150 (MQ=255)
                                     aTGACTGGAAACACAACAGCGTTGCCGCCTACGCCCGTGATGACCTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGtt                                                                                                               >  1:1369842/1‑150 (MQ=255)
                                         cTGGAAACACAACAGCGTTGCCGCCTACGCCCGTGATGACCTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTACGTCATTTCCTGGGTGGCTCCACCCTGTTAGCTGCGCTGGGCGGTGGTTTTCt                                                                                                           <  1:178836/150‑1 (MQ=255)
                                          tGGAAACACAACAGCGTTGCCGCCTACGCCCGTGATGACCTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCACTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTg                                                                                                          >  2:317485/1‑150 (MQ=255)
                                           ggAAACACAACAGCGTTGCCGCCTACGCCCGTGATGACCTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCTGTGGTTTTCTgg                                                                                                         <  1:642224/150‑1 (MQ=255)
                                               acacAACAGCGTTGCCGCCTACGCCCGTGATGACCTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGc                                                                                                     <  1:869901/150‑1 (MQ=255)
                                               acacAACAGCGTTGCCGCCTACGCCCGTGATGACCTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGc                                                                                                     >  1:1225213/1‑150 (MQ=255)
                                                cacaACAGCGTTGCCGCCTACGCCCGTGATGACCTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGcc                                                                                                    <  1:832482/150‑1 (MQ=255)
                                                  caacaGCGTTGCCGCCTACGCCCGTGATGACCTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGa                                                                                                  <  1:643200/150‑1 (MQ=255)
                                                   aacaGCGTTGCCGCCTACGCCCGTGATGACCTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAg                                                                                                 >  2:21020/1‑150 (MQ=255)
                                                     caGCGTTGCCGCCTACGCCCGTGATGACCTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCa                                                                                               >  2:776521/1‑150 (MQ=255)
                                                         gTTGCCGCCTACGCCCGTGATGACCTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACGGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTg                                                                                           <  1:949275/150‑1 (MQ=255)
                                                           tGCCGCCTACGCCCGTGATGACCTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGc                                                                                         <  1:181386/150‑1 (MQ=255)
                                                            gccgccTACGCCCGTGATGACCTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGtttt                                                                                                             >  1:188093/1‑129 (MQ=255)
                                                            gccgccTACGCCCGTGATGACCTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGtttt                                                                                                             <  2:188093/129‑1 (MQ=255)
                                                             ccgccTACGCCCGTGATGACCTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCa                                                                                       <  1:21020/150‑1 (MQ=255)
                                                             ccgccTACGCCCGTGATGACCTGCCGCTGATCAACGCCGTCTTGCTCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCa                                                                                       <  2:1225213/150‑1 (MQ=255)
                                                                   aCGCCCGTGATGACCTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTgg                                                                                 <  1:628925/150‑1 (MQ=255)
                                                                     gCCCGTGATGACCTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGgcg                                                                               <  1:505352/150‑1 (MQ=255)
                                                                       ccGTGATGACCTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCt                                                                             >  1:411528/1‑150 (MQ=255)
                                                                       ccGTGATGACCTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCt                                                                             >  1:450550/1‑150 (MQ=255)
                                                                       ccGTGATGACCTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCt                                                                             >  1:584159/1‑150 (MQ=255)
                                                                          tgatgaCCTGCCGCTGATCAACGACTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGCGTGGCCCCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGtt                                                                          <  1:6333/150‑1 (MQ=255)
                                                                                ccTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGtt                                                                          <  1:1345658/144‑1 (MQ=255)
                                                                                ccTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGtt                                                                          >  2:1345658/1‑144 (MQ=255)
                                                                                  tGCCGCTGATCAACGCCTTCGTGCGCGAGGAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGtt                                                                          <  1:74195/142‑1 (MQ=255)
                                                                                  tGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGtt                                                                          >  2:74195/1‑142 (MQ=255)
                                                                                       cTGATCAACGCCTTCGTGGGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGTTCGGTACCCAGATc                                                             >  1:1047484/1‑150 (MQ=255)
                                                                                            cAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGTTCGGTACCCa                                                                 >  1:695113/1‑141 (MQ=255)
                                                                                            cAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGTTCGGTACCCa                                                                 <  2:695113/141‑1 (MQ=255)
                                                                                             aaCGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGTTCGGTACCCAGATCAGCCgt                                                       >  2:1087112/1‑150 (MQ=255)
                                                                                             aaCGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGTTCGGTACCCAGATCAGCCgt                                                       >  1:826089/1‑150 (MQ=255)
                                                                                                gCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGTTCGGTACCCAGATCAGCCGTGTc                                                    >  2:800995/1‑150 (MQ=255)
                                                                                                  cTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCt                                                                                                           <  1:687876/93‑1 (MQ=255)
                                                                                                  cTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCt                                                                                                           >  2:687876/1‑93 (MQ=255)
                                                                                                     cGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGTTCGGTACCCAGATCAGCCGTGTCTATTg                                               <  2:935046/150‑1 (MQ=255)
                                                                                                          gcgAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGTTCGGTACCCAGATCAGCCGTGTCTATTGGCCGt                                          <  2:1030892/150‑1 (MQ=255)
                                                                                                           cgAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGTTCGGTACCCAGATCAGCCGTGTCTATTGGCCGtt                                         <  2:1369842/150‑1 (MQ=255)
                                                                                                                aaTCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGTTCGGTACCCAGATCAGCCGTGTCTATTGGCCGTTGAAAg                                    >  2:641097/1‑150 (MQ=255)
                                                                                                                 aTCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGTTCGGTACCCAGATCAGCCGTGTCTATTGGCCGTTGAAAGt                                   >  1:725897/1‑150 (MQ=255)
                                                                                                                   cGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGTTCGGTACCCAGATCAGCCGTGTCTATTGGCCGTTGAAAGTgc                                 <  1:559489/150‑1 (MQ=255)
                                                                                                                           gcgcCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGTTCGGTACCCAGATCAGCCGTGTCTATTGGCCGTTGAAAGTGCCGCCGCTg                         >  2:903791/1‑150 (MQ=255)
                                                                                                                           gcgcCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGTTCGGTACCCAGATCAGCCGTGTCTATTGGCCGTTGAAAGTGCCGCCGCTg                         <  1:209336/150‑1 (MQ=255)
                                                                                                                              ccgcACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGTTCGGTACCCAGATCAGCCGTGTCTATTGGCCGTTGAAAGTGCCGCCGCTGGCt                      <  1:641097/150‑1 (MQ=255)
                                                                                                                                  aCTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTg                                                                                           >  1:636405/1‑77 (MQ=255)
                                                                                                                                  aCTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTg                                                                                           <  2:636405/77‑1 (MQ=255)
                                                                                                                                  aCTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGTTCGGTACCCAGATCAGCCGTGTCTATTGGCCGTTGAAAGTGCCGCCGCTGGCTTggg                  >  1:1220994/1‑150 (MQ=255)
                                                                                                                                       gTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGTTCGGTACCCAGATCAGCCGTGTCTATTGGCCGTTGAAAGTGCCGCCGCTGGCTTGGGGcgcg             >  1:1038953/1‑150 (MQ=255)
                                                                                                                                        ttGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGTTCGGTACCCAGATCAGCCGTGTCTATTGGCCGTTGAAAGTGCCGCCGCTGGCTTg                    >  1:877083/1‑142 (MQ=255)
                                                                                                                                        ttGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGTTCGGTACCCAGATCAGCCGTGTCTATTGGCCGTTGAAAGTGCCGCCGCTGGCTTg                    <  2:877083/142‑1 (MQ=255)
                                                                                                                                          ggTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGTTCGGTACCCAGATCAGCCGTGTCTATTGGCCGTTGAAAGTGCCGCCGCTGGCTTGGGGCGCGAAg          >  2:382634/1‑150 (MQ=255)
                                                                                                                                            tCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGTTCGGTACCCAGATCAGCCGTGTCTATTGGCCGTTGAAAGTGCCGCCGCTGGCTTGGGGCGCGAAGCt        >  2:648756/1‑150 (MQ=255)
                                                                                                                                             cATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGTTCGGTACCCAGATCAGCCgtgt                                                     >  2:575947/1‑104 (MQ=255)
                                                                                                                                             cATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGTTCGGTACCCAGATCAGCCgtgt                                                     <  1:575947/104‑1 (MQ=255)
                                                                                                                                             cATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGTTCGGTACCCAGATCAGCCGTGTCTATTGGCCGTTGAAAGTGCCGCCGCTGGCTTGGGGCGCGAAGCTg       <  1:800995/150‑1 (MQ=255)
                                                                                                                                              aTTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGTTCGGTACCCAGATCAGCCGTGTCTATTGGCCGTTGAAAGTGCCGCCGCTGGCTTGGGGCGCGAAGCTgg      >  2:693061/1‑150 (MQ=255)
                                                                                                                                               ttCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGTTCGGTACCCAGATCAGCCGTGTCTATTGGCCGTTGAAAGTGCCGCCGCTGGCTTGGGGCGCGAAGCTGGt     <  2:826089/150‑1 (MQ=255)
                                                                                                                                               ttCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGTTCGGTACCCAGATCAGCCGTGTCTATTGGCCGTTGAAAGTGCCGCCGCTGGCTTGGGGCGCGAAGCTGGt     <  1:903791/150‑1 (MQ=255)
                                                                                                                                                 cccTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGTTCg                                                                        >  2:946789/1‑81 (MQ=255)
                                                                                                                                                 cccTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGTTCg                                                                        <  1:946789/81‑1 (MQ=255)
                                                                                                                                                 cccTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGTTCGGTACCCAGATCAGCCGTGTCTATTGGCCGTTGAAAGTGCCGCCGCTGGCTTGGGGCGCGAAGCTGGTGc   <  2:411528/150‑1 (MQ=255)
                                                                                                                                                  ccTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGTTCGGTACCCAGATCAGCCGTGTCTATTGGCCGTTGAAAGTGCCGCCGCTGGCTTGGGGCGCGAAGCTGGTGCt  >  1:818816/1‑150 (MQ=255)
                                                                                                                                                  ccTGGGAGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTCTTCCGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGTTCGGTACCCAGCTCAGCCGTGTCTATTGGCCGTTGAAAGTGCCGCCGCTGGCGTGGGGCGCGACGCTGGTGCt  >  2:522872/1‑150 (MQ=255)
                                                                                                                                                    tGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGTTCGGTAccc                                                                  >  2:1423290/1‑84 (MQ=255)
                                                                                                                                                    tGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGTTCGGTAccc                                                                  <  1:1423290/84‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCGGAAATGCGAGGCCACGGCCTGTCGCCGCGTAACCATGACTGGAAACACAACAGCGTTGCCGCCTACGCCCGTGATGACCTGCCGCTGATCAACGCCTTCGTGCGCGAGCAATCGGGCCAGGCGCCGCACTGGGTTGGTCATTCCCTGGGTGGCACCACCCTGGTAGCTGCGCTGGGCGGTGGTTTTCTGGCTGCCGAGCAGTTGGCCAGCGTGGCGCTGTTCGGTACCCAGATCAGCCGTGTCTATTGGCCGTTGAAAGTGCCGCCGCTGGCTTGGGGCGCGAAGCTGGTGCT  >  NC_002947/2375966‑2376261

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