Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,008,361 C→T 5.0% S259S (TCG→TCA PP_2628 ← ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,008,3610CT 5.0% 236.7 / 2.1 79S259S (TCG→TCAPP_2628ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (35/40);  new base T (2/2);  total (37/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 = 5.51e-01

GCTCACCAGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAGGGCAATGATCGCGATGAAAT  >  NC_002947/3008217‑3008509
                                                                                                                                                |                                                                                                                                                    
gCTCACCAGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCg                                                                                                                                                 >  2:55158/1‑150 (MQ=255)
 cTCACCAGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCgg                                                                                                                                                <  1:411651/150‑1 (MQ=255)
  tCACCAGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGCCGCCCgcg                                                                                                                                               >  1:1817028/1‑148 (MQ=255)
  tCACCAGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCggg                                                                                                                                               >  1:349578/1‑150 (MQ=255)
     ccAGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCGTCGATGACCGGGcgg                                                                                                                                            <  2:1038100/150‑3 (MQ=255)
      cAGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCAtt                                                                                                                                           >  1:1266574/1‑150 (MQ=255)
      cAGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCAtt                                                                                                                                           <  1:55158/150‑1 (MQ=255)
         gcgcATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGc                                                                                                                                        <  2:1492562/150‑1 (MQ=255)
         gcgcATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGc                                                                                                                                        <  2:1136796/150‑1 (MQ=255)
            cATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTg                                                                                                                                     <  2:1101824/150‑1 (MQ=255)
               cTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTg                                                                                                                                          <  2:277366/142‑1 (MQ=255)
               cTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTg                                                                                                                                          >  1:277366/1‑142 (MQ=255)
                 ggCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCt                                                                                                                                >  1:1350027/1‑150 (MQ=255)
                  gCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCATGCTCCTg                                                                                                                               >  2:851114/1‑150 (MQ=255)
                    gAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGcc                                                                                                                             >  2:622307/1‑150 (MQ=255)
                    gAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGCTCATCGACGACCGGGCATTGGCCTGCTCCTGcc                                                                                                                             >  1:456530/1‑150 (MQ=255)
                        aCCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGAc                                                                                                                         >  1:517186/1‑150 (MQ=255)
                             tCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGcg                                                                                                                                             >  1:1655693/1‑124 (MQ=255)
                             tCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGcg                                                                                                                                             <  2:1655693/125‑2 (MQ=255)
                             tCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGc                                                                                                                    >  1:387996/1‑150 (MQ=255)
                                 cTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAg                                                                                                                >  2:1049424/1‑150 (MQ=255)
                                   tGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGcc                                                                                                              >  1:129271/1‑150 (MQ=255)
                                      aaTCTTCCATGGCTGCGCACGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGAcc                                                                                                                        <  1:1083740/137‑1 (MQ=255)
                                      aaTCTTCCATGGCTGCGCACGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGAcc                                                                                                                        >  2:1083740/1‑137 (MQ=255)
                                              aTGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGaaa                                                                                                   <  2:129271/150‑1 (MQ=255)
                                              aTGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGaaa                                                                                                   <  2:517186/150‑1 (MQ=255)
                                               tGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGAcc                                                                                                                        <  2:1352289/128‑1 (MQ=255)
                                               tGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGAcc                                                                                                                        >  1:1352289/1‑128 (MQ=255)
                                                  cTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAATCCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACa                                                                                               <  1:22950/150‑1 (MQ=255)
                                                   tGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAg                                                                                              <  2:149673/150‑1 (MQ=255)
                                                   tGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAg                                                                                              <  2:483011/150‑1 (MQ=255)
                                                      gcgcgcAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCAt                                                                                           <  2:1151837/150‑1 (MQ=255)
                                                      ccgcgcACAGTGCGCTTCGCGGTTGGAGTGCGAAACCAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCCACCTTGCCCAGCCGGGGGACGAAATACCGCAt                                                                                           <  2:1155464/149‑1 (MQ=255)
                                                       cgctcAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTAATATCGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATg                                                                                          <  1:1180571/150‑1 (MQ=255)
                                                       cgcgcAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGCGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATg                                                                                          <  2:530175/150‑1 (MQ=255)
                                                        gcgcAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATgg                                                                                         <  1:851114/150‑1 (MQ=255)
                                                        gcgcAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATgg                                                                                         >  1:1535080/1‑150 (MQ=255)
                                                         cgcAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGCTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGc                                                                                        <  1:1666540/150‑1 (MQ=255)
                                                            aaaGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCaa                                                                                     <  2:402596/150‑1 (MQ=255)
                                                            aaaGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCaa                                                                                     >  2:1010997/1‑150 (MQ=255)
                                                              ggtggggttCGCGGTTGGAGTGGGAAAACAGGTTTGCCGTGGTGATGTTGGTGGAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGc                                                                                   <  2:1817028/143‑1 (MQ=255)
                                                                     ttCGCTGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCggg                                                                                                           >  2:1935248/1‑119 (MQ=255)
                                                                     ttCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCggg                                                                                                           <  1:1935248/119‑1 (MQ=255)
                                                                     ttCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACa                                                                                               <  1:797643/131‑1 (MQ=255)
                                                                     ttCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACa                                                                                               >  2:797643/1‑131 (MQ=255)
                                                                                 gTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGc                                                                                                                                        <  2:243261/78‑1 (MQ=255)
                                                                                 gTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGc                                                                                                                                        >  1:243261/1‑78 (MQ=255)
                                                                                    cGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAAt                                                             <  2:1535080/150‑1 (MQ=255)
                                                                                    cGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCAGAACAGCACGACCCAGGCAAt                                                             <  1:44766/150‑1 (MQ=255)
                                                                                     gAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAGACCCGCCCGGTCATCGACGACCGGGCAtt                                                                                                                                           >  2:1622860/1‑71 (MQ=255)
                                                                                     gAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAGACCCGCCCGGTCATCGACGACCGGGCAtt                                                                                                                                           <  1:1622860/71‑1 (MQ=255)
                                                                                         aCAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAgg                                                        <  2:349578/150‑1 (MQ=255)
                                                                                           aGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGgc                                                      <  1:1444667/150‑1 (MQ=255)
                                                                                             ttttACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGAcc                                                                     >  1:1716370/1‑133 (MQ=255)
                                                                                             ttttACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGAcc                                                                     <  2:1716370/133‑1 (MQ=255)
                                                                                                    gtggtgATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGt                                             >  1:329282/1‑150 (MQ=255)
                                                                                                     tggtgATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGtt                                            <  2:1350027/150‑1 (MQ=255)
                                                                                                        tgATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCa                                         >  1:1480315/1‑150 (MQ=255)
                                                                                                          aTGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAgg                                       <  2:329282/150‑1 (MQ=255)
                                                                                                            gTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAg                                        <  2:1347086/147‑1 (MQ=255)
                                                                                                            gTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAg                                        >  1:1347086/1‑147 (MQ=255)
                                                                                                              tGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGACGTGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCaa                                   <  1:622307/150‑1 (MQ=255)
                                                                                                                    aaGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGTGCAAAgccgc                             >  2:1584391/1‑150 (MQ=255)
                                                                                                                      gCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTCGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGccgccg                           >  1:1601704/1‑150 (MQ=255)
                                                                                                                           aGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTc                                                                                                                                  >  1:1350411/1‑42 (MQ=255)
                                                                                                                           aGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTc                                                                                                                                  <  2:1350411/42‑1 (MQ=255)
                                                                                                                             aaaCCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAggg                    >  1:1301031/1‑150 (MQ=255)
                                                                                                                               aCCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGc                                                                                             <  2:1582062/75‑1 (MQ=255)
                                                                                                                               aCCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGc                                                                                             >  1:1582062/1‑75 (MQ=255)
                                                                                                                                   gcccgGTCGTCGATGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAGGGCAATGa              >  1:919048/1‑150 (MQ=255)
                                                                                                                                   gcccgGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGttt                                           >  2:468385/1‑121 (MQ=255)
                                                                                                                                   gcccgGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGttt                                           <  1:468385/121‑1 (MQ=255)
                                                                                                                                    cccgGTCGTCGATGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAGGGCAATGAt             >  2:20924/1‑150 (MQ=255)
                                                                                                                                     ccgGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGaa                                                                               >  2:956562/1‑83 (MQ=255)
                                                                                                                                     ccgGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGaa                                                                               <  1:956562/83‑1 (MQ=255)
                                                                                                                                      cgGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGACAGGCCGGAAGACCCGTGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAGGGCAATGATcg           >  1:1491593/1‑150 (MQ=255)
                                                                                                                                          cgtCGATGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGCAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAGGGCAATGATCGCGAt       <  1:20924/148‑1 (MQ=255)
                                                                                                                                          cATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAGGGCAATGATCGCGAt       <  1:1049424/150‑1 (MQ=255)
                                                                                                                                           aTCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAGGGCAATGATCGCGATg      >  1:1354624/1‑150 (MQ=255)
                                                                                                                                              gacgacCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGc                                                                >  2:1270280/1‑89 (MQ=255)
                                                                                                                                              gacgacCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGc                                                                <  1:1270280/89‑1 (MQ=255)
                                                                                                                                              gacgacCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAGGGCAATGATCGCGATGaaa   <  2:260458/150‑1 (MQ=255)
                                                                                                                                              gacgacCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAGGGCAATGATCGCGATGaaa   >  1:260458/1‑150 (MQ=255)
                                                                                                                                               acgacCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAGGGCAATGATCGCGATGAAAt  <  2:1601704/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GCTCACCAGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGTCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGTTTCAGGGCAAAGCCGCCGGCCAGGGCAATGATCGCGATGAAAT  >  NC_002947/3008217‑3008509

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: