Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,414,637 C→T 6.9% A535V (GCC→GTC)  PP_3026 → recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,414,6370CT 6.9% 235.9 / 4.3 83A535V (GCC→GTC) PP_3026recombinase
Reads supporting (aligned to +/- strand):  ref base C (33/44);  new base T (3/3);  total (36/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGATCGATACTTTCGATCGCATCGTGATTTACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTCCTGTACCTCGACGCATACACGCCAGGCGGCCCATGCGAC  >  NC_002947/3414490‑3414785
                                                                                                                                                   |                                                                                                                                                    
tgatCGATACTTTCGATCGCATCGTGATTTACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGcggccg                                                                                                                                                    >  2:1160705/1‑150 (MQ=255)
tgatCGATACTTTCGATCGCATCGTGATTTACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGcggccg                                                                                                                                                    >  2:2378891/1‑150 (MQ=255)
 gatCGATACTTTCGATCGCATCGTGATTTACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGcggccgg                                                                                                                                                   >  2:484646/1‑150 (MQ=255)
   tCGATACTTTCGATCGCATCGTGATTTACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCTGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCg                                                                                                                                                 >  1:2309879/1‑150 (MQ=255)
     gATACTTTCGATCGCATCGTGATTTACATGCGTGGTTAAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGCTTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTc                                                                                                                                               >  1:194775/1‑150 (MQ=255)
      aTACTTTCGATCGCATCGTGATTTACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCg                                                                                                                                              >  2:1567598/1‑150 (MQ=255)
       tACTTTCGATCGCATCGTGATTTACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGa                                                                                                                                             <  2:1505293/150‑1 (MQ=255)
          tttCGATCGCATCGTGATTTACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCg                                                                                                                                          >  1:2190450/1‑150 (MQ=255)
                  gCATCGTGATTTACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGAGAGCGGAGATGGCTCAGGGTAGAAAGGAAGAGCGGAGAGTGGGGGGCCGGCGTCGACCGcccggcct                                                                                                                                  >  2:742068/1‑145 (MQ=255)
                    aTCGTGATTTACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGa                                                                                                                                <  1:1514484/150‑1 (MQ=255)
                     tCGTGATTTACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGAt                                                                                                                               <  1:1567598/150‑1 (MQ=255)
                     tCGTGATTTACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGAt                                                                                                                               <  1:134100/150‑1 (MQ=255)
                     tCGTGATTTACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGAt                                                                                                                               >  2:1691736/1‑150 (MQ=255)
                     tCGTGATTTACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGAAACAGGGTAGAAAGGAAGACCGGAGAGTGTGCGGCCGGCGTCGACCGCCCAGACTGAt                                                                                                                               >  1:24619/1‑150 (MQ=255)
                       gTGATTTACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGt                                                                                                                             <  2:2381834/150‑1 (MQ=255)
                          aTTTACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTa                                                                                                                          >  2:1065767/1‑150 (MQ=255)
                           tttACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTAc                                                                                                                         <  1:351461/150‑1 (MQ=255)
                           tttACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTAc                                                                                                                         <  1:2043791/150‑1 (MQ=255)
                           tttACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGCACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTAc                                                                                                                         <  2:1650925/150‑1 (MQ=255)
                           tttACATGCGGGGTGTAGTGCCAGCGGGGCGCAGGGCCAGCTACCTCGCTGTGTTGATGATTTCCAAGGCAGGACTGCGGCGGTGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTAc                                                                                                                         <  2:274178/150‑1 (MQ=255)
                            ttACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACg                                                                                                                        <  2:2043002/150‑1 (MQ=255)
                              aCATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCt                                                                                                                      <  1:1071861/150‑1 (MQ=255)
                                aTGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCCGGCACATGGGTGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGa                                                                                                                    <  2:194775/150‑1 (MQ=255)
                                aTGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGa                                                                                                                    <  2:2190450/150‑1 (MQ=255)
                                aTGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGa                                                                                                                    >  1:9504/1‑150 (MQ=255)
                                  gCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGAct                                                                                                                  <  2:2357719/150‑1 (MQ=255)
                                  gCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGAct                                                                                                                  >  1:2357719/1‑150 (MQ=255)
                                    gtggtgTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGActtc                                                                                                                <  1:1691736/150‑1 (MQ=255)
                                      ggtgTAGGGACAGCGGGGCGCAGGGCCAGGTACATCGAGGTGTTGCTGATTTCCAAGGAAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGTCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGActtctt                                                                                                              <  2:1036485/150‑1 (MQ=255)
                                        tgtAGTGCCAGCGGGGCGCAGCTCCACGGGCATCGATGTCTTGCTGGTTTCCAAGGCAGGACAGCGGGGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGTCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTgt                                                                                                            <  2:260838/150‑1 (MQ=255)
                                         gtAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTgtc                                                                                                           <  2:1861383/150‑1 (MQ=255)
                                          tAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTgtcg                                                                                                          <  1:827459/150‑1 (MQ=255)
                                           aGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTgtcgt                                                                                                         <  1:2378891/150‑1 (MQ=255)
                                            gTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTgtcgtc                                                                                                        >  1:221815/1‑150 (MQ=255)
                                                  gCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGAACGGCTGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTgg                                                                                                  >  1:1956759/1‑150 (MQ=255)
                                                   cGGAGCGCAGTGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTggg                                                                                                 <  1:1307206/150‑1 (MQ=255)
                                                    ggggCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGa                                                                                                <  1:2284228/150‑1 (MQ=255)
                                                         gcAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGccc                                                                                                                                       <  2:127846/106‑1 (MQ=255)
                                                         gcAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGccc                                                                                                                                       >  1:127846/1‑106 (MQ=255)
                                                         gcAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCg                                                                                           >  2:2225853/1‑150 (MQ=255)
                                                            gggCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGt                                                                                        <  2:221815/150‑1 (MQ=255)
                                                                cAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGCCTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGc                                                                                    <  1:1160705/150‑1 (MQ=255)
                                                                   gTACATCGATGTCTTGCTGATTTCCAGAAAATGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAg                                                                                 >  1:442542/1‑150 (MQ=255)
                                                                    tACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAgt                                                                                <  1:484646/150‑1 (MQ=255)
                                                                     aCATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGTCGACGCCGGCCGCCCa                                                                                                                                      <  2:364377/95‑1 (MQ=255)
                                                                     aCATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGTCGACGCCGGCCGCCCa                                                                                                                                      >  1:364377/1‑95 (MQ=255)
                                                                       aTCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGtgtg                                                                             <  2:2309879/150‑1 (MQ=255)
                                                                                  gCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGt                                                                  >  1:244208/1‑150 (MQ=255)
                                                                                   gtGATTGCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTg                                                                 <  2:2315544/149‑1 (MQ=255)
                                                                                   cTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTg                                                                 >  1:1175808/1‑150 (MQ=255)
                                                                                    tGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGGAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTgg                                                                >  1:1999820/1‑150 (MQ=255)
                                                                                     gATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTggg                                                                                                 <  2:938032/116‑1 (MQ=255)
                                                                                     gATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTggg                                                                                                 >  1:938032/1‑116 (MQ=255)
                                                                                     gATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGc                                                               >  2:1786783/1‑150 (MQ=255)
                                                                                      aTTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTggg                                                                                     <  2:1360483/127‑1 (MQ=255)
                                                                                      aTTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTggg                                                                                     >  1:1360483/1‑127 (MQ=255)
                                                                                      aTTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGcc                                                              >  1:297228/1‑150 (MQ=255)
                                                                                      aTTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTACTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGcc                                                              <  1:2069307/150‑1 (MQ=255)
                                                                                          ccAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGt                                                          <  2:1175808/150‑1 (MQ=255)
                                                                                           cAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGtc                                                         >  2:1226008/1‑150 (MQ=255)
                                                                                           cAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGtc                                                         <  1:2225853/150‑1 (MQ=255)
                                                                                             aGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGtccc                                                       >  2:1077397/1‑150 (MQ=255)
                                                                                              ggCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTAc                                                                                                                         <  1:843268/83‑1 (MQ=255)
                                                                                              ggCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTAc                                                                                                                         >  2:843268/1‑83 (MQ=255)
                                                                                                 aGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTggg                                                                                                 >  1:393213/1‑104 (MQ=255)
                                                                                                 aGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTggg                                                                                                 <  2:393213/104‑1 (MQ=255)
                                                                                                  ggACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCa                                                  <  1:1065767/150‑1 (MQ=255)
                                                                                                    aCAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGGGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGa                                                <  2:1770103/150‑1 (MQ=255)
                                                                                                       gcgccgCTGGCTCAGGGTACAAAGGAGGACCGGCGAGCGGGCGGCCGGCGTCGACCGCCCAGACTGAAGGTTACGGGGACTTCTTGTCGTCAATTGGGAGGCCGGGTGGGCTACTCCGCCCCACCAGGTGGCCTTAACCCTCACAGAAAt                                             >  2:762177/1‑150 (MQ=255)
                                                                                                        cggcgATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGTCGGCGTCGACCGCCCAGACTGATGTTTACGCTGActtc                                                                                                                >  2:2170443/1‑82 (MQ=255)
                                                                                                        cggcgATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGTCGGCGTCGACCGCCCAGACTGATGTTTACGCTGActtc                                                                                                                <  1:2170443/82‑1 (MQ=255)
                                                                                                           cgATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTc                                         <  2:244208/150‑1 (MQ=255)
                                                                                                             aTGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTCCt                                       <  2:9504/150‑1 (MQ=255)
                                                                                                                      ggTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTCCTGTACCTCGa                              >  1:810746/1‑150 (MQ=255)
                                                                                                                       gTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTCCTGTACCTCGAc                             <  2:1896835/150‑1 (MQ=255)
                                                                                                                           aaaGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGGTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTCCTGTACCTCGACGCAt                         <  1:344989/150‑1 (MQ=255)
                                                                                                                             aGGAAGACCGGAGAGTGGGCGGCCGTCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGAAGGCCGGGTGGGCTAGTGTGCACAACCAGGTGGCCTTGTCCCTCCCAGAAATTTTCCTGTACCTCGACGCATac                       >  2:864379/1‑150 (MQ=255)
                                                                                                                             aGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTgtc                                                                                                           <  1:2390951/66‑1 (MQ=255)
                                                                                                                             aGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTgtc                                                                                                           >  2:2390951/1‑66 (MQ=255)
                                                                                                                              ggAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTCCTGTACCTCGACGCATaca                      >  2:1026169/1‑150 (MQ=255)
                                                                                                                                     cGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCt                                                                                                                      <  1:662663/47‑1 (MQ=255)
                                                                                                                                     cGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCt                                                                                                                      >  2:662663/1‑47 (MQ=255)
                                                                                                                                        agtctGGGCGGTCGACGCCGGCCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAgt                                                                                <  1:74941/78‑1 (MQ=25)
                                                                                                                                        agtctGGGCGGTCGACGCCGGCCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAgt                                                                                >  2:74941/5‑82 (MQ=25)
                                                                                                                                            tGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCtt                                                            <  1:731748/98‑1 (MQ=255)
                                                                                                                                            tGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCtt                                                            >  2:731748/1‑98 (MQ=255)
                                                                                                                                               gcggcCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTCCTGTACCTCGACGCATACACGCCAGGCGGCCCATGc     >  2:355182/1‑150 (MQ=255)
                                                                                                                                                 ggccggcGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTCCTGTACCTCGACGCATACACGCCAGGCGGCCCATGCGa   <  1:1786783/150‑1 (MQ=255)
                                                                                                                                                  gccggcGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTCCTGTACCTCGACGCATACACGCCAGGCGGCCCATGCGAc  <  2:810746/150‑1 (MQ=255)
                                                                                                                                                  gccggcGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTCCTGTACCTCGACGCATACACGCCAGGCGGCCCATGCGAc  <  1:1026169/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGATCGATACTTTCGATCGCATCGTGATTTACATGCGTGGTGTAGTGCCAGCGGGGCGCAGGGCCAGGTACATCGATGTCTTGCTGATTTCCAAGGCAGGACAGCGGCGATGGCTCAGGGTAGAAAGGAAGACCGGAGAGTGGGCGGCCGGCGTCGACCGCCCAGACTGATGTTTACGCTGACTTCTTGTCGTCTATTGGGAGGCCGGGTGGGCTAGTGTGCACCACCAGGTGGCCTTGTCCCTCCCAGAAATTTTCCTGTACCTCGACGCATACACGCCAGGCGGCCCATGCGAC  >  NC_002947/3414490‑3414785

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: