Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,826,147 G→A 5.6% A1696A (GCC→GCT)  pvdL ← non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,826,1470GA 5.6% 268.8 / 3.8 89A1696A (GCC→GCT) pvdLnon‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base G (38/46);  new base A (2/3);  total (40/49)
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.53e-01

TGCCGCCGACGTTGTACGCCGGGCTGTCCGGCTCCATCTGCCAGAGGAACCACATTCGCTGCTGCGAGTACGACAGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGG  >  NC_002947/4825999‑4826294
                                                                                                                                                     |                                                                                                                                                   
tGCCGCCGACGTTGTACGCCGGGCTGTCCGGCTCCATCTGCCAGAGGAACCACATTCGCTGCTGCGAGTACGACAGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGg                                                                                                                                                    <  1:1061626‑M2/150‑2 (MQ=255)
     ccgACGTTGTACGCCGGGCTGTCCGGCTCCATCTGCCAGAGGAACCACATTCGCTGCTGCGAGTCCGACAGCGGCACTGCCTGACGGCGGTCGACACAGGTGATTGCACCGTGCAGTTTGCGGTCGCCGGC‑CGCGCGCAGGCGGGCGAtc                                                                                                                                               <  2:711196/150‑1 (MQ=255)
        aCGTTGTACGCCGGGCTGTCCGGCTCCATCTGCCAGAGGAACCACATTCGCTGCTGCGAGTACGACAGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCg                                                                                                                                            <  2:315007/150‑1 (MQ=255)
                cGCCGGGCTGTCCGGCTCCATCTGCCAGAGGAACCACATTCGCTGCTGCGAGTACGACAGCGGCACTGCCTGACGGCGGTGGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGaa                                                                                                                                    <  1:688769/150‑1 (MQ=255)
                    gggCTGTCCGGCTCCATCTGCCAGAGGAACCACATTCGCTGCTGCGAGTACGACAGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGgcg                                                                                                                                >  1:281748/1‑150 (MQ=255)
                       cTGTCCGGCTCCATCTGCCAGAGGAACCACATTCGCTGCTGCGAGTACGACAGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACGTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGccc                                                                                                                             <  1:344522/150‑1 (MQ=255)
                             ggCTCCATCTGCCAGAGGAACCACATTCGCTGCTGCGAGTACGACAGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCg                                                                                                                       >  2:874511/1‑150 (MQ=255)
                               cTCCATCTGCCAGAGGAACCACATTCGCTGCTGCGAGTACGACAGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCt                                                                                                                     <  2:293721/150‑1 (MQ=255)
                                tCCATCTGCCAGAGGAACCACATTCGCTGCTGCGAGTACGACAGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTg                                                                                                                    <  2:261545/150‑1 (MQ=255)
                                tCCATCTGCCAGAGGAACCACATTCGCTGCTGCGAGTACGACAGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCAGTCTCGTCACAGAAGGCGCCCAGCTCGCTg                                                                                                                    <  1:543630/150‑1 (MQ=255)
                                        ccAGAGGAACCACATTCGCTGCTGCGAGTACGACAGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGaa                                                                                                            >  2:533335/1‑150 (MQ=255)
                                         cAGTGGAAACACATTCGCTGTTGCGAGTACGACAGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAAc                                                                                                           >  2:1044203/1‑150 (MQ=255)
                                           gagGAACCACATTCGCTGCTGCGAGTACGACAGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAg                                                                                                         <  1:994840/150‑1 (MQ=255)
                                           gagGAACCACATTCGCTGCTGCGAGTACGACAGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCCCGCCTGCGCGCG‑GCCGGCGAgcgcaacctgcaaggt                                                                                                                                     >  1:1037655‑M2/1‑106 (MQ=255)
                                           gagGAACCACATTCGCTGCTGCGAGTACGACAGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCCCGCCTGCGCGCG‑GCCGGCGAgcgcaacctgcaaggt                                                                                                                                     <  2:1037655‑M2/122‑17 (MQ=255)
                                            agGAACCACATTCGCTGCTGCGAGTACGACAGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAgg                                                                                                        <  2:1061626/150‑1 (MQ=255)
                                             ggAACCACATTCGCTGCTGCGAGTACGACAGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAggg                                                                                                       >  2:1449350/1‑150 (MQ=255)
                                               aaCCACATTCGCTGCTGCGAGTACGACAGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCCCGCCTGCGCGCG‑GCCGGCGAgcgcaacctgc                                                                                                                                          >  1:1770641‑M2/1‑102 (MQ=255)
                                               aaCCACATTCGCTGCTGCGAGTACGACAGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCCCGCCTGCGCGCG‑GCCGGCGAgcgcaacctgc                                                                                                                                          <  2:1770641‑M2/113‑12 (MQ=255)
                                                aCCACATTCGCTGCTGCGAGTACGACAGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCt                                                                                                    <  1:680336/150‑1 (MQ=255)
                                                 ccACATTCGCTGCTGCGAGTACGACAGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCtt                                                                                                   >  2:1235026/1‑150 (MQ=255)
                                                   acaTTCGCTGCTGCGAGTACGACAGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGgc                                                                                                                                 <  1:683572/118‑1 (MQ=255)
                                                   acaTTCGCTGCTGCGAGTACGACAGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGgc                                                                                                                                 >  2:683572/1‑118 (MQ=255)
                                                     aTTCGCTGCTGCGAGTACGACAGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGc                                                                                               <  1:386767/150‑1 (MQ=255)
                                                      ttCGCTGCTGCGAGTACGACAGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCg                                                                                              >  1:806464/1‑150 (MQ=255)
                                                      ttCGCTGCTGCGAGTACGACAGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCg                                                                                              >  2:1702404/1‑150 (MQ=255)
                                                          ctgctgCGAGTACGACAGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAg                                                                                          >  2:42385/1‑150 (MQ=255)
                                                          ctgctgCGAGTACGACAGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGAGCTCGGCACAGAAGGCGCCAAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAg                                                                                          >  2:1175856/1‑150 (MQ=255)
                                                               gcGAGTACGACAGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGcaca                                                                                                                                       >  2:1485867/1‑100 (MQ=255)
                                                               gcGAGTACGACAGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGcaca                                                                                                                                       <  1:1485867/100‑1 (MQ=255)
                                                                cGAGTACGACAGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGc                                                                                                                  <  1:109779/120‑1 (MQ=255)
                                                                cGAGTACGACAGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGc                                                                                                                  >  2:109779/1‑120 (MQ=255)
                                                                 gAGTACGACAGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACa                                                                                   <  1:533335/150‑1 (MQ=255)
                                                                  aGTACGACAGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTAGACAt                                                                                  >  1:455314/1‑150 (MQ=255)
                                                                       gACAGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAg                                                                             >  1:1442536/1‑150 (MQ=255)
                                                                        aCAGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAgg                                                                            >  1:1852965/1‑150 (MQ=255)
                                                                        aCAGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAg                                                                             >  1:1209385/1‑149 (MQ=255)
                                                                        aCAGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAg                                                                             <  2:1209385/149‑1 (MQ=255)
                                                                         cAGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGc                                                                           <  1:1931064/150‑1 (MQ=255)
                                                                          aGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGcc                                                                          >  2:1697694/1‑150 (MQ=255)
                                                                           gCGGCACTGCCTGACGTCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCt                                                                         >  1:521755/1‑150 (MQ=255)
                                                                             ggCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTgg                                                                       >  2:1387008/1‑150 (MQ=255)
                                                                              gCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGc                                                                      <  1:1145722/150‑1 (MQ=255)
                                                                                aCTGCCTGACGGCGGTCTACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCATCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCgg                                                                    <  1:1162621/150‑1 (MQ=255)
                                                                                 cTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCggg                                                                   <  1:1784648/150‑1 (MQ=255)
                                                                                   gCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTc                                                                                       <  2:1924178/128‑1 (MQ=255)
                                                                                   gCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTc                                                                                       >  1:1924178/1‑128 (MQ=255)
                                                                                   gCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTg                                                                 >  2:45465/1‑150 (MQ=255)
                                                                                    ccTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTgc                                                                <  1:45465/150‑1 (MQ=255)
                                                                                    ccTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTgc                                                                >  2:1245034/1‑150 (MQ=255)
                                                                                      tGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTgcgc                                                              <  2:1852965/150‑1 (MQ=255)
                                                                                      tGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTgcgc                                                              <  1:154057/150‑1 (MQ=255)
                                                                                          ggcggTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGaaa                                                          <  1:449095/150‑1 (MQ=255)
                                                                                          ggcggTCGACACGCGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCCGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTTAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGaaa                                                          <  1:1851046/150‑1 (MQ=255)
                                                                                             ggTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGa                                                       >  1:961260/1‑150 (MQ=255)
                                                                                              gTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGAt                                                      <  1:1702404/150‑1 (MQ=255)
                                                                                               tCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATc                                                     >  2:1157804/1‑150 (MQ=255)
                                                                                                cGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCt                                                    <  1:1044203/150‑1 (MQ=255)
                                                                                                  acacGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTgg                                                  <  2:521755/150‑1 (MQ=255)
                                                                                                   cacGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTggg                                                 >  1:1204436/1‑150 (MQ=255)
                                                                                                    acGGGTGATTGCACCTTGCAGGTTGCGCTCGCCAGC‑CGCGCTCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGt                                                >  2:560655/1‑150 (MQ=255)
                                                                                                      gggTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGaaa                                                          >  1:1404655/1‑138 (MQ=255)
                                                                                                      gggTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGaaa                                                          <  2:1404655/138‑1 (MQ=255)
                                                                                                      gggTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAg                                              >  1:525400/1‑150 (MQ=255)
                                                                                                       ggTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGc                                             >  1:1224321/1‑150 (MQ=255)
                                                                                                         tGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCa                                           <  1:1157804/150‑1 (MQ=255)
                                                                                                           aTTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGc                                         >  1:1872115/1‑150 (MQ=255)
                                                                                                            ttGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGCCTGGGGTAGGGGCCTTGAGCGGCAGCTCGACATAGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAgaa                                        >  1:337593/1‑148 (MQ=255)
                                                                                                              gCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAAc                                      <  2:525400/150‑1 (MQ=255)
                                                                                                                    tGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGa                                <  2:806464/150‑1 (MQ=255)
                                                                                                                     gCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCt                                                                                                    >  2:892437/1‑81 (MQ=255)
                                                                                                                     gCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCt                                                                                                    <  1:892437/81‑1 (MQ=255)
                                                                                                                             cgcTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAgg                       <  1:42385/150‑1 (MQ=255)
                                                                                                                              gcTCGCCGGC‑CGCGCGCAGGCGGgcgagcgcaacctgcaaggtgcaatcacccgtgtcgaccgccgtcaggcagtgccgctgtc                                                                                        >  1:3420‑M2/1‑23 (MQ=255)
                                                                                                                              gcTCGCCGGC‑CGCGCGCAGGCGGgcgagcgcaacctgcaaggtgcaatcacccgtgtcgaccgccgtcaggcagtgccgctgtc                                                                                        <  2:3420‑M2/84‑62 (MQ=255)
                                                                                                                              gcTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAATGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAggg                      >  1:146844/1‑150 (MQ=255)
                                                                                                                              gcTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCg                                                                     >  2:1450101/1‑103 (MQ=255)
                                                                                                                              gcTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCg                                                                     <  1:1450101/103‑1 (MQ=255)
                                                                                                                              gcTCGCCGGC‑CGCGCGCAGGCGAgcgcaacctgcaaggtgcaatcacccgtgtcgaccgccgtcaggcagtgcccctgtcgtactcgcagccgcgaatgtggttcctctggcagatggagccggacagcccggtgtacaacgtcggcggc                      <  2:157003‑M2/150‑128 (MQ=255)
                                                                                                                                tCGCCGGC‑CGCGCGCAGGCGGgcgagcgcaacctgcaaggtgcaatcacccgtgtcgaccgccgtcaggcagtgccgctgtcgtactcgcagcagcgaatgtggttcctctggcagatggagccggacagcccggcgtacaacgtcggcg                    <  1:571165‑M2/150‑130 (MQ=255)
                                                                                                                                 cGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGc                                                                      <  2:1362097/99‑1 (MQ=255)
                                                                                                                                 cGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGc                                                                      >  1:1362097/1‑99 (MQ=255)
                                                                                                                                     ggc‑cgCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCgaagaa               <  2:1872115/150‑1 (MQ=255)
                                                                                                                                     ggc‑cgCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCACTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCgaagaa               >  1:192845/1‑150 (MQ=255)
                                                                                                                                      gc‑cgCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGaagaag              <  1:1616047/150‑1 (MQ=255)
                                                                                                                                          gcgcgcAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAgtcg           <  1:874511/150‑1 (MQ=255)
                                                                                                                                            gcgcAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAgtcgtc         <  1:144183/150‑1 (MQ=255)
                                                                                                                                             cgcAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTcc        <  1:1235026/150‑1 (MQ=255)
                                                                                                                                             cgcAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTcc        <  2:1204436/150‑1 (MQ=255)
                                                                                                                                                 ggcgggcgATCTCGGCACAGAAGGCGCCCGGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCa    <  1:912570/150‑1 (MQ=255)
                                                                                                                                                 ggcgggcgATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAACCGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCa    <  1:48213/150‑1 (MQ=255)
                                                                                                                                                  gcgggcgATCTCGGCACAGAGGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAg                                                                             <  1:1338558/76‑1 (MQ=255)
                                                                                                                                                  gcgggcgATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAg                                                                             >  2:1338558/1‑76 (MQ=255)
                                                                                                                                                   cgggcgATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACAGAGGGCGAATAATTCTGCCTGCAgg  <  2:1442536/150‑1 (MQ=255)
                                                                                                                                                     ggcgATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCCACGAGTGACCACCGAgg                       >  2:902122/1‑127 (MQ=255)
                                                                                                                                                     ggcgATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAgg                       <  1:902122/127‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TGCCGCCGACGTTGTACGCCGGGCTGTCCGGCTCCATCTGCCAGAGGAACCACATTCGCTGCTGCGAGTACGACAGCGGCACTGCCTGACGGCGGTCGACACGGGTGATTGCACCTTGCAGGTTGCGCTCGCCGGC‑CGCGCGCAGGCGGGCGATCTCGGCACAGAAGGCGCCCAGCTCGCTGGCTTCGAACAGGGCCTTGAGCGGCAGCTCGACATCGCAGGCCTGGCGGGTGCGCGAAACGATCTGGGTAGCCAGCAACGAGTGACCACCGAGGGCGAAGAAGTCGTCCTGCAGG  >  NC_002947/4825999‑4826294

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