Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,256,787 A→T 6.6% L55Q (CTG→CAG)  PP_3729 ← amino acid‑binding protein‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,256,7870AT 6.6% 228.0 / 7.0 77L55Q (CTG→CAG) PP_3729amino acid‑binding protein‑like protein
Reads supporting (aligned to +/- strand):  ref base A (33/38);  new base T (2/3);  total (35/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.59e-01

CCGGCTCCGTGGAAAAACGCATACGCTGGTCACGTTGGTCGCTGTTGTAGGCCAGGACCACGTCCACCCGCCCGGCGGCGGCATCCATCAGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTACACCGGGTAGTCGCAGTAGCGCACCGGCTCACTGGCACTGGCGTGCAGGGTGGGCAACGTCAAGGCCAGCGCCAACTGCCAG  >  NC_002947/4256638‑4256936
                                                                                                                                                     |                                                                                                                                                     
ccGGCTCCGTGGAAAAACGCATACGCTGGTCACGTTGGTCGCTGTTGTGATGTTGGACCACGTCCACCCGCCCGGCGGCGGCATCCATCAGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACa                                                                                                                                                       >  1:706335/1‑150 (MQ=255)
ccGGCTCCGTGGAAAAACGCATACGCTGGTCACGTTGGTCGCTGTTGTAGGCCAGGACCACGTCCACCCGCCCGGCGGCGGCATCCATCAGGCAGCGCTTCCAGTTGCCCCAAGGCGTGTTTTGCACCTCGTAGCCCATGCGCGCGAACa                                                                                                                                                       <  1:1207254/150‑1 (MQ=255)
     tCCGTGGAAAAACGCATACGCTGGTCACGTTGGTCGCTGTTGTAGGCCAGGACCACGTCCACCCGCCCGGCGGCGGCATCCATCAGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTcg                                                                                                                                                  >  2:1071005/1‑150 (MQ=255)
          ggAAAAACGCATACGCTGGTCACGTTGGTCGCTGTTGTAGGCCAGGACCACGTCCACCCGCCCGGCGGCGGCATCCATCAGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGcac                                                                                                                                             >  2:990004/1‑150 (MQ=255)
           gAAAAACGCATACGCTGGTCACGTTGGTCGCTGTTGTAGGCCAGGACCACGTCCACCCGCCCGGCGGCGGCATCCATCAGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACTGTTCgcgc                                                                                                                                               <  2:512624/147‑1 (MQ=255)
           gAAAAACGCATACGCTGGTCACGTTGGTCGCTGTTGTAGGCCAGGACCACGTCCACCCGCCCGGCGGCGGCATCCATCAGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACTGTTCgcgc                                                                                                                                               >  1:512624/1‑147 (MQ=255)
             aaaaCGCATACGCTGGTCACGTTGGTCGCTGTTGTAGGCCAGGACCACGTCCACCCGCCCGGCGGCGGCATCCATCAGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACTGTTCgcgcgcatg                                                                                                                                          <  2:1187938/150‑6 (MQ=255)
                aCGCATACGTTGGTTATTTTTGTCGCTGTTGTAGGCCAGGACCACGTCCACCCGCCCGGCGGCAGCATCCATCAGGTAGCGCTTACAGTTGCCCAGCACCACGGTGTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGt                                                                                                                                       <  2:706335/150‑1 (MQ=255)
                aCGCATACGCTGGTCACGTTGGTCGCTGTTGTAGGCCAGGACCACGTCCACCCGCCCGGCGGCGGCATCCATCAGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGt                                                                                                                                       >  2:1202168/1‑150 (MQ=255)
                  gCATACGCTGGTCACGTTGGTCGCTGTTGTAGGCCAGGACCACGTCCACCCGCCCGGCGGCGGCATCCATCAGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACTGTTCgcgcgca                                                                                                                                            >  2:225546/1‑140 (MQ=255)
                  gCATACGCTGGTCACGTTGGTCGCTGTTGTAGGCCAGGACCACGTCCACCCGCCCGGCGGCGGCATCCATCAGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACTGTTCgcgcgca                                                                                                                                            <  1:225546/143‑4 (MQ=255)
                   cATACGCTGGTCACGTTGGTCGCTGTTGTAGGCCAGGACCACGTCCACCCGCCCGGCGGCGGCATCCATCAGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTggg                                                                                                                                    <  2:1112976/150‑1 (MQ=255)
                   cATACGCTGGTCACGTTGGTCGCTGTTGTAGGCCAGGACCACGTCCACCCGCCCGGCGGAGACATCCGTCAGGCAGCGCGTCCAGTTGCCCAGCACCACGGGTTGCACCTCGTAGCCCATGCGCGAGAGCACAGGGCGCACCACGGTggg                                                                                                                                    >  1:498924/1‑150 (MQ=255)
                    aTACGCTGGTCACGTTGGTCGCTGTTGTAGGCCAGGACCACGTCCCCCCGCCCGGCGGCGGCATCCATCAGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTGGCCCATGCGCGCGAACAGTTCGCGCACCACGGTgggg                                                                                                                                   <  2:1349210/150‑1 (MQ=255)
                       cGCTGGTCACGTTGGTCGCTGTTGTAGGCCAGGACCACGTCCACCCGCCCGGCGGCGGCATCCATCAGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGcc                                                                                                                                <  1:1626327/150‑1 (MQ=255)
                           ggTCACGTTGGTCGCTGTTGTAGGCCAGGACCACGTCCACCCGCCCGGCGGCGGCATCCATCAGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGc                                                                                                                            <  1:1071005/150‑1 (MQ=255)
                               aCGTTGGTCGCTGTTGTAGGCCAGGACCACGTCCACCCGCCCGGCGGCGGCATCCATCAGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCcgcg                                                                                                                        >  1:465725/1‑150 (MQ=255)
                               aCGTTGGTCGCTGTTGTAGGCCAGGACCACGTCCACCCGCCCGGCGGCGGCATCCATCAGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCcgcg                                                                                                                        >  1:613674/1‑150 (MQ=255)
                                   tGGTCGCTGTTGTAGGCCAGGACCACGTCCACCCGCCCGGCGGCGGCATCCATCAGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCAcc                                                                                                                    <  2:272257/150‑1 (MQ=255)
                                      tCGCTGTTGTAGGCCAGGACCACGTCCACCCGCCCGGCGGCGGCATCCATCAGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGa                                                                                                                 >  1:1576815/1‑150 (MQ=255)
                                       cGCTGTTGTAGGCCAGGACCACGTCCACCCGCCCGGCGGCGGCATCCATCAGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGAt                                                                                                                >  2:1034218/1‑150 (MQ=255)
                                         cTGTTGTAGGCCAGGACCACGTCCACCCGCCCGGCGGCGGCATCCATCAGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATgg                                                                                                              <  1:196211/150‑1 (MQ=255)
                                          tgttgtAGGCCAGGACCACGTCCACCCGCCCGGCGGCGGCATCCATCAGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGc                                                                                                             >  2:136243/1‑150 (MQ=255)
                                          tgttgtAGGCCAGGACCACGTCCACCCGCCCGGCGGCGGCATCAATCAGGCAGCGCTTCCAGTTGCCGAGCACCACGGTGTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGc                                                                                                             >  2:598280/1‑150 (MQ=255)
                                                  gCCAGGACCACGTCCACCCGCCCGGCGGCGGCATCCATCAGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCAcc                                                                                                                    >  2:1123108/1‑135 (MQ=255)
                                                  gCCAGGACCACGTCCACCCGCCCGGCGGCGGCATCCATCAGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCAcc                                                                                                                    <  1:1123108/135‑1 (MQ=255)
                                                  gCCAGGACCACGTCCACCCGCCCGGCGGCGGCATCCATCAGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTc                                                                                                     >  1:1405509/1‑150 (MQ=255)
                                                          cACGTCCACCCGCCCGGCGGCGGCATCCATCAGGCAGCGCTTCCAGTTGCCCAGCACTACGGTTTGCACCTCGTAGCCCATGCGCGCAAAAAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCAACTGATGGCAGTCGCTCCAGGAAAt                                                                                             >  1:77229/1‑150 (MQ=255)
                                                           aCGTCCACCCGCCCGGCGGCGGCATCCATCAGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATg                                                                                            >  1:625552/1‑150 (MQ=255)
                                                                cACCCGCCCGGCGGCGGCATCCATCAGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCgg                                                                                       >  2:807208/1‑150 (MQ=255)
                                                                 aCCCGCCCGGCGGCGGCATCCATCAGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGaaa                                                                                              >  1:592965/1‑142 (MQ=255)
                                                                 aCCCGCCCGGCGGCGGCATCCATCAGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGaaa                                                                                              <  2:592965/142‑1 (MQ=255)
                                                                   ccgcccgGCGGCGGCATCCATCAGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTa                                                                                    >  1:814699/1‑150 (MQ=255)
                                                                   ccgcccgGCGGCGGCATCCATCAGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTa                                                                                    >  2:297515/1‑150 (MQ=255)
                                                                   ccgcccgGCGGCGGCATCCATCAGGCAGCGCTTCCAGTGGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTa                                                                                    >  2:948032/1‑150 (MQ=255)
                                                                       ccgGCGGCGGCATCCATCAGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTACAcc                                                                                <  1:807208/150‑1 (MQ=255)
                                                                         ggcggcggcATCCATCAGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTACACCgg                                                                              <  2:613674/150‑1 (MQ=255)
                                                                         ggcggcggcATCCATCAGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTACACCgg                                                                              <  1:136243/150‑1 (MQ=255)
                                                                         ggcggcggcATCCATCAGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTACACCgg                                                                              <  1:600643/150‑1 (MQ=255)
                                                                         ggcggcggcATCCATCAGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGAGGCCATGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTACACCgg                                                                              <  1:611128/150‑1 (MQ=255)
                                                                          gcggcggcATCCATCAGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTACACCggg                                                                             <  2:1070009/150‑1 (MQ=255)
                                                                             gcggcATCCATCAGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTACACCGGGTAg                                                                          <  1:506533/150‑1 (MQ=255)
                                                                                gcATCCATCAGTCAGCGTTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTACACCGGGTAGTCg                                                                       <  2:77229/150‑1 (MQ=255)
                                                                                       tcAGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTACACCGGGTAGTCGCAGTAgc                                                                <  1:948032/150‑1 (MQ=255)
                                                                                         aGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCggg                                                                                      >  1:1002155/1‑126 (MQ=255)
                                                                                         aGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCggg                                                                                      <  2:1002155/126‑1 (MQ=255)
                                                                                           gCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTACACCGGGTa                                                                           >  2:1315991/1‑135 (MQ=255)
                                                                                           gCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTACACCGGGTa                                                                           <  1:1315991/135‑1 (MQ=255)
                                                                                             aGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCggg                                                                                      <  1:271052/122‑1 (MQ=255)
                                                                                             aGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCggg                                                                                      >  2:271052/1‑122 (MQ=255)
                                                                                                    ccAGTTGCCCAGCACTACGGTTTGCACCTCGTAGCCCATGCGCGCGAACCGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAATTTGGCGGGTACACCGGGTAGTCGCAGTAGCGCACCGGCTcact                                                   <  1:659081/150‑1 (MQ=255)
                                                                                                    ccAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTACACCGGGTAGTCGCAGTAGCGCACCGGCTcact                                                   <  1:990004/150‑1 (MQ=255)
                                                                                                     cAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTACACCGGGTAGTCGCAGTAGCGCACCGGCTcactg                                                  <  2:551782/150‑1 (MQ=255)
                                                                                                      aGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTACACCGGGTAGTCGCAGTAGCGCACCGGCTcactgg                                                 <  1:1179018/150‑1 (MQ=255)
                                                                                                      aGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTACACCGGGTAGTCGCAGTAGCGCACCGGCTcactgg                                                 <  1:1565131/150‑1 (MQ=255)
                                                                                                       gTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTACACCGGGTAGTCGCAGTAGCGCAc                                                            >  1:553963/1‑138 (MQ=255)
                                                                                                       gTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTACACCGGGTAGTCGCAGTAGCGCAc                                                            <  2:553963/138‑1 (MQ=255)
                                                                                                            ccAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTACACCGGGTAGTCGCAGTAGCGCACCGGCTcactggcactgg                                           >  2:52808/1‑150 (MQ=255)
                                                                                                             ccataccacGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTACACCGGGTAGTCGCAGTAGCGCACCGGCTCactggcactggc                                          >  1:315705/5‑150 (MQ=255)
                                                                                                               gCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTACACCGGGTAGTCGCAGTAGCGCACCGGCTCACTGGCACTGGCGt                                        <  2:1405509/150‑1 (MQ=255)
                                                                                                                  ccacGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGTCGTCGCTCCAGGAAATGGGCGGGTACACCGGGTAGTCGCAGTAGCGCACCGGCTCACTGGCACTGGCGTGCa                                     <  2:1194213/150‑1 (MQ=255)
                                                                                                                  ccacGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTACACCGGGTa                                                                           <  2:156526/112‑1 (MQ=255)
                                                                                                                  ccacGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTACACCGGGTa                                                                           >  1:156526/1‑112 (MQ=255)
                                                                                                                   cacGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTACACCGGGTAGTGATAGGAGCGCACCGGCTCACTGGCACTGGCGTGCAg                                    >  2:794183/1‑150 (MQ=255)
                                                                                                                    acGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTACACCGGGTAGTCGCAGTAGCGCACCGGCTCACTGGCACTCGCGTgattt                                   <  2:814699/150‑5 (MQ=255)
                                                                                                                        tttGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGTGCGGGTACACCGGGTAGTCGCAGTAGCGCACCGGCTCACTGGCACTGGCGTGCAGGGTgg                               <  1:1300064/150‑1 (MQ=255)
                                                                                                                          tGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTACACCGGGTAGTCGCAGTAGCGCACCGGCTCACTGGCACTGGCGTGCAGGGTGGGc                             <  2:465725/150‑1 (MQ=255)
                                                                                                                            cACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTACACCGGGTAGTCGCAGTAGCGCACCGGCTCACTGGCACTGGCGTGCAGGGTgg                               <  1:370011/146‑1 (MQ=255)
                                                                                                                            cACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTACACCGGGTAGTCGCAGTAGCGCACCGGCTCACTGGCACTGGCGTGCAGGGTgg                               >  2:370011/1‑146 (MQ=255)
                                                                                                                             aCCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTACACCGGGAAGTCGCAGTAGCGCACCGGCTCACTGGCACTGGCGTGCAGGGTGGGCagc                          >  2:5471/1‑148 (MQ=255)
                                                                                                                              ccTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTACACCGGGTAGTCGCAGTAGCGCACCGGCTCACTGGCACTGGCGTGCAGGGTGGGCAACg                         >  1:671124/1‑150 (MQ=255)
                                                                                                                                  gTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTACACCGGGTAGTCGCAGTAGCGCACCGGCTcactgg                                                 >  2:664666/1‑122 (MQ=255)
                                                                                                                                  gTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTACACCGGGTAGTCGCAGTAGCGCACCGGCTcactgg                                                 <  1:664666/122‑1 (MQ=255)
                                                                                                                                  gTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTACACCGGGTAGTCGCAGTAGCGCACCGGCTCACTGGCACTGGCGTGGAGGGTGGGCAACGTCaa                     <  2:498924/150‑1 (MQ=255)
                                                                                                                                   tAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTACACCGGGTAGTCGCAGTAGCGCACCGGCTCACTGGCACTGGCGTGCAGGGTGGGCAACGTCAAg                    <  1:598280/150‑1 (MQ=255)
                                                                                                                                     gCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTACACCGGGTAGTCGCAGTAGCGCACCGGCTCACTGGCACTGGCGTGCAGGGTGGGCAACGTCAAGGc                  >  1:991706/1‑150 (MQ=255)
                                                                                                                                       ccATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGAAAATGGGCGGGTACACCGGGTAGTCGCAGTAGCGCACCGGATCAACGGCACTGGCGTGCAGGGTGGGCAACGTCAAGGCCa                <  1:794183/150‑1 (MQ=255)
                                                                                                                                        cATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTACACCGGGTAGTCGCAGTAGCGCACCGGCTCACTGGCACTGGCGTGCAGGGTGGGCAACGTCAAGGCCAg               >  1:1216363/1‑150 (MQ=255)
                                                                                                                                                cgAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTACACCGGGTAGTCGCAGTAGAGCACCGGCTCACTGGCACTGGCGTGCAGGGTGGGCAACGTCAAGGCCAGCGCCAACt       <  2:625552/150‑1 (MQ=255)
                                                                                                                                                     aGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTACACCGGGTAGTCGCAGTAGCGCACCGGCTCACTGGCACTGGCGTGCAGGGTGGGCAACGTCAAGGCCAGCGCCAACTGCCAg  <  2:671124/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGGCTCCGTGGAAAAACGCATACGCTGGTCACGTTGGTCGCTGTTGTAGGCCAGGACCACGTCCACCCGCCCGGCGGCGGCATCCATCAGGCAGCGCTTCCAGTTGCCCAGCACCACGGTTTGCACCTCGTAGCCCATGCGCGCGAACAGTTCGCGCACCACGGTGGGGGCCAGGCCGCGCACCTGATGGCCGTCGCTCCAGGAAATGGGCGGGTACACCGGGTAGTCGCAGTAGCGCACCGGCTCACTGGCACTGGCGTGCAGGGTGGGCAACGTCAAGGCCAGCGCCAACTGCCAG  >  NC_002947/4256638‑4256936

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