Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,068,184 T→A 6.9% L116Q (CTG→CAG)  PP_0926 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,068,1840TA 6.9% 224.0 / 6.4 72L116Q (CTG→CAG) PP_0926membrane protein
Reads supporting (aligned to +/- strand):  ref base T (32/35);  new base A (2/3);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

GGTCAAGCCCGGTGTCATCCTGTATTTCATCGTCGCCACCCTGGTCGGTTTCGCCATTGCCTGGGGCATGGCCATCTGGCGCAGCCCGATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTGGTCATCCTCATGTACAACTCGCTGTCAGCGGTGGTACTGGCGGTGTACAGCCCGGACCTCAAGTCCGATCCGTGGAGCATC  >  NC_002947/1068035‑1068329
                                                                                                                                                     |                                                                                                                                                 
ggTCAAGCCCGGTGTCATCCTGTATTTCATCGTCGCCACCCTGGTCGGTTTCGCCATTGCCTGGGGCATGGCCATCTGGCGCAGCCCGATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGcct                                                                                                                                                   <  2:1204389/150‑1 (MQ=255)
           gtgtCATCCTGTAGTTCATCGTCGCCACCCTGGTCGGTTTCGCCATTGCCTGGGGCATGGCCATCTGGCGCAGCCCGATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGgccg                                                                                                                                        >  2:956317/1‑150 (MQ=255)
            tgtCATCCTGTATTTCATCGTCGCCACCCTGGTCGGTTTCGCCATTGCCTGGGGCATGGCCATCTGGCGCAGCCCGATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGgccgc                                                                                                                                       >  2:142915/1‑150 (MQ=255)
             gtCATCCTGTATTTCATCGTCGCCACCCTGGTCGGTTTCGCCATTGCCTGGGGCATGGCCATCTGGCGCAGCCCGATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCAGGGCCAGgccgat                                                                                                                                      <  2:58706/150‑3 (MQ=255)
                 tCCTGTATTTCATCGTCGCCACCCTGGTCGGTTTCGCCATTGCCTGGGGCATGGCCATCTGGCGCAGCCCGATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGcc                                                                                                                                  >  2:153242/1‑150 (MQ=255)
                  ccTGTATTTCATCGTCGCCACCCTGGTCGGTTTCGCCATTGCCTGGGGCATGGCCATCTGGCGCAGCCCGATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCt                                                                                                                                 >  1:775422/1‑150 (MQ=255)
                    tGTATTTCATCGTCGCCACCCTGGTCGGTTTCGCCATTGCCTGGGGCATGGCCATCTGGCGCAGCCCGATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCAGGGCCAGgccgat                                                                                                                                      >  1:477390/1‑141 (MQ=255)
                    tGTATTTCATCGTCGCCACCCTGGTCGGTTTCGCCATTGCCTGGGGCATGGCCATCTGGCGCAGCCCGATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCAGGGCCAGgccgat                                                                                                                                      <  2:477390/143‑3 (MQ=255)
                         ttCATCGTCGCCACCCTGGTCGGTTTCGCCATTGCCTGGGGCATGGCCATCTGGCGCAGCCCGATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGt                                                                                                                          <  1:345836/150‑1 (MQ=255)
                         ttCATCGTCGCCACCCTGGTCGGTTTCGCCATTGCCTGGGGCATGGCCATCTGGCGCAGCCCGATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCATGtcaggt                                                                                                                          <  2:777214/150‑6 (MQ=255)
                             tcgtcgCCACCCTGGTCGGTTTCGCCATTGCCTGGGGCATGGCCATCTGGCGCAGCCCGATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACt                                                                                                                      >  1:498541/1‑150 (MQ=255)
                               gtcgCCACCCTGGTCGGTTTCGCCATTGCCTGGGGCATGGCCATCTGGCGCAGCCCGATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTAt                                                                                                                    <  1:1049675/150‑1 (MQ=255)
                               gtcgCCACCCTGGTCGGTTTCGCCATTGCCTGGGGCATGGCCATCTGGCGCAGCCCGATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTAt                                                                                                                    <  1:153242/150‑1 (MQ=255)
                                  gCCACCCTGGTCGGTTTCGCCATTGCCTGGGGCATGGCCATCTGGCGCAGCCCGATGGCCGACCTGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGc                                                                                                                 >  1:33103/1‑150 (MQ=255)
                                            tCGGTTTCGCCATTGCCTGGGGCATGGCCATCTGGCGCAGCCCGATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTgg                                                                                                       <  1:376883/150‑1 (MQ=255)
                                               gTTTCGCCATTGCCTGGGGCATGGCCATCTGGCGCAGCCCGATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTggggg                                                                                                    <  2:775422/150‑1 (MQ=255)
                                                tttCGCCATTGCCTGGGGCATGGCCATCTGGCGCAGCCCGATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGc                                                                                                   >  2:72669/1‑150 (MQ=255)
                                                tttCGCCATTGCCTGGGGCATGGCCATCTGGCGCAGCCCGATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCATGGGGGc                                                                                                   <  1:343999/150‑1 (MQ=255)
                                                      cATTGCCTGGGGCATGGCCATCTGGCGCAGCCCGATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACt                                                                                             <  1:72669/150‑1 (MQ=255)
                                                           ccTGGGGCATGGCCATCTGGCGCAGCCCGATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAg                                                                                        <  1:896521/150‑1 (MQ=255)
                                                              ggggCATGGCCATCTGGCGCAGCCCGATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACt                                                                                                                      >  2:809595/1‑117 (MQ=255)
                                                              ggggCATGGCCATCTGGCGCAGCCCGATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACt                                                                                                                      <  1:809595/117‑1 (MQ=255)
                                                               gggCATGGCCATCTGGCGCAGCCCGATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCt                                                                                    <  1:113680/150‑1 (MQ=255)
                                                                ggCATGGCCATCTGGCGCAGCCCGATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTg                                                                                   <  1:956317/150‑1 (MQ=255)
                                                                 gCATGGCCATCTGGGGCAGCCCGATGGCCGACCGGGGTATCTATATCCAGGGCGCCTTCCGCGGCAACAACGGTGTTATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTgg                                                                                  <  2:362504/150‑1 (MQ=255)
                                                                 gCATGGCCATCTGGCGCAGCCCGATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTgg                                                                                  >  1:1128756/1‑150 (MQ=255)
                                                                 gCATGGCCATCTGGCGCAGCCCGATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCAGGGCCAGgccgat                                                                                                                                      <  2:111041/98‑3 (MQ=255)
                                                                 gCATGGCCATCTGGCGCAGCCCGATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCAGGGCCAGgccgat                                                                                                                                      >  1:111041/1‑96 (MQ=255)
                                                                     ggCCATCTGGCGCAGCCCGATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTg                                                                                   >  2:992905/1‑145 (MQ=255)
                                                                     ggCCATCTGGCGCAGCCCGATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTg                                                                                   <  1:992905/145‑1 (MQ=255)
                                                                      gCCATCTGGCGCAGCCCGATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTGGTCATc                                                                             >  1:1020491/1‑150 (MQ=255)
                                                                       ccATCTGGCGCAGCCCGATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTGGTCATcc                                                                            >  2:933957/1‑150 (MQ=255)
                                                                          tCTGGCGCAGCCCGATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTGGTCATCCTCa                                                                         >  2:399004/1‑150 (MQ=255)
                                                                                gcAGCCCGATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGcc                                                                                     <  2:350019/132‑1 (MQ=255)
                                                                                gcAGCCCGATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGcc                                                                                     >  1:350019/1‑132 (MQ=255)
                                                                                gcAGCCCGATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTGGTCATCCTCATGTACa                                                                   >  2:51540/1‑150 (MQ=255)
                                                                                  aGCCCGATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTGGTCATCCTCATGTACAAc                                                                 >  1:826700/1‑150 (MQ=255)
                                                                                      cGATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTGGTCATCCTCATGTACAACTCGc                                                             >  2:1150309/1‑150 (MQ=255)
                                                                                       gATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACACCGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTGGTCATCCTCATGTACAACTCGCt                                                            <  1:977971/150‑1 (MQ=255)
                                                                                       gATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTGGTCATCCTCATGTACAACTCGCt                                                            <  2:1128756/150‑1 (MQ=255)
                                                                                           gCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTGGTCATCCTCATGTACAACTCGCTGTCa                                                        <  2:498541/150‑1 (MQ=255)
                                                                                                 cGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGCGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTGGTCATCCTCATGTACAACTCGCTGTCAGCggta                                                  <  2:826700/150‑2 (MQ=255)
                                                                                                   gggCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATgg                                                                                                                  >  1:291051/1‑84 (MQ=255)
                                                                                                   gggCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATgg                                                                                                                  <  2:291051/84‑1 (MQ=255)
                                                                                                   gggCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTg                                                                                   <  2:375552/115‑1 (MQ=255)
                                                                                                   gggCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTg                                                                                   >  1:375552/1‑115 (MQ=255)
                                                                                                       aTCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTGGTCATCCTCATGTACAACTCGCTGTCAGCGGTGGTACTg                                            >  2:231512/1‑150 (MQ=255)
                                                                                                       aTCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTGGTCATCCTCATGTACAACTCGCTGTCAGCGGTGGTACTg                                            >  2:231539/1‑150 (MQ=255)
                                                                                                        tCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTGGTCATCCTCATGTACAACTCGCTGTCAGCGGTGGTACTgg                                           >  2:175619/1‑150 (MQ=255)
                                                                                                          tataCCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTGGTCATCCTCATGTACAACTCGCTGTCAGCGGTGGTACTGGCg                                         <  1:231512/150‑1 (MQ=255)
                                                                                                          tataCCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTGGTCATCCTCATGTACAACTCGCTGTCAGCGGTGGTACTGGCg                                         <  1:231539/150‑1 (MQ=255)
                                                                                                          tataCCCAGGGCGCCTTCCGCGGCAACAACGCTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTGGTCATCCTCATGTACAACTCGCTGTCAGCGGTGGTACTGGCg                                         >  2:743908/1‑150 (MQ=255)
                                                                                                           ataCCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTGGTCATCCTCATGTACAACTCGCTGTCAGCGGTGGTACTGGCgg                                        <  1:933957/150‑1 (MQ=255)
                                                                                                               ccAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTCCTCGCAGGCCTgg                                                                                  <  1:818340/104‑1 (MQ=255)
                                                                                                               ccAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTgg                                                                                  >  2:818340/1‑104 (MQ=255)
                                                                                                               ccAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTGGTCATCCTCATGTACAACTCGCTGTCAGCGGTGGTACTGGCGGTGTa                                    >  2:1178421/1‑150 (MQ=255)
                                                                                                                        cTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTGGTCATCCTCATGTACAACTCGCTGTCAGCGGTGGTACTGGCGGTGTACAGCCCGGa                           >  1:245654/1‑150 (MQ=255)
                                                                                                                             gcgGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTGGTCATCCTCATGTACAACTCGCTGTCAGCGGTGGTACTGGCGGTGTACAGCCCGGACCTCa                      >  2:656040/1‑150 (MQ=255)
                                                                                                                                gCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTGGTCATCCTCATGTACAACTCGCTGTCAGCGGTGGTACTGGCGGTGTACAGCCCGGACCTCAAGt                   <  2:245654/150‑1 (MQ=255)
                                                                                                                                  aacaacGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCg                                                                                           <  1:324829/76‑1 (MQ=255)
                                                                                                                                  aacaacGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCg                                                                                           >  2:324829/1‑76 (MQ=255)
                                                                                                                                    caacGGTGTCATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTGGTCATCCTCATGTACAACTCGCTGTCAGCGGTGGTACTGGCGGTGTACAGCCCGGACCTCAAGTCCGa               <  1:142915/150‑1 (MQ=255)
                                                                                                                                     aacGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTGGTCATCCTCATGTACAACTCGCTGTCAGCGGTGGTACTGGCGGTGTACAGccc                              >  1:596724/1‑134 (MQ=255)
                                                                                                                                     aacGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTGGTCATCCTCATGTACAACTCGCTGTCAGCGGTGGTACTGGCGGTGTACAGccc                              <  2:596724/134‑1 (MQ=255)
                                                                                                                                     aacGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTGGTCATCCTCATGTACAACTCGCTGTCAGCGGTGGTACTGGCGGTGTACAGCCCGGACCTCAAGTCCGAt              >  1:1046976/1‑150 (MQ=255)
                                                                                                                                     aacGGTGTGATCGGCCTGGCCCCGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTGGTCATCCTCATGTACAACTCGCTGTCAGCGGTGGTACTGGCGGTGTACAGCCCGGACCTCAAGTCCGAt              >  1:513556/1‑150 (MQ=255)
                                                                                                                                         gtgtGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTGGTCATCCTCATGTACAACTCGCTGTCAGCGGTGGTACTGGCGGTGTACAGCCCGGACCTCAAGTCCGATCCGt          <  2:165080/150‑1 (MQ=255)
                                                                                                                                         gtgtGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTGGTCATCCGCATGTACAACTCGCTGTCAGCGGTGGTACTGGCGGTGTACAGCCCGGAACTCAAGTCCGATCCGt          <  2:33103/150‑1 (MQ=255)
                                                                                                                                          tgtgATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTGGTCATCCTCATGTACAACTCGCTGTCAGCGGTGGTACTGGCGGTGTACAGCCCGGACCTCAAGTCCGATCCGTg         <  1:862376/150‑1 (MQ=255)
                                                                                                                                              aTCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTGGTCATCCTCATGTACAACTCGCTGTCAGCGGTGGTACTGGCGGTGTACAGCCCGGACCTCAAGTCCGATCCGTGGAGc     <  2:751828/150‑1 (MQ=255)
                                                                                                                                                 ggCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTGGTCATCCTCATGTACAACTCGCTGTCAGCGGTGGTACTGGCGGTGTACAGCCCGGACCTCAAGTCCGATCCGTGGAGCATc  <  2:1020491/150‑1 (MQ=255)
                                                                                                                                                  gCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTGGTCATCCTCATGTACAACTCGCTg                                                           >  1:249055/1‑92 (MQ=255)
                                                                                                                                                  gCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTGGTCATCCTCATGTACAACTCGCTg                                                           <  2:249055/92‑1 (MQ=255)
                                                                                                                                                     tggccctggccGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTGGTCATCCTCATGTACAACTCGCTGTCAGCGGTGGTACTGGCg                                         >  1:18850/1‑107 (MQ=255)
                                                                                                                                                     tggccctggccGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTGGTCATCCTCATGTACAACTCGCTGTCAGCGGTGGTACTGGCg                                         <  2:18850/107‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GGTCAAGCCCGGTGTCATCCTGTATTTCATCGTCGCCACCCTGGTCGGTTTCGCCATTGCCTGGGGCATGGCCATCTGGCGCAGCCCGATGGCCGACCGGGGCATCTATACCCAGGGCGCCTTCCGCGGCAACAACGGTGTGATCGGCCTGGCCCTGGCCGCCAGCCTGTACGGTGACTATGGCATTTCCCTGGGGGCGGTACTCGCAGGCCTGGTCATCCTCATGTACAACTCGCTGTCAGCGGTGGTACTGGCGGTGTACAGCCCGGACCTCAAGTCCGATCCGTGGAGCATC  >  NC_002947/1068035‑1068329

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